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Abstract 

Present-day  expeditionary  oceanography  is  beginning  to  shift  from  a  focus  on  short¬ 
term  ship  and  submersible  deployments  to  an  ocean  observatory  mode  where  long¬ 
term  temporally-focused  studies  are  feasible.  As  a  result,  a  critical  need  for  in  situ 
chemical  sensors  is  evolving.  New  sensors  take  a  significant  amount  of  time  to  develop; 
thus,  the  evaluation  of  techniques  in  the  laboratory  for  use  in  the  ocean  environment 
is  becoming  increasingly  important.  Laser-induced  breakdown  spectroscopy  (LIBS) 
possesses  many  of  the  characteristics  required  for  such  in  situ  chemical  sensing,  and 
is  a  promising  technique  for  field  measurements  in  extreme  environments.  Although 
many  LIBS  researchers  have  focused  their  work  on  liquid  jets  or  surfaces,  little  at¬ 
tention  has  been  paid  to  bulk  liquid  analysis,  and  especially  to  the  effect,  of  oceanic 
pressures  on  LIBS  signals.  In  this  work,  laboratory  experiments  validate  the  LIBS 
technique  in  a  simulated  deep  ocean  environment  to  pressures  up  to  2.76  x  10'  Pa.  A 
key  focus  of  this  work  is  the  validation  that  select  elements  important  for  understand¬ 
ing  hydrothermal  vent  fluid  chemistry  (Na,  Ca,  Mn,  Mg,  K,  and  Li)  are  detectable 
using  LIBS.  A  data  processing  scheme  that  accurately  deals  with  the  extreme  nature 
of  laser-induced  plasma  formation  was  developed  that  allows  for  statistically  accu¬ 
rate  comparisons  of  spectra.  The  use  of  both  single  and  double  pulse  LIBS  for  high 
pressure  bulk  aqueous  solutions  is  explored  and  the  system  parameters  needed  for 
the  detection  of  the  key  analytes  are  optimized.  Using  both  single  and  double  pulse 
LIBS,  the  limits  of  detection  were  found  to  be  higher  than  expected  as  a  result  of  the 
spectrometer  used  in  this  experimentation.  However,  the  results  of  this  validation 
show  that  LIBS  possesses  the  characteristics  to  be  a  viable  chemical  sensing  method 
for  in  situ  analyte  detection  in  high  pressure  environments  like  the  deep  ocean. 
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Chapter  1 
Introduction 


The  development  of  in  situ  chemical  sensors  is  needed  to  explore  and  understand 
the  ocean  and  its  processes.  In  present  day  oceanographic  work,  new  sensors  are 
required  for  expeditionary  science  with  underwater  vehicles  such  as  ROVs,  AUVs, 
and  submersibles.  A  new  paradigm  for  ocean  study  has  begun  with  the  implemen¬ 
tation  of  ocean  observatories.  As  these  permanent  observatories  become  the  new 
mode  of  oceanography,  there  will  be  a  critical  need  for  chemical  sensors  capable  of 
long-term  deployment  for  ocean  observatories  to  reach  their  full  potential.  The  time 
required  to  transform  a  bench-top  laboratory  technique  into  a  full  ocean-going  system 
is  significant.  The  development,  phase  initially  requires  validation  that  an  analytical 
technique  will  work  under  in  situ  conditions.  In  this  thesis,  laser-induced  breakdown 
spectroscopy  (LIBS)  is  evaluated  in  the  laboratory  to  determine  if  this  bench-top  an¬ 
alytical  technique  is  viable  for  development  into  a  field-going  oceanographic  chemical 
sensor. 


1.1  Laser-Induced  Breakdown  Spectroscopy 

Laser-induced  breakdown  spectroscopy  is  a  type  of  atomic  emission  spectroscopy  that 
was  first  reported  in  the  literature  in  1962  and  has  since  evolved  into  a  technique  for 
laboratory  chemical  analysis  [1].  LIBS  technology  is  currently  undergoing  transfor¬ 
mation  from  a  benchtop  analytical  technique  into  a  viable  tool  for  field  measurements, 
and  is  emerging  as  a  tool  for  chemical,  geochemical,  and  environmental  analysis  in 
extreme  and  hostile  environments  [2].  LIBS  has  been  used  in  the  analysis  of  a  wide 
variety  of  sample  types  including  soils  [3],  archaeological  artifacts  [4],  metal  alloys 
[5],  bacterial  spores,  molds,  pollens,  and  proteins  [6],  pharmaceuticals  [7],  glass  [8], 
nuclear  power  station  steam  generator  tubes  [9],  and  pigments  in  artwork  [10].  A 
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mobile  instrument  has  been  designed  for  studying  polluted  soils  [11,  12].  LIBS  is  also 
capable  of  stand-off,  non-contact  measurements,  and  a  field  deployable  system  has 
been  proposed  by  Palanco  et  al.  with  a  detection  range  on  the  order  of  hundreds  of 
meters  [13].  A  field-portable  LIBS  system  has  been  developed  for  landmine  detection 
[14].  Along  with  land-based  environmental  applications,  LIBS  is  finding  applicability 
to  space  exploration  [15  20].  LIBS  will  be  part  of  the  ChemCam  instrument  package 
for  the  Mars  Science  Laboratory  Rover  that  will  be  launched  in  2009  for  geological 
analysis. 

Although  there  are  many  proposed  land  and  space  based  applications  of  LIBS 
field-going  sensors,  it  also  holds  promise  for  in  situ  ocean  use.  There  are  numerous 
possible  applications  in  the  areas  of  chemical,  geological,  and  biological  oceanography 
ranging  from  laboratory  experimentation  to  fieldwork  in  environments  from  estuaries 
to  the  coastal  zone  to  the  deep  ocean.  The  development  of  an  oceanographic  LIBS 
sensor  could  allow  scientists  to  determine  the  chemical  composition  of  sediments, 
rocks,  or  ocean  fluids  in  an  in  situ,  real-time  mode.  Several  researchers  have  already 
applied  the  LIBS  technology  to  marine-related  applications.  For  example,  Niu  ct  al. 
successfully  used  LIBS  to  determine  strontium  levels  in  marine  algae  [21]  and  Barbini 
et  al.  used  LIBS  shipboard  to  analyze  marine  sediments  [22].  In  the  laboratory,  Do 
Giacomo  et  al.  have  explored  the  use  of  LIBS  for  the  detection  of  submerged  solid 
targets  [23  26]. 

1.1.1  How  LIBS  Works 

The  LIBS  technique  is  based  on  the  analysis  of  the  spectral  emission  from  laser- 
induced  plasmas,  and  is  a  type  of  atomic  emission  spectroscopy.  Atomic  emission 
spectroscopy  uses  an  external  energy  source  to  excite  ground  state  atoms.  The  atoms 
spontaneously  emit  radiation  as  they  revert  back  to  a  lower  energy  state,  with  the 
emission  intensity  being  proportional  to  the  concentration  of  atoms  in  the  ground 
state  [27].  In  the  LIBS  technique,  one  or  two  high-power,  pulsed  lasers,  typically 
Q-switched  NdrYAG  units,  are  used  as  the  excitation  source.  The  lasers  are  usually 
operated  at  the  fundamental  wavelength  of  1064  nm,  although  other  wavelengths  can 
be  used;  for  example,  the  lasers  can  be  frequency  doubled  to  deliver  532  nm.  Addi¬ 
tional  components  of  a  typical  LIBS  system  include  focusing  and  collection  optics,  a 
spectrometer,  and  a  data  acquisition  computer.  Commercially  available  broadband 
spectrometers  span  the  200  -  1000  nm  wavelength  region,  allowing  for  the  simultane¬ 
ous  detection  of  multiple  elements  [28]. 


20 


To  generate  a  plasma  spark,  a  short  duration,  high  power  laser  beam  is  focused 
onto  or  into  a  sample.  An  optically-induced  plasma  or  spark  is  formed  on  the  surface 
(of  a  solid  or  on  a  liquid)  or  in  the  sample  (in  bulk  liquid  or  in  a  gaseous  medium) 
when  the  laser  power  density  or  irradiance  exceeds  the  breakdown  threshold  of  the 
sample  [27].  Breakdown  threshold  irradiances  using  a  7-ns  pulse  width,  1064  nm 
Nd:YAG  laser  have  been  reported  by  Kennedy  et  al.  to  be  5.6  x  109  W/cm2  for 
tap  water  and  8.31  x  109  W/cm2  for  saline  solutions  [29].  The  plasma  radiates  both 
a  continuum  component  due  to  inverse  Bremsstrahlung  radiation  from  electron-ion 
collisions,  which  decays  rapidly,  and  an  emission  line  component  that  decays  more 
slowly.  Therefore,  the  plasma  emission  can  be  analyzed  by  spectroscopic  methods  by 
time  gating.  The  plasma  light  is  initially  dominated  by  a  white  light  continuum  which 
contains  little  intensity  variation  as  a  function  of  wavelength.  After  breakdown,  the 
plasma  expands  outwards  and  back  towards  the  focusing  lens.  The  expansion  occurs 
at  105  m/s  and  creates  an  audible  shock  wave  [30].  Plasma  temperatures  in  bulk 
liquids  are  in  the  range  of  7,000  -  12,000  K  [30].  Plasma  decay  occurs  by  radiative, 
quenching,  and  electron-ion  recombination  processes  that,  result  in  the  formation  of 
neutral  species  [31]. 

A  gated  spectrometer  covering  part  or  all  of  the  ultraviolet  through  near  infrared 
range  (200  nm  -  1000  nm)  is  used  to  capture  the  plasma  spectrum.  For  simultane¬ 
ous  multi-element,  analysis,  an  Echelle  spectrometer  can  be  used  which  contains  an 
Echelle  diffraction  grating  with  coarse  grooves  and  large  blaze  angles.  The  grooves 
have  steep  sides  to  cover  the  full  range  of  wavelengths,  and  a  prism  is  needed  to 
separate  overlapping  orders  of  the  grating  [27,  32].  CCDs  (charge  coupled  devices)  or 
ICCDs  (intensified  charge  coupled  devices)  serve  as  the  detector  devices  within  the 
spectrometers.  The  spectral  line  wavelengths  and  intensities  obtained  from  plasma 
ablation  can  be  compared  with  a  standard  atomic  line  reference  and/or  calibrated 
against,  samples  of  known  makeup  to  determine  the  elemental  composition  of  the 
sample.  The  intensity  of  the  spectral  lines  provides  a  quantitative  description  of  the 
elemental  concentration  [28,  33]. 


1.1.2  Advantages  and  Disadvantages  of  the  LIBS  Technique 

LIBS  as  a  spectrochemical  technique  possesses  many  advantages  which  make  it  es¬ 
pecially  attractive  for  development,  into  an  in  situ  sensor  for  oceanographic  research. 
However,  the  technique  also  possesses  several  disadvantages  that,  must,  be  addressed. 
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LIBS  can  provide  rapid  multi-element  analysis  and  has  the  capability  to  detect 
elements  in  the  parts-per-million  (ppm)  range.  All  elements  emit  in  the  200  -  940  nm 
wavelength  region  and  every  element  has  its  own  unique  spectral  signature;  therefore, 
all  elements  can  theoretically  be  detected  with  LIBS.  The  emission  lines  obtained 
from  plasma  ablation  can  be  compared  with  a  standard  reference  of  atomic  lines  to 
determine  the  chemical  composition  of  the  sample  by  comparison  of  wavelengths.  For 
example,  the  National  Institute  of  Standards  and  Technology  (NIST)  maintains  an  on¬ 
line  atomic  spectrum  database  (http:/ /physics. nist.gov/PhysRefData/ASD/index.litml). 
Unlike  other  techniques  that  are  useful  for  analyzing  one  form  of  material,  LIBS  can 
chemically  analyze  all  three  states  of  matter  (solid,  liquid,  and  gas).  The  sample  size 
for  LIBS  is  minimal  with  typically,  0.1  jig  -  0.1  mg  of  material  ablated  if  a  solid  sample 
is  analyzed.  The  sample  size  required  minimizes  the  destructiveness  and  invasiveness 
of  the  technique.  LIBS  does  not  require  any  preparation  of  a  sample  before  analysis 
and  unlike  many  analytical  techniques,  does  not  require  chemical  reagents  to  dissolve 
the  sample  [27].  This  avoids  contamination  problems  and  reduces  the  time  for  anal¬ 
ysis.  These  characteristics  makes  LIBS  suitable  for  real-time,  in  situ  analysis.  Most 
analytical  methods  (e.g.,  wet  chemistry  techniques)  cannot  produce  instant  results  in 
a  Held  environment.  Other  methods  can  be  time-intensive,  with  several  days  needed 
before  results  are  available.  Rapid  analysis  is  possible  with  LIBS  as  it  only  requires 
one  laser  pulse  to  generate  a  plasma,  although  typically  numerous  accumulat  ions  of 
plasmas  are  obtained.  A  large  number  of  measurements  can  be  taken  quickly,  making 
the  elemental  composition  of  the  sample  identifiable  on  a  nearly  real-time  basis. 

LIBS  has  been  identified  as  a  viable  technique  for  use  in  extreme  environments 
because  direct  contact  with  the  sample  is  not  necessary:  only  optical  access  is  required. 
Fiber  optics  can  be  used  to  reach  distances  far  from  the  laser.  The  laser  can  also  be 
focused  on  a  sample  at  a  distance,  making  LIBS  useful  for  remote  measurements.  The 
stand-off  analysis  capabilities  of  LIBS  make  it  a  viable  in  situ  analytical  technique. 
The  components  of  the  system  can  also  be  made  small  and  rugged.  Unlike  many 
traditional  analytical  techniques  that  require  collection  of  a  sample  followed  by  the 
transport  of  the  sample  to  a  laboratory,  LIBS  measurements  can  be  made  directly 
in  the  field.  This  eliminates  the  need  to  store  the  sample  in  conditions  that  differ 
from  the  natural  environment.  This  is  an  important  advantage  for  oceanography  due 
to  pressure  and  temperature  effects.  LIBS  can  also  be  used  for  surface  cleaning  of 
a  sample.  Repetitive  sampling  at  the  same  spot  can  be  used  for  ablation  through 
weathered  surfaces  to  reach  underlying  material.  This  could  be  useful  for  rocks  which 
could  be  “drilled"  through  to  areas  beneath  fouling  and  to  determine  chemical  makeup 
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of  the  fouling.  These  benefits  make  LIBS  useful  for  chemical  analysis  in  extreme 
environments  suggesting  that  it  is  a  viable  technique  for  use  in  the  deep  ocean  [28, 
31,  34-36]. 

Although  there  are  many  advantages  to  the  LIBS  technique,  there  are  several 
disadvantages  that  must  be  acknowledged.  LIBS  is  a  semi-quantitative  technique  as 
it  is  often  difficult  to  obtain  suitable  standards.  LIBS  does  not  have  the  sensitivity 
and  precision  of  many  state-of-the-art  laboratory  analytical  techniques.  Furthermore, 
measurements  are  highly  variable.  There  are  also  chemical  interference  (matrix)  ef¬ 
fects  between  analytes  that  must  be  addressed  [27]. 

1.1.3  LIBS  in  Liquids 

Several  researchers  have  been  successful  in  causing  plasma  ablation  to  occur  on  ma¬ 
terials  submerged  in  water  and  on  liquid  surfaces.  However,  only  a  few  researchers 
have  focused  their  work  on  LIBS  of  dissolved  analytes  within  bulk  aqueous  solut  ions 
[23,  37-42]  due  to  the  inherent  difficulty  of  detection.  If  the  liquid  is  transparent  at 
the  laser  wavelength,  a  plasma  can  be  formed  in  the  bulk  liquid  below  the  surface 
[37].  The  plasma  formed  in  a  bulk  liquid  compared  to  that  formed  on  a  solid  or  in  air 
displays  reduced  light  intensity  and  emission  lifetime  due  to  quenching  [37,  38,  43,  44]. 
This  emission  time  interval  is  on  the  order  of  1  /rs  or  less  in  liquid,  which  is  signifi¬ 
cantly  lower  than  at  the  air-liquid  interface  where  the  time  interval  averages  5  -  20 
/ts.  In  addition,  spectral  lines  are  broadened  by  the  Stark  effect  which  results  from 
ion  and  electron  collisions  [37,  45].  Furthermore,  “moving  breakdown”  occurs  that 
randomly  changes  the  distance  between  the  plasma  and  the  collection  optics,  a  phe¬ 
nomenon  that  is  not  relevant  to  solids  in  air.  The  plasma  expands  along  the  beam 
path,  resulting  in  an  elongated  plasma  that  cavitates  cylindrically  [46].  For  many 
aqueous  applications,  these  issues  can  be  avoided  by  analysis  on  a  liquid  surface,  jet, 
or  film;  however,  for  the  development  of  an  in  situ  oceanic  system,  it  is  necessary  to 
work  directly  with  bulk  liquids. 

As  early  as  1984,  Cremers  et  al.  [37]  were  able  to  identify  many  elements  in 
an  aqueous  solution.  This  work  showed  that  the  use  of  two  laser  sparks,  (double 
or  dual  pulse  LIBS),  with  a  short  interpulse  delay,  improved  the  detection  limit  for 
dissolved  species.  Using  two  laser  pulses  allowed  lower  concentrations  of  elements  to 
be  detected.  The  first  spark  creates  a  cavitation  bubble.  The  second  pulse  is  fired 
into  the  bubble,  forming  a  plasma  within  it,  allowing  the  measurement  to  be  made 
in  a  gaseous  environment. 
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Few  LIBS  researchers  have  focused  on  liquids  and  no  work  has  been  done  at  high 
pressures,  and  therefore  little  is  known  about  the  pressure  effect.  Relatively  few 
studies  have  focused  on  dissolved  species  in  liquids  and  therefore  much  work  on  t  he 
experimental  parameters  for  measurements  of  such  species  is  necessary. 


1.2  Scientific  Application:  Hydrothermal  Vent 
Chemistry 

Study  of  in  situ  hydrothermal  vent  chemistry  could  benefit  greatly  from  the  develop¬ 
ment  of  a  LIBS  sensor.  At  mid-ocean  ridges,  seawater  circulates  through  the  fractured 
and  permeable  crust.  Interactions  with  the  surrounding  rock  induces  major  chemical 
changes  to  the  fluid.  At  vent  orifices,  exit  temperatures  reach  200  -  405°C  at  ambient 
pressures  of  8.1  x  10°  Pa  to  3.6  x  10'  Pa  corresponding  to  ocean  depths  of  800  m 
to  3600  m  [47].  The  circulation  is  driven  by  the  direct  or  indirect  thermal  effects  of 
magma  at  sub-seafloor  depths  of  up  to  a  few  km;  additionally,  substantial  changes 
in  fluid  composition  occur  due  to  interaction  with  the  host  rock,  phase  separation, 
and  possibly  magma  degassing  [48].  For  example,  many  alkalis  (e.g.,  Li.  Na,  and  Ca) 
and  transition  metals  (e.g..  Fe,  Mn,  Cu,  and  Zn)  are  leached  from  the  host  rock  and 
concentrated  to  varying  degrees  in  the  fluid,  while  Mg  and  SOj  are  largely  removed 
from  the  fluid  [47]. 

As  vent  fluids  exit  the  seafloor,  the  very  hot  fluids  mix  with  ambient  seawater 
which  rapidly  changes  the  chemistry  (Figure  1-1).  Measuring  the  fluid  properties 
in  situ  is  very  difficult.  Collecting  samples  for  measurement  shipboard  or  back  in  a 
shore  laboratory  is  usually  done  by  using  non-reactive  titanium  samplers  to  extract 
water,  which  is  then  brought  to  the  surface.  Some  elements  remain  in  solution  as  the 
temperature  and  pressure  changes,  however,  others  precipitate  out  [49].  In  addition, 
the  chemistry  of  some  vents  are  known  to  change  over  short  (days  to  years)  time 
scales  [50].  The  use  of  a  method  like  LIBS,  therefore,  is  attractive  for  obtaining  an 
understanding  of  the  chemistry  of  vents  that  has  thus  far  been  impossible  to  achieve. 

Six  elements  (sodium,  calcium,  manganese,  magnesium,  potassium,  and  lithium) 
were  selected  as  the  primary  focus  of  this  work  because  of  their  key  importance  at 
hydrothermal  vents.  Sodium  (Na)  is  the  most  abundant  cation  in  vent  fluids  and  its 
study  is  important  for  gaining  an  understanding  of  phase  separation  processes  [48]. 
Manganese  (Mn)  exists  as  a  trace  metal  in  seawater  but  has  a  higher  concentration 
in  vent  fluids  due  to  leaching  from  the  host  rock  [48].  Mn  can  also  be  studied  si- 
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multaneously  with  Fe  as  an  indication  of  subsurface  deposition  as  Fe  precipitates  out 
while  Mn  stays  in  solution.  Calcium  (Ca)  is  the  second  most  abundant  cation  in 
vent  fluids,  and  is  typically  enriched  in  vent  fluids  compared  to  seawater  [50].  Ca  is 
released  into  vent  fluids  when  sodium  is  taken  up  during  albitization  reactions  with 
the  host  rock  [50].  Magnesium  (Mg)  is  very  low  to  nonexistent  in  hydrothermal  vent 
fluids;  however,  it  is  present  in  seawater  [48].  If  Mg  is  detected  in  vent  fluid  sam¬ 
ples,  it  is  evidence  for  contamination;  thus,  a  sensor  that  can  detect  Mg  is  desirable. 
Potassium  (K)  and  Lithium  (Li)  are  typically  highly  enriched  in  vent  fluids  due  to 
leaching  from  basalts  [48].  In  vent  fluids,  concentrations  range  from  approximately 
250  -  23,163  ppm  for  Na,  0.6  -  399  ppm  for  Mn,  -54  -  4477  ppm  for  Ca,  -47  -  3166 
ppm  for  K,  0  ppm  for  Mg.  and  0.7  -  1073  ppm  for  Li  [47].  In  seawater,  concentrations 
are  approximately  10933  ppm  Na,  <0.001  Mn,  419  ppm  Ca,  405  ppm  K,  1300  ppm 
Mg,  and  0.2  ppm  Li  [47]. 
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Figure  1-1:  At  hydrothermal  vents,  the  cold  seawater  seeps  down  through  the  perme¬ 
able  seafloor  where  it  undergoes  chemical  changes.  The  hot  vent  fluid  finally  vents  at 
the  seafloor.  Illustration  by  E.  Paul  Oberlander.  Reprinted  from  Oceanus,  Dec.  1, 
1998  with  permission  from  Susan  Humphris,  WHOI. 
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1.3  Thesis  Overview 


New  sensors  take  a  significant  amount  of  time  to  develop;  thus,  the  evaluation  of 
techniques  in  the  laboratory  for  use  in  the  ocean  environment  is  becoming  increas¬ 
ingly  important.  This  thesis  focuses  on  this  proof-of-concept  phase,  in  which  the 
LIBS  analytical  technique  is  evaluated  in  the  laboratory  under  in  situ  conditions.  It 
is  divided  into  five  chapters  that  cover  single  and  double  pulse  LIBS  and  delves  into 
the  parameters  that  must  be  optimized  for  the  detection  of  elements  in  high  pressure 
aqueous  solutions.  A  new  data  processing  scheme  for  dealing  with  the  inherent  vari¬ 
ability  of  laser-induced  plasmas  is  developed  in  this  thesis.  This  processing  scheme  is 
applied  to  all  data  presented  in  Chapters  4-6. 

Chapter  Two,  “Laser-induced  breakdown  spectroscopy  of  bulk  aqueous  solutions 
at  oceanic  pressures:  Evaluation  of  key  measurement  parameters,”  is  a  manuscript 
that  appeared  in  the  1  May  2007  issue  of  Applied  Optics  [40].  It  presents  preliminary 
investigations  on  the  feasibility  of  using  LIBS  to  detect  analytes  in  bulk  liquids  at 
oceanic  pressures.  This  work  was  completed  as  part  of  two  extensive  research  visits 
to  the  University  of  South  Carolina. 

Chapter  Three,  “Analysis  of  laser-induced  breakdown  spectroscopy  (LIBS)  spec¬ 
tra:  The  case  for  extreme  value  statistics,”  is  a  manuscript  that  has  been  accepted  for 
publication  by  Spectrochimica  Acta:  Part  B  [51].  It  presents  a  new  data  processing 
approach  for  LIBS  spectra. 

Chapters  Four  and  Five  are  complementary  chapters  that  look  at  the  detection  of 
analytes  in  bulk  aqueous  solutions  at  oceanic  pressures  using  single  pulse  (Chapter 
4)  and  double  pulse  (Chapter  5)  LIBS.  These  two  chapters  focus  on  the  optimization 
of  the  key  experimental  parameters  for  the  detection  of  analytes.  Chapter  Four  deals 
with  the  detection  of  three  elements:  sodium,  calcium,  and  manganese  and  the  in¬ 
terrelationship  of  pressure,  gate  delay,  and  pulse  energy.  Chapter  Five  concentrates 
on  the  detection  of  sodium,  calcium,  manganese,  potassium,  and  magnesium  and  the 
interrelationship  of  pressure,  gate  delay,  pulse  energies,  and  interpulse  delay.  In  both 
chapters,  calibration  curves  and  limits  of  detection  are  presented. 

Chapter  Six  presents  preliminary  investigations  into  matrix  effects  for  three  ele¬ 
ments:  sodium,  calcium,  and  potassium.  This  chapter  also  explores  the  effect  of  a 
chloride  versus  sulfate  background  matrix  on  the  detection  of  sodium  and  potassium. 


27 


Bibliography 

[1]  F.  Brech  and  L.  Cross.  Optical  microemission  stimulated  by  a  a  ruby  maser. 
Applied  Spectroscopy ,  16:59,  1962. 

[2]  R.  S.  Harmon,  F.  C.  DeLucia,  C.  E.  McManus,  N.  J.  McMillan,  T.  F.  Jenkins, 
M.  E.  Walsh,  and  A.  Miziolek.  Laser-induced  breakdown  spectroscopy  -  an 
emerging  chemical  sensor  technology  for  real-time  field-portable,  geochemical, 
mineralogical,  and  environmental  applications.  Applied  Geochemistry ,  21(5):730 
747,  May  2006. 

[3]  A.  S.  Eppler,  D.  A.  Cremers,  D.  D.  Hickmott,  M.  J.  Ferris,  and  A.  C.  Koskelo. 
Matrix  effects  in  the  detection  of  Pb  and  Ba  in  soils  using  laser-induced  break¬ 
down  spectroscopy.  Applied  Spectroscopy ,  50:1175  1181,  1996. 

[4]  D.  Anglos.  Laser-induced  breakdown  spectroscopy  in  art  and  archaeology.  Ap¬ 
plied  Spectroscopy ,  55:186A-205A,  2001. 

[5]  G.  Galbacs,  I.  B.  Gornushkin,  B.  W.  Smith,  and  .J.  D.  Winefordner.  Semi- 
quantitative  analysis  of  binary  alloys  using  laser-induced  breakdown  spectroscopy 
and  a  new  calibration  approach  based  on  linear  correlation.  Spectrochimica  Acta 
Part  B,  56:1159  1173,  2001. 

[6]  A.  C.  Samuels,  F.  C.  DeLucia,  Jr.,  K.  L.  McNesby,  and  A.  W.  Miziolek.  Laser- 
induced  breakdown  spectroscopy  of  bacterial  spores,  molds,  pollens,  and  protein: 
initial  studies  of  discrimination  potential.  Applied  Optics ,  42(30):6205  6209. 
2003. 

[7]  L.  St-Onge,  E.  Kwong,  M.  Sabsabi,  and  E.  B.  Vadas.  Quantitative  analysis  of 
pharmaceutical  products  by  laser-induced  breakdown  spectroscopy.  Spectrochim¬ 
ica  Acta  Part  B,  57(7) :  1 131-1 140,  2002. 

[8]  N.  Carmona,  M.  Oujjab,  E.  Rebollar,  H.  Romich,  and  M.  Castillejo.  Analysis  of 
corroded  glasses  by  laser  induced  breakdown  spectroscopy.  Spectrochimica  Acta 
Part  B,  2005. 

[9]  A.  I.  Whitehouse,  J.  Young,  I.  M.  Botheroyd,  S.  Lawson,  C.  P.  Evans,  and 
J.  Wright.  Remote  material  analysis  of  nuclear  power  station  steam  generator 
tubes  by  laser-induced  breakdown  spectroscopy.  Spectrochimica  Acta  Part  B, 
56:821  830,  2001. 

[10]  D.  Anglos,  S.  Couris,  and  C.  Fotakis.  Artworks:  laser-induced  breakdown  spec¬ 
troscopy  in  pigment  identification.  Applied  Spectroscopy ,  51:1025-1030,  1997. 

[11]  A.  Bertolini,  G.  Carelli,  F.  Francesconi,  M.  Francesconi,  L.  Marchesini,  P.  Mar- 
sili,  F.  Sorrentino,  G.  Cristoforetti,  S.  Legnaioli,  V.  Palleschi,  L.  Pardini,  and 
A.  Salvetti.  Modi:  a  new  mobile  instrument  for  in  situ  double-pulse  LIBS  anal¬ 
ysis.  Analytical  and  Bioanalytical  Chemistry,  V385(2):240-247,  May  2006. 


28 


[12]  M.  Corsi,  G.  Cristoforetti,  M.  Hidalgo,  S.  Legnaioli,  V.  Palleschi,  A.  Salvetti, 
E.  Tognoni,  and  C.  Vallebona.  Double  pulse,  calibration- free  laser-induced  break¬ 
down  spectroscopy:  A  new  technique  for  in  situ  standard-less  analysis  of  polluted 
soils.  Applied  Geochemistry ,  21  (5):748— 755,  2006. 

[13]  S.  Palanco,  C.  Lopez-Moreno,  and  J.  J.  Laserna.  Design,  construction  and  as¬ 
sessment  of  a  field-deployable  laser-induced  breakdown  spectrometer  for  remote 
elemental  sensing.  Spectrochimica  Acta  Part  B ,  61(1  ):88— 95,  January  2006. 

[14]  R.  S.  Harmon,  F.  C.  DeLucia,  A.  LaPointe,  R.  J.  Winkel,  and  A.  W.  Miziolek. 
LIBS  for  landmine  detection  and  discrimination.  Analytical  and  Bioanalytical 
Chemistry,  V385(6):  1140-1 148,  2006. 

[15]  Z.  A.  Arp,  D.  A.  Cremers,  R.  C.  Wiens,  D.  M.  Wayne,  B.  Salle,  and  S.  Mau¬ 
rice.  Analysis  of  water  ice  and  water  ice/soil  mixtures  using  laser-induced  break¬ 
down  spectroscopy:  Application  to  Mars  polar  exploration.  Applied  Spectroscopy, 
58:897-909,  2004. 

[16]  Z.  A.  Arp,  D.  A.  Cremers,  R.  D.  Harris,  D.  M.  Oschwald,  G.  R.  Parker  Jr., 
and  D.  M.  Wayne.  Feasibility  of  generating  a  useful  laser-induced  breakdown 
spectroscopy  plasma  on  rocks  at  high  pressure:  preliminary  study  for  a  Venus 
mission.  Spectrochimica  Acta  Part  B,  59:987-999,  2004. 

[17]  G.  B.  Courreges-Lacoste,  B.  Ahlers,  and  F.  R.  Perez.  Combined  Raman 
spectrometer/laser-induced  breakdown  spectrometer  for  the  next  ESA  mission 
to  Mars.  Spectrochimica  Acta  Part  A,  In  Press,  2007. 

[18]  R.  Brennetot.,  J.  L.  Lacour,  E.  Vors,  A.  Rivoallan,  D.  Vailhen,  and  S.  Maurice. 
Mars  analysis  by  laser-induced  breakdown  spectroscopy  (MALIS):  Influence  of 
Mars  atmosphere  on  plasma  emission  and  study  of  factors  influencing  plasma 
emission  with  the  use  of  Doehlert  designs.  Applied  Spectroscopy,  57(7):744  752, 
2003. 

[19]  A.  Knight,  N.  Scherbarth,  D.  Cremers,  and  M.  Ferris.  Characterization  of  laser- 
induced  breakdown  spectroscopy  (LIBS)  for  application  to  space  exploration. 
Applied  Spectrosopy,  54:331-340,  2000. 

[20]  B.  Salle,  J.-L.  Lacour,  P.  Mauchien,  P.  Fichet,  S.  Maurice,  and  G.  Manhes. 
Comparative  study  of  different  methodologies  for  quantitative  rock  analysis  by 
laser-induced  breakdown  spectroscopy  in  a  simulated  Martian  atmosphere.  Spec¬ 
trochimica  Acta  Part  B,  61:301-313,  2006. 

[21]  L.  Niu,  H.  Cho,  K.  Song,  H.  Cha,  Y.  Kim,  and  Y.  Lee.  Direct  determination 
of  strontium  in  marine  algae  samples  by  laser-induced  breakdown  spectrometry. 
Applied  Spectroscopy,  56:1511-1514,  2002. 

[22]  R.  Barbini,  F.  Colao,  V.  Lazic,  R.  Fantoni,  A.  Palucci,  and  M.  Angelone.  On 
board  LIBS  analysis  of  marine  sediments  collected  during  the  XVI  Italian  cam¬ 
paign  in  Antarctica.  Spectrochimica  Acta  B,  57:1203-1218,  2002. 


29 


[23]  A.  De  Giacomo,  M.  Dell’Aglio,  F.  Colao,  R.  Fantoni,  and  V.  Lazic.  Double-pulse 
LIBS  in  bulk  water  and  on  submerged  bronze  samples.  Applied  Surface  Science , 
247:157-162,  2005. 

[24]  A.  De  Giacomo,  M.  Dell’Aglio,  F.  Colao,  and  R.  Fantoni.  Double  pulse  laser 
produced  plasma  on  metallic  target  in  seawater:  basic  aspects  and  analytical 
approach.  Spectrochirnica  Acta  B ,  59:1431-1438,  2004. 

[25]  A.  De  Giacomo,  M.  Dell’Aglio,  A.  Casavola,  G.  Colonna,  0.  De  Pascale,  and 
M.  Capitelli.  Elemental  chemical  analysis  of  submerged  targets  by  double-pulse 
laser-induced  breakdown  spectroscopy.  Analytical  and  Bioanalytical  Chemistry , 
V385(2):303-311,  2006. 

[26]  A.  De  Giacomo,  M.  Dell’Aglio,  O.  De  Pascale,  and  M.  Capitelli.  From  single  pulse 
to  double  pulse  ns-laser  induced  breakdown  spectroscopy  under  water:  elemental 
analysis  of  aqueous  solutions  and  submerged  solid  samples.  Spectrochirnica  Acta 
Part  B ,  in  press,  2007. 

[27]  Y.  Lee,  K.  Song,  and  J.  Sneddon.  Laser-induced  breakdown  spectrometry.  Nova 
Science  Publishers,  Inc,  Huntington,  NY,  2000. 

[28]  X.  Hou  and  B.  T.  Jones.  Field  instrumentation  in  atomic  spectroscopy.  Micro- 
chemical  Journal  66:115  145,  2000. 

[29]  P.  K.  Kennedy,  S.  A.  Boppart,  D.  X.  Hammer,  B.  A.  Rockwell,  G.  D.  Noojin.  and 
W.  P.  Roach.  A  first-order  model  for  computation  of  laser-induced  breakdown 
thresholds  in  ocular  and  aqueous  media:  Part  II  -  comparison  to  experiment. 
IEEE  Journal  of  Quantum  Electronics ,  31:2250-2257,  1995. 

[30]  D.  A.  Cremers  and  L.  J.  Radziemski.  Laser-Induced  Breakdown  Spectroscopy 
(LIBS):  Fundamentals  and  Applications,  chapter  History  and  Fundamentals  of 
LIBS,  pages  1-39.  Cambridge  University  Press,  2006. 

[31]  D.  A.  Rusak,  B.  C.  Castle,  B.  W.  Smith,  and  J.  D.  Winefordner.  Recent  trends 
and  the  future  of  laser-induced  plasma  spectroscopy.  Trends  in  Analytical  Chem¬ 
istry,  17:453-461,  1998. 

[32]  J.  Ingle  Jr.  and  S.  Crouch.  Spectrochemical  Analysis.  Prentice  Hall,  New  Jersey, 
1988. 

[33]  R.  Nyga  and  W.  Neu.  Double-pulse  technique  for  optical  emission  spectroscopy 
of  ablation  plasmas  of  samples  in  liquids.  Optics  Letters,  18:747  749,  1993. 

[34]  R.  T.  Wainner,  R.  S.  Harmon,  A.  W.  Miziolek,  K.  L.  McNesby,  and  P.  D.  French. 
Analysis  of  environmental  lead  contamination:  Comparison  of  LIBS  field  and 
laboratory  instruments.  Spectrochirnica  Acta  Part  B,  56:777  793,  2001. 


30 


[35]  D.  A.  Rusak,  B.  C.  Castle,  B.  W.  Smith,  and  J.  D.  Winefordner.  Fundamentals 
and  applications  of  laser-induced  breakdown  spectroscopy.  Critical  Reviews  in 
Analytical  Chemistry ,  27:257-290,  1997. 

[36]  E.  Tognoni,  V.  Palleschi,  M.  Corsi,  and  G.  Cristoforetti.  Quantitative  micro¬ 
analysis  by  laser-induced  breakdown  spectroscopy:  A  review  of  the  experimental 
approaches.  Spectrochimica  Acta  Part  B,  57:1115-1130,  2002. 

[37]  D.  A.  Cremers,  L.  J.  Radziemski,  and  T.  R.  Loree.  Spectrochemical  analysis  of 
liquids  using  the  laser  spark.  Applied  Spectroscopy ,  38:721-729,  1984. 

[38]  R.  Knopp,  F.  J.  Scherbaum,  and  J.  I.  Kim.  Laser  induced  breakdown  spec¬ 
troscopy  (LIBS)  as  an  analytical  tool  for  the  detection  of  metal  ions  in  aqueous 
solutions.  Fresenius’  Journal  of  Analytical  Chemistry ,  355:16-20,  1996. 

[39]  W.  Pearman,  J.  Scaffidi,  and  S.  M.  Angel.  Dual-pulse  laser-induced  breakdown 
spectroscopy  in  bulk  aqueous  solution  with  an  orthogonal  beam  geometry.  Ap¬ 
plied  Optics ,  42:6085-6093,  2003. 

[40]  A.  P.  M.  Michel,  M.  Lawrence-Snyder,  S.  M.  Angel,  and  A.  D.  Chave.  Laser- 
induced  breakdown  spectroscopy  of  bulk  aqueous  solutions  at  oceanic  pressures: 
evaluation  of  key  measurement  parameters.  Applied  Optics ,  46,  2007. 

[41]  M.  Lawrence-Snyder,  J.  Scaffidi,  S.  M.  Angel,  A.  P.M.  Michel,  and  A.  D.  Chave. 
Sequential- pulse  laser-induced  breakdown  spectroscopy  of  high-pressure  bulk 
aqueous  solutions.  Applied  Spectroscopy,  61:171-176,  2007. 

[42]  M.  Lawrence-Snyder,  J.  Scaffidi,  S.  M.  Angel,  A.  P.  M.  Michel,  and  A.  D.  Chave. 
Laser-induced  breakdown  spectroscopy  of  high-pressure  bulk  aqueous  solutions. 
Applied  Spectroscopy,  60:786-790,  2006. 

[43]  A.  E.  Pichahchy,  D.  A.  Cremers,  and  M.  J.  Ferris.  Elemental  analysis  of  metals 
under  water  using  laser-induced  breakdown  spectroscopy.  Spectrochimica  Acta 
Part  B ,  52:25-39,  1997. 

[44]  C.  Haisch,  J.  Liermann,  U.  Panne,  and  R.  Niessner.  Characterization  of  colloidal 
particles  by  laser-induced  plasma  spectroscopy  (LIPS).  Analytica  Chimica  Acta, 
346:23-25,  1997. 

[45]  D.  A.  Cremers  and  L.  J.  Radziemski.  Handbook  of  Laser-Induced  Breakdown 
Spectroscopy.  John  Wiley  and  Sons,  Ltd.,  2006. 

[46]  P.  K.  Kennedy,  D.  X.  Hammer,  and  B.  A.  Rockwell.  Laser-induced  breakdown 
in  aqueous  media.  Progress  in  Quantum  Electronics,  21:155-248,  1997. 

[47]  C.  R.  German  and  K.  L.  Von  Damm.  Treatise  on  Geochemistry,  chapter  Hy¬ 
drothermal  Processes,  pages  181-222.  Elsevier,  2003. 


31 


[48]  K.  L.  Von  Damm.  Controls  on  the  chemisty  and  temporal  variability  of  seafloor 
hydrothermal  fluids.  Seafloor  hydrothermal  systems:  physical ,  chemical ,  biologi¬ 
cal,  and  geological  interactions:  Geophysical  Monograph  91,  pages  222  247,  1995. 

[49]  .1.  H.  Trefry,  D.  B.  Butterfield,  S.  Metz,  G.  J.  Massoth,  R.  R  Trocine,  and  R.  A. 
Feely.  Trace  metals  in  hydrothermal  solutions  from  Cleft  segment  on  the  southern 
Juan  de  Fuca  Ridge.  Journal  of  Geophysical  Research,  99:4925  4935,  1994. 

[50]  K.  L.  Von  Damm.  Chemistry  of  hydrothermal  vent  fluids  from  9°  -  10°  N, 
East  Pacific  Rise:  ‘Time  zero,’  The  immediate  posteruptive  period.  Journal  of 
Geophysical  Research,  105:11203-11222,  2000. 

[51]  A.  P.  M.  Michel  and  A.  D.  Chave.  Analysis  of  laser-induced  breakdown  spec¬ 
troscopy  (LIBS)  spectra:  The  case  for  extreme  value  statistics.  Spectrochimica 
Acta  Part  D,  In  Press. 


32 


Chapter  2 


Laser-induced  breakdown 
spectroscopy  of  bulk  aqueous 
solutions  at  oceanic  pressures: 
evaluation  of  key  measurement 
parameters 


The  work  in  this  chapter  is  published  in  the  1  May  2007  issue  of  Applied  Optics  as  A. 
P.  M.  Michel,  M.  Lawrence-Snyder,  S.  M.  Angel,  and  A.  D.  Chave,  “Laser-induced 
breakdown  spectroscopy  of  bulk  aqueous  solutions  at  oceanic  pressures:  evaluation  of 
key  measurement,  parameters.” 


2.1  Abstract 

The  development  of  in  situ  chemical  sensors  is  critical  for  present  day  expeditionary 
oceanography  and  the  new  mode  of  ocean  observing  systems  that  we  are  entering. 
New  sensors  take  a  significant  amount  of  time  to  develop;  therefore,  validation  of  tech¬ 
niques  in  the  laboratory  for  use  in  the  ocean  environment  is  necessary.  Laser-induced 
breakdown  spectroscopy  (LIBS)  is  a  promising  in  situ  technique  for  oceanography. 
Laboratory  investigations  on  the  feasibility  of  using  LIBS  to  detect  analytes  in  bulk 
liquids  at  oceanic  pressures  were  carried  out.  LIBS  was  successfully  used  to  detect 
dissolved  Na,  Mn,  Ca,  K,  and  Li  at  pressures  up  to  2.76  x  10'  Pa.  The  effects 
of  pressure,  laser  pulse  energy,  interpulse  delay,  gate  delay,  temperature,  and  NaCl 
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concentration  on  the  LIBS  signal  were  examined.  An  optimal  range  of  laser  pulse 
energies  was  found  to  exist  for  analyte  detection  in  bulk  aqueous  solutions  at  both 
low  and  high  pressures.  No  pressure  effect  was  seen  on  the  emission  intensity  for 
Ca  and  Na  and  an  increase  in  emission  intensity  with  increased  pressure  was  seen 
for  Mn.  Using  the  dual  pulse  technique  for  several  analytes,  a  very  short  interpulse 
delay  resulted  in  the  greatest  emission  intensity.  The  presence  of  NaCl  enhanced  the 
emission  intensity  for  Ca,  but  had  no  effect  on  peak  intensity  of  Mn  or  K.  Overall,  in¬ 
creased  pressure,  the  addition  of  NaCl  to  a  solution,  and  temperature  did  not  inhibit 
detection  of  analytes  in  solution  and  sometimes  even  enhanced  the  ability  to  detect 
the  analytes.  The  results  suggest  that  LIBS  is  a  viable  chemical  sensing  method  for 
in  situ  analyte  detection  in  high  pressure  environments  like  the  deep  ocean. 


2.2  Introduction 

Since  laser-induced  breakdown  spectroscopy  (LIBS)  was  first  reported  in  1962,  the 
technique  has  evolved  into  a  widely  used  method  for  laboratory  analytical  chemist  ry 
[1  -8].  Due  to  several  advantages  over  other  methods,  LIBS  has  been  identified  as 
a  viable  tool  for  in  situ  measurements,  especially  in  extreme  environments  [9,  10]. 
The  technique  yields  simultaneous  sensitivity  to  virtually  all  elements  in  the  parts- 
per-million  (ppm)  or  better  range  in  solids,  gases,  aerosols,  and  at  the  gas-liquid 
interface.  LIBS  is  effectively  non-invasive,  requiring  only  a  small  sample  (typically,  pg 
to  ng  of  material  are  ablated).  Unlike  for  many  analysis  techniques,  the  sample  does 
not  need  to  be  prepared.  LIBS  requires  only  optical  access  to  a  sample  and  therefore 
can  be  used  in  a  stand-off  mode  without  perturbing  the  sample  environment.  LIBS 
measurements  are  essentially  real-time,  with  typical  sampling  rates  of  less  than  one 
second  per  cycle.  These  characteristics  are  all  required  for  in  situ  chemical  sensing 
in  the  ocean  [1115]. 

Although  researchers  have  been  successful  at  inducing  plasma  ablation  on  sub¬ 
merged  materials  [16],  on  a  water  surface  or  film  [17  22],  and  in  liquid  jets,  droplets, 
and  flowing  solutions,  [23  29]  only  limited  LIBS  work  has  focused  on  analyte  de¬ 
tection  within  bulk  aqueous  solutions  [30-32].  Furthermore,  the  work  within  bulk 
aqueous  solutions  has  been  at  atmospheric  pressure.  Pioneering  work  by  Cremers 
et  al.  [30]  showed  that  LIBS  could  identify  Li,  Na,  K,  Rb,  Cs,  Be,  Ca,  B,  and  A1 
in  aqueous  solutions  with  varying  detection  limits,  but  typically  at  the  ppm  level. 
Several  studies  in  bulk  liquids  have  displayed  a  reduction  in  the  time  during  which 
plasma  emission  can  be  observed  as  compared  to  that  in  air  [16,  30,  31,  33].  The 
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plasma  lifetime  is  typically  <  1  /xs  in  bulk  liquids,  whereas  at  an  air-liquid  interface 
it  averages  5-20  fis.  Laser-induced  plasmas  formed  in  solution  are  also  characterized 
by  a  reduction  in  plasma  light  intensity. 

The  effects  of  elevated  pressure  and  temperature  on  LIBS  spectra  have  received 
limited  attention.  Although  a  few  researchers  report  on  LIBS  at  super-atmospheric 
pressures,  they  do  not  extend  beyond  1  x  10'  Pa  (note:  1  Pa  =  1  x  10~  ’  bars),  which 
is  well  below  the  ambient  pressure  in  the  deep  ocean  [9,  34];  yet,  none  of  these  studies 
were  for  liquids.  Although,  we  have  previously  reported  preliminary  findings  that 
show  the  ability  to  detect  analytes  in  high  pressure  bulk  aqueous  solutions  [35],  we 
now  focus  on  the  key  measurement  parameters  that  are  needed  for  analyte  detection. 
The  influence  of  in  situ  temperature  is  anticipated  to  be  weak  because  of  the  high 
plasma  temperature  («  8000K  at  early  times)  [36-39]. 

For  many  years,  oceanography  has  been  in  an  expeditionary  mode  where  research 
vessels  are  used  for  short  term  instrument  deployments  with  limited  resolution  in 
time.  Although  oceanographers  will  continue  to  study  the  ocean  in  this  way,  a  new 
paradigm  using  ocean  observatories  for  long  term  in  situ  observing  is  upon  us.  As 
this  shift  towards  long  term  ocean  observing  systems  becomes  recognized,  we  must 
also  acknowledge  the  need  for  in  situ  sensors;  especially  those  capable  of  temporal 
studies.  A  major  need  is  for  chemical  sensors.  The  development,  of  new  sensors  for 
oceanography  takes  a  significant  amount  of  time,  and  hence  laboratory  validation  of 
techniques  such  as  LIBS  is  necessary  now  to  identify  techniques  that  are  viable  for 
chemical  detection  in  high  pressure,  high  salinity,  aqueous  environments. 

Although  LIBS  has  the  potential  for  use  in  numerous  ocean  environments,  and 
has  applicability  to  solids  and  liquids,  we  have  focused  on  the  feasibility  of  detect¬ 
ing  elements  at  one  extreme  ocean  environment,  hydrothermal  vents.  Hydrothermal 
venting  occurs  on  mid-ocean  ridges  where  seawater  circulates  through  the  fractured 
and  permeable  oceanic  crust.  Exit  temperatures  at  discrete  (orifice  diameters  of  a 
few  centimeters)  high-temperature  vents  range  from  200-405°C  at  ambient  pressures 
of  1.5  x  107  to  3.7  x  107  Pa.  Low-temperature  (usually  <  35°C)  diffuse  flow  seeping 
from  porous  surfaces  or  cracks  is  frequently  observed  [40].  The  circulation  is  driven 
by  the  direct  or  indirect  thermal  effects  of  magma  at  sub-seafloor  depths  of  up  to 
a  few  km.  Substantial  changes  in  fluid  composition  occur  due  to  interaction  with 
the  host  rock,  phase  separation  into  a  mixed  liquid-vapor  form,  and  possibly  magma 
degassing.  Many  alkalis  (e.g.,  Li,  Na,  and  Ca)  and  transition  metals  (e.g.,  Fe,  Mn, 
Cu,  and  Zn)  are  leached  from  the  host  rock  and  concentrated  to  varying  degrees  in 
the  hydrothermal  fluid,  while  Mg  and  S04  are  largely  removed  from  the  fluid  by  pre- 
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cipitation  into  Mg-OH-Si  minerals  and  anhydrite,  respectively  [40].  Von  Dannn  and 
Butterfield  et  al.  provide  comprehensive  reviews  of  the  chemistry  of  hydrothermal 
vent  fluids  [40,  41]. 

In  this  paper,  we  explore  the  effect  of  vent  system  environmental  factors  such 
as  pressure,  temperature,  and  NaCl  concentration  on  the  LIBS  signal  to  assess  the 
feasibility  of  developing  LIBS  for  in  situ  chemical  sensing  in  the  ocean.  In  addition, 
several  system  parameters  (laser  energy  per  pulse,  interpulse  spacing,  and  gate  delay) 
are  optimized  for  high  pressures  for  the  first  time. 

2.3  Experiment 

A  laboratory  LIBS  system  was  designed  to  operate  with  a  high  pressure  cell  (Figure  2- 
1).  For  single  pulse  experiments,  a  Continuum  Surelite  III  laser  (5-ns  pulse  width)  was 
utilized.  For  dual  pulse  experiments,  a  Quantel  Nd  580  (9-ns  pulse  width)  was  used  for 
the  first  laser  pulse  followed  by  a  second  pulse  from  the  Surelite  laser.  Both  lasers  were 
Q-switched  Nd:YAG  types  operated  at  the  fundamental  wavelength  with  a  repetition 
rate  of  5  Hz.  For  dual  pulse  experiments,  a  variable  clock  (Stanford  Instruments 
Model  SR250)  and  a  delay  generator  (Stanford  Instruments  Model  DG535)  controlled 
laser  triggering. 

The  laser  pulses  were  focused  into  a  high  pressure  cell,  designed  to  reach  pressures 
of  3.45  x  10'  Pa  and  constructed  of  stainless  steel  Swagelok  fittings  with  six  1"-ID 
(1  in.  =  2.54  cm)  and  l-^-OD  ports.  Stainless  steel  tubing  (k")  connected  one  port 
to  a  pump  (Isco  Syringe  Pump  Model  260D,  Teledyne  Technologies  Incorporated) 
that  allowed  aqueous  solutions  to  flow  into  the  cell  and  the  cell  to  be  pressurized. 
A  second  port  was  equipped  with  the  same  tubing  and  a  regulating  valve  for  cell 
drainage.  Two  ports  were  fitted  with  1"  diameter,  j"  thick  circular  sapphire  windows 
(MSW100/125,  Meller  Optics  Incorporated)  held  in  place  by  hex  nuts  and  sealed 
with  rubber  washers,  allowing  of  each  window  to  be  visible  outside  the  cell.  The 
remaining  two  ports  were  sealed  with  Swagelok  plugs  (SS-1610-P). 

Three  different,  optical  arrangements  for  focusing  the  laser  pulses  into  the  cell  and 
for  collection  of  the  plasma  emission  were  used  in  these  experiments,  as  detailed  in 
Figure  2-2.  Because  the  purpose  of  these  experiments  was  initial  investigation  into  the 
feasibility  of  using  LIBS  for  ocean  applications,  one  of  the  goals  was  determining  the 
best  optical  set-up.  For  dual  pulse  operation,  the  lasers  were  collinear.  In  some  single 
pulse  configurations,  light  collection  was  collinear  to  the  laser  pulse,  while  in  others  it 
was  orthogonal  for  ease  of  alignment.  All  optics  were  mounted  on  micrometer  stages, 
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Figure  2-1:  Schematic  of  the  laboratory  LIBS  apparatus.  Note  that  in  the  drawing, 
the  laser  pulses  are  simply  represented  by  arrows  as  their  optical  paths  are  described 
in  Figure  2-2. 
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Figure  2-2:  Optical  arrangements  used  in  experiments  showing  the  high  pressure  cell 
with  respect  to  incoming  laser  pulses  (signified  by  a  dashed  line).  FO  =  Optical 
Fiber,  (a)  L[,  L2,  and  L3  =  // 4  lenses;  Mi  =  dielectric  coated  mirror,  (b)  L|=//4,  To 
study  the  effect  of  NaCl  concentration  on  spectra:  L2  =  // 3  lens,  L.{  =  // 2  lens.  To 
study  the  detection  of  Ca  at  varying  concentrations:  L_>  =  // 4  lens,  L3  =  // 3  lens,  (c) 
Li=//4,  Mi  and  M2  =  parabolic  off-axis  mirrors 


enabling  precise  control  of  beam  overlap  and  collection  field  of  view  within  the  high 
pressure  cell.  All  lenses  were  made  of  fused  silica.  In  all  optical  configurations, 
the  plasma  emission  was  focused  onto  a  2-mm-core-diameter,  0.51-N.A.  light  guide 
(Edmund  Scientific  Co.  Model  02551).  The  light  guide  was  connected  to  a  0.25-m, 
// 4  spectrograph  (Chromex  model  250is/RF)  with  a  1200-groove/ mm  grating  blazed 
at  500  nm.  The  slit  width  (W)  ranged  from  25  -  250  /mi.  Data  were  collected 
on  an  intensified  CCD  detector  (Princeton  Instruments,  I-Max  1024E)  and  acquired 
with  a  computer  running  WinSpec/32  software.  All  spectra  were  accumulations  of 
250  shots  at  the  maximum  gain  setting  of  255.  All  error  bars  represent  ±\cr.  A 
similar  apparatus  and  set-up  was  previously  used  to  demonstrate  the  feasibility  of 
high  pressure  LIBS  [35]. 

The  key  LIBS  timing  parameters  have  been  previously  described  [10,  32],  The 
first  and  second  laser  pulse  energies  are  referred  to  as  Ei  and  E2.  For  dual  pulse 
experiments,  the  time  interval  between  the  two  pulses  or  interpulse  delay,  is  referred 
to  as  AT.  The  gate  delay,  t^,  is  the  time  between  the  last  laser  pulse  and  turn-on  of 
the  detector.  The  plasma  emission  is  recorded  by  the  detector  for  the  length  of  time 
set  by  the  gate  width,  t.6,  which  was  set  at  1  /rs  for  all  experiments  reported  here. 
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Laser  beam  waist  width  d„o can  be  estimated  by: 


d<r0 


4/AM2 
7 tD 


(2.1) 


where  f  is  the  focal  length  of  the  focusing  lens  (100  mm),  A  is  the  laser  wavelength 
(1064  nm),  M2  is  the  beam  propagation  ratio  which  is  typically  2-10  for  Nd:YAG 
lasers  (we  use  a  value  of  6),  and  D  is  the  diameter  of  the  illuminated  aperture  of  the 
focusing  lens  («  25  mm)  [42].  The  beam  waist  width  for  the  system  is  approximately 
0.03  mm.  The  average  irradiance  (I/)  at  the  beam  waist  can  be  estimated  using: 

'  4r,,/'-'A2A/'  (  '  1 

where  is  the  laser  pulse  energy  and  tl,  is  the  pulse  duration  at  the  full  peak  width 
at  half  of  the  maximum  intensity  (FWHM)[42]  (for  the  Continuum  laser,  =  5  ns, 
and  for  the  Quantel  laser,  tl  =  9  ns).  The  pulse  energies  of  the  Continuum  laser 
used  vary  between  «  10  -  100  mJ.  The  irradiance  of  the  beam  at  the  beam  waist  thus 
varies  from  ss  2.4  x  1011  to  ~  2.4  x  1012  W/cm2.  The  pulse  energies  of  the  Quantel 
laser  used  vary  between  «  10  -  125  mJ  with  the  irradiance  thus  varying  from  1.3  x 
1011  to  1.7  x  1012  W/cm2. 

Sample  solutions  were  made  by  dissolving  NaCl,  CaCL,  LiCl,  and  MnSOi-H^O  in 
deionized  water.  Where  noted,  NaCl  was  added  to  the  solutions  to  simulate  a  seawater 
environment.  All  concentrations  are  listed  in  parts  per  million  (ppm  wt./vol.). 


2.4  Results  and  Discussion 

2.4.1  The  Effect  of  Pulse  Energy  on  LIBS  Emission 

Single-Pulse  LIBS 

Two  key  constraints  on  the  design  of  an  oceanographic  sensor  system  are  instrument 
power  consumption  and  form  factor,  both  of  which  must  be  minimized.  LIBS  opera¬ 
tion  with  a  small,  low  power  laser  would  simplify  the  design  of  an  oceanographic  LIBS 
instrument.  The  effect  of  pulse  energy  on  signal  intensity  for  analytes  in  solution  at 
elevated  pressure  was  investigated  with  the  goal  of  minimizing  power  consumption. 
The  peak  signal  intensity  for  four  analytes  (Li,  Ca,  Na,  and  Mn)  was  measured  at 
laser  pulse  energies  ranging  from  11  to  91  mJ  at  both  low  (7  x  105  Pa)  and  high  (2.76 
x  107  Pa)  pressures  using  the  collinear  optical  configuration  shown  in  Figure  2-2(a) 
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Figure  2-3:  Effect  of  laser  pulse  energy  on  the  LIBS  signal  intensity  of  100  ppm  Na(I) 
(588.995  nm). 

(td  =  350  ns,  W  =  75  /mi  for  Na,  Mn,  and  Ca  studies,  and  W  =  250  nm  for  Li).  Ten 
spectra  were  recorded  and  averaged  for  each  condition. 

Figure  2-3  shows  the  dependence  of  the  Na(I)  (588.995  nm)  emission  line  on  laser 
pulse  energy  for  100  ppm  Na.  In  both  low  and  high  pressure  experiments,  as  pulse 
energy  increases,  a  corresponding  increase  in  peak  intensity  occurs  until  a  maximum 
intensity  is  reached  at  22  mJ  (Figure  2-3(a)).  Above  this  value,  emission  intensity 
decreases  sharply  up  to  ~  50  mJ,  above  which  a  more  gradual  decrease  with  energy 
is  observed.  These  data  suggest  that,  independent  of  pressure,  a  low  laser  pulse 
energy  yields  greater  emission  intensity  providing  the  energy  exceeds  a  threshold 
value.  Figure  2-3(b)  compares  spectra  taken  at  laser  pulse  energies  below,  above, 
and  in  the  optimal  energy  range  for  Na.  The  top  trace  (22  mJ)  shows  a  significantly 
greater  intensity  than  at  either  a  very  low  (middle  trace,  11  mJ)  or  a  high  (bottom 
trace,  88  mJ)  pulse  energy. 
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The  effect  of  laser  pulse  energy  on  Ca  (422.673  nm)  and  Li  (670.776  nm  and 
670.791  nm,  unresolved  doublet)  emission  displayed  similar  trends.  When  less  than 
14  mJ  was  used,  Ca  was  virtually  undetectable.  As  the  pulse  energy  was  increased 
above  this  level,  emission  intensified  until  a  maximum  was  achieved  at  36  mJ  for  low 
pressure  (7  x  105  Pa)  and  at  29  mJ  for  high  pressure  (2.76  x  107  Pa).  This  range 
for  both  the  low  and  high  pressure  environments  was  ~  25  to  50  mJ.  At  energy  levels 
beyond  the  optimal  range,  intensity  decreased  slowly  with  increasing  pulse  energy, 
possibly  due  to  plasma  shielding.  Plasma  shielding  occurs  when  the  plasma  itself 
reduces  the  transmission  of  the  laser  pulse  energy  along  the  beam  path.  Calcium 
displayed  a  more  gradual  increase  and  then  decrease  in  intensity  and  a  wider  range 
of  optimal  energy  compared  to  Na.  Similar  trends  were  observed  for  Li.  At  both  low 
and  high  pressures,  plasma  emission  was  not  detectable  below  11  mJ.  At  higher  pulse 
energies  and  both  pressures,  the  emission  maximum  was  recorded  at  27  mJ,  above 
which  a  sharp  decrease  in  intensity  to  46  mJ  was  observed,  followed  by  flattening  to 
72  mJ. 

The  relationship  between  emission  intensity  and  laser  pulse  energy  for  the  unre¬ 
solved  403  nm  Mn(I)  triplet  was  slightly  different  than  for  the  other  three  analytes. 
Figure  2-4(a)  shows  that  the  lowest  laser  pulse  energy  (11  mJ)  resulted  in  the  high¬ 
est  emission  intensity.  At  pulse  energies  greater  than  11  mJ,  the  emission  intensity 
gradually  decreased  until  it  was  no  longer  detectable  above  ~  40  mJ  and  ~  70  mJ 
for  low  and  high  pressures,  respectively.  The  peak  intensity  was  greater  at  high  than 
at  low  pressure.  Figure  2-4(b)  compares  spectra  taken  at  11,  22,  and  88  mJ  at  2.76 
x  107  Pa. 

The  data  for  Na,  Ca,  Li,  and  Mn  suggest  that  the  pulse  energy  required  to  op¬ 
timize  the  LIBS  signal  is  analyte-dependent  due  to  different  ionization  energies,  but 
is  minimally  pressure  dependent.  A  pulse  energy  threshold  is  also  observed.  For  the 
four  analytes  studied,  a  relatively  low  laser  pulse  energy  (less  than  50  mJ)  produced 
the  greatest  signal  intensity.  A  low  energy  optimal  range  may  exist  due  to  effects 
from  plasma  shielding  or  moving  breakdown.  Plasmas  can  expand  back  along  the 
laser  beam  path  towards  the  laser  resulting  in  elongated  plasmas  [43].  A  higher  en¬ 
ergy  pulse  may  form  a  more  elongate  plasma  or  a  series  of  plasmas  as  the  breakdown 
threshold  of  the  liquid  is  exceeded  before  the  pulse  reaches  the  focal  point.  This  may 
result  in  non-optimal  collection  of  the  plasma  emission.  Further  studies  using  imaging 
techniques  are  needed  to  elucidate  the  effect  of  pulse  energy  on  the  plasma. 

Figure  2-5  shows  the  effect  of  pressure  on  the  LIBS  signal  for  Na  (588.995  nm),  Mn 
(403  nm  unresolvable  triplet)  and  Ca  (422.673  nm)  using  a  low  energy  single  pulse. 
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(b)  Mn(I)  spectra  taken  at  2.76  x  107  Pa.  Spectra  offset  for  clarity 


Figure  2-4:  Effect  of  laser  pulse  energy  on  the  LIBS  emission  intensity  of  the  unre- 
solvable  Mn  (I)  triplet  (403  run)  (5,000  ppm  Mil  in  2,540  ppm  NaCl). 


42 


Figure  2-5:  Effect  of  pressure  on  LIBS  emission  intensity.  □  =  100  ppm  Na  (588.995 
urn)  with  E  =  22  mJ;  •  =  5,000  ppm  Mn  (403  nm  unresolvable  triplet)  with  2,540 
ppm  NaCl,  E  =  14  mJ;  A  =  500  ppm  Ca  (422.673  nm)  with  2,540  ppm  NaCl,  E  = 
20  mJ. 

The  gate  delay  was  fixed  at  350  ns  and  the  slit  width  was  fixed  at  75/tm.  Na  and 
Ca  display  no  change  in  signal  intensity  with  increasing  pressure,  but  Mn  shows  an 
increase.  For  all  analytes  examined,  the  FWHM  did  not  change  with  pressure.  Pres¬ 
sure  under  oceanic  conditions  does  not  induce  a  deleterious  effect  on  signal  intensity 
or  on  FWHM. 

In  these  single  pulse  energy  experiments,  the  same  gate  delay  and  gate  width  were 
used  for  all  energy  levels  and  pressures.  As  discussed  later  in  this  paper,  optimal  gate 
delay  may  be  energy  dependent.  Optimal  gate  width  was  not  investigated  and  may 
be  pressure  and/or  energy  level  dependent.  As  a  result,  the  selected  gate  width  and 
gate  delay  may  influence  the  measured  emission  intensity.  Optimal  gate  delay  could 
also  be  analyte  dependent,  and  hence  a  different  gate  delay  could  yield  another  trend 
with  pulse  energy.  However,  the  selected  conditions  demonstrate  that  low  energy 
single  laser  pulses  at  high  pressures  are  viable  for  measuring  analytes  in  bulk  aqueous 
liquids.  This  is  promising  toward  the  development  of  an  ocean-going  instrument 
where  a  small,  low  power  laser  will  be  critical. 

Dual  Pulse  LIBS  at  High  Pressure 

An  evaluation  of  dual  pulse  LIBS  for  high  pressure  bulk  solutions  shows  that  analyte 
detection  is  highly  dependent  on  the  interpulse  delay.  If  the  interpulse  delay  is  short 
(<C  1  ps),  signal  intensity  is  greatly  enhanced  when  compared  to  that  measured  using 
longer  delay  times.  However,  such  a  small  interpulse  delay  may  not  be  sufficient 
for  a  cavitation  bubble  to  fully  form  before  the  second  laser  pulse  creates  a  spark. 
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Dual  pulse  LIBS  has  been  shown  to  enhance  the  signal  intensity  for  some  analytes  in 
bulk  aqueous  solutions  at  atmospheric  pressure  [30,  32].  However,  such  enhancements 
using  longer  interpulse  delay  times  do  not  occur  for  high  pressure  liquids. 

To  demonstrate  the  coupled  effect  of  interpulse  delay  and  pulse  energy  on  emission 
intensity,  four  energy  level  conditions  were  compared  for  four  analytes  at  high  pressure 
(2.76  x  10'  Pa)  over  a  range  of  interpulse  delay  times.  The  four  conditions  were:  1) 
low  Ei,  low  E2  (low-low),  2)  high  Elt  high  E2  (high-high),  3)  low  Ei,  high  E2  (low- 
high),  and  4)  high  Ej ,  low  E2  (high-low),  and  are  detailed  in  Table  2.1,  (t d  =  350  ns, 
4  —  1  pis).  These  experiments  were  completed  using  the  optical  configuration  shown 
in  Figure  2-2(a). 

For  Ca  (W  =  l()0//m),  using  a  low  E|  followed  by  a  low  E2  resulted  in  the  highest 
peak  intensity,  possibly  because  when  summed  they  give  a  low  total  energy  (Figure 
2-6).  The  greatest  emission  is  observed  for  Ei  =  13  mJ  and  E2  =  6  mJ  and  yields  the 
ionic  Ca  peaks  (393.366  mil  and  396.847  nm  in  addition  to  the  atomic  peak  422.673 
nm).  For  this  low-low  condition,  Figure  2-7(a)  shows  the  emission  intensity  change 
with  AT.  For  AT  greater  than  1  pis,  the  intensity  remained  stable  at  a  value  of  1.5 
x  1()4  a.u. 

For  AT  less  than  1  /ts,  the  low-low  configuration  yielded  intensities  between  2.5 
x  104  and  8.7  x  1()4  a.u.  Figure  2-7(b)  compares  spectra  at  very  short  (30  ns  -  upper 
trace)  to  long  (30  /./.s  -  lower  trace)  AT  values.  When  a  short  AT  is  used,  three  Ca 
peaks  (Ca(II)  393.366  nm,  Ca(II)  396.847  nm,  and  Ca(I)  422.673  nm)  are  visible, 
while  for  long  AT,  only  the  Ca(I)  peak  is  present  with  a  much  lower  intensity.  When 
AT  is  30  ns,  the  low-low  configuration  yields  significantly  greater  emission  intensity 


Table  2.1:  Conditions  used  to  study  the  effect  of  dual  pulse  energies  on  LIBS  emission 


Low  Ei 
Low  E2 

High  Ei 
High  E2 

Low  Ei 
High  E2 

High  E, 
Low  E2 

Ei 

e2 

E, 

e2 

Ei 

E2 

E, 

e2 

(mJ) 

(mJ) 

(mJ) 

(mJ) 

(mJ) 

(mJ) 

(mJ) 

(mJ) 

1,000  ppm  Ca, 
2,540  ppm  NaCl 

13 

6 

105 

84 

13 

84 

105 

6 

100  ppm  Li 

31 

20 

105 

84 

31 

84 

105 

20 

100  ppm  Na 

13 

6 

105 

84 

13 

84 

105 

6 

5,000  ppm  Mil, 
2,540  ppm  NaCl 

13 

6 

105 

84 

13 

84 

105 

6 
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Figure  2-6:  Spectra  of  1,000  ppm  Ca  with  2,540  ppm  NaCl  at  2.76  x  10'  Pa  under 
four  dual  pulse  conditions. 

than  the  other  pulse  energy  configurations.  For  Li  (W  =  250  /tm),  a  low  Ej  followed 
by  a  low  E2  resulted  in  the  greatest  emission  intensity.  Table  2.2  shows  peak  emission 
for  Li  for  four  different  dual  pulse  conditions  for  AT  between  50  ns  and  1  /ts. 

A  small  delay  time  (<  1  /ts)  enhanced  the  emission  as  compared  to  a  longer  delay 
time  when  the  low-high  and  low-low  energy  levels  were  used.  For  Na(I),  the  low-high 
and  low-low  conditions  yielded  similar  intensities  at  all  delay  times,  with  maximum 
values  of  9.3  x  105  a.u.  and  8.3  x  105  a.u.,  respectively,  (W  =  75  /tm).  After  these 
four  runs  were  compared,  an  additional  configuration  consisting  of  a  13  mJ  first  pulse 
followed  by  a  22  mJ  second  pulse  was  tested  as  a  low-low  dual  pulse  condition  with  a 
slightly  increased  second  pulse  energy.  This  resulted  in  peaks  with  intensities  of  2.4 
x  106  -  2.9  x  106  a.u.  for  all  AT  values  between  10  ns  and  100  /ts,  suggesting  again 
that  a  low- low  energy  condition  produces  the  greatest  emission  intensity.  For  Mn  (W 
=  250  /tm),  at  all  interpulse  delay  times  between  20  11s  and  100  /ts,  a  low  E|  followed 
by  a  high  E2,  resulted  in  the  highest  emission  intensity  (Figure  2-8). 

These  results  show  that  the  best  dual  pulse  conditions  vary  by  analyte.  However, 


Table  2.2:  Dual  pulse  emission  intensity  (a.u.) 


Condition  1 

Condition  2 

Condition  3 

Condition  4 

Low  El 
Low  E2 

High  El 
High  E2 

Low  El 
High  E2 

High  El 
Low  E2 

100  ppm  Li 

2.5  x  105  - 
3.7  x  105 

2  x  104  - 

7  x  104 

5  x  104  - 

1.5  x  105 

1.5  x  104  - 

5.3  x  104 

5,000  ppm  Mn, 
2,540  ppm  NaCl 

4  x  103  - 
4.2  x  105 

2.1  x  105  - 

5  x  105 

7  x  105 - 

8.3  x  105 

1.6  x  103  - 

2.6  x  105 
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(a)  Effect  of  dual  laser  pulse  energies  on  emission  intensity  at  2.76  x  10'  Pa  for  1,000 
ppm  Ca  in  2,540  ppm  NaCl  at  various  interpulse  delays.  Each  data  point  is  the 
average  of  5  spectra. 


(b)  Spectra  of  Ca  showing  the  enhancement  in  signal  for  AT  =  60  ns  over  AT  =  30 
/x  s. 


Figure  2-7:  Ca  emission  at  different  interpulse  delays.  (Ei  =  13  mJ,  E2  =  6  mJ.) 
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Figure  2-8:  Spectra  of  5,000  ppm  Mn  with  2,540  ppm  NaCl  at  2.76  x  10'  Pa  under 
four  dual  pulse  conditions.  The  highest  emission  intensity  is  observed  for  a  low-high 
pulse  combination. 

it  is  important  to  note  that  at  high  pressure,  very  short  interpulse  spacing  results  in 
a  higher  signal  intensity  than  when  dual  pulses  are  separated  by  a  more  significant 
delay  in  time.  The  need  for  such  rapid  firing  of  the  two  pulses  is  only  accomplished 
using  two  independent  lasers  instead  of  firing  one  laser  rapidly.  Two  pulses  separated 
by  a  short  AT  approaches  single  pulse  conditions,  suggesting  that  dual  pulse  LIBS 
may  not  be  advantageous  at  elevated  pressure. 

2.4.2  Interrelationship  of  pulse  energy,  gate  delay,  and  pres¬ 
sure  for  Lithium 

Emission  intensity  was  recorded  for  the  unresolved  Li(I)  doublet  (670.776  mn  and 
670.791  nm)  at  two  single  pulse  energies  (27  mJ  and  68  mJ)  at  both  low  (7  x  105 
Pa)  and  high  pressure  (2.76  x  107  Pa)  over  a  range  of  gate  delays  (0.1  -  3.7  ^s),  (4 
=  1  fj, s,  W  =  25  nm)  using  the  optical  configuration  of  Figure  2-2(a).  Comparing  the 
two  curves  in  Figure  2-9,  it  is  clear  that  a  short  gate  delay  should  be  used  to  enhance 
emission  intensity.  These  results  also  suggest  that  the  optimal  gate  delay  may  be 
pulse  energy  but  not  pressure  dependent. 

2.4.3  Effect  of  NaCl  Concentration  on  LIBS  Spectra 

Understanding  how  pervasive  Na  and  Cl  ions  in  solution  affect  the  detection  of  other 
analytes  is  important  for  assessing  the  feasibility  of  using  LIBS  in  the  ocean,  where 
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Figure  2-9:  Effect  of  gate  delay  on  the  LIBS  signal  for  1,000  ppm  Li  (070  nm  unre- 
solvable  doublet).  7  x  10  5  Pa,  A=  2.57  x  107  Pa. 
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Figure  2-10:  Effect  of  the  addition  of  NaCl  in  solution  on  spectra  of  1,000  ppm  Ca  at 
2.57  x  107  Pa. 

the  nominal  concentrations  of  Na  and  Cl  are  1.08  x  104  ppm  and  1.95  x  104  ppm, 
respectively  [44],  Cremers  et.  al.  previously  reported  a  decrease  in  the  intensity  ratio 
of  Ca  Il/Ca  I  with  the  addition  of  NaCl  [30].  The  peak  signal  intensity  for  three 
analytes  (1,000  ppm  Ca,  100  ppm  Mn,  and  1,000  ppm  K)  was  compared  in  three 
solutes:  1)  deionized  water,  2)  2,540  ppm  NaCl  dissolved  in  deionized  water,  and 
3)  25,400  ppm  NaCl  dissolved  in  deionized  water  using  the  optical  configuration  of 
Figure  2-2(b)  and  for  a  range  of  pressures  (3  x  105  ,  7  x  105,  1.7  x  106  ,  3.4  x  10(',  6.9 
x  106,  1.38  x  107,  2.07  x  107,  2.76  x  107  Pa).  These  studies  were  carried  out  with  Ei 
=  40  mJ,  E2  =  125  mJ,  AT  —  46  ns,  W  =35  //ra,  and  t<*  =  100  ns  for  Ca  and  K  and  trf 
=  200  ns  for  Mn.  The  addition  of  NaCl  significantly  increased  the  emission  intensity 
of  the  422.673  nm  Ca(I)  atomic  line  whereas  no  significant  effect  was  seen  on  the 
393.366  nm  and  396.847  nm  Ca(II)  ionic  lines  (Figure  2-10).  The  signaknoise  ratio 
for  the  Ca  (II)  ionic  lines  showed  no  significant  change  with  the  addition  of  NaCl, 
whereas  the  signaknoise  ratio  of  Ca(I)  increased  from  22  to  30  with  the  addition  of 
NaCl.  The  same  increase  was  seen  with  the  addition  of  254  ppm  NaCl  and  25,400 
ppm  NaCl.  In  atomic  emission,  adding  an  easily  ionizable  element,  for  example  Na, 
can  suppress  ionization  of  other  elements.  This  suggests  that  ionization  suppression 
by  Na  increases  the  Ca(I)  emission  relative  to  the  Ca(II)  lines.  No  intensity  change 
was  seen  for  Mn(I)  (403  nm  unresolvable  triplet)  or  K(I)  (766.491  nm  and  769.897 
nm).  However,  since  only  atomic  lines  were  detectable  for  Mn  and  K,  the  relative 
increase  of  atomic  to  ionic  lines  could  not  be  compared. 

These  two  outcomes  (enhancement  of  the  signal  or  no  change  to  the  signal)  suggest 
that  the  high  NaCl  concentration  in  the  ocean  will  not  have  a  deleterious  effect  on 
the  ability  to  detect  Ca,  Mn,  and  K  analytes.  It  also  suggests  that  further  work  is 
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Figure  2-11:  Detection  of  Ca  (422.673  nm)  in  a  simulated  vent  fluid  at  varying  pres¬ 
sures  and  concentrations.  7  x  10  5  Pa,  A  =  7  x  106  Pa,  □  =  2.76  x  10' 

Pa. 


needed  to  elucidate  the  effect  NaCl  has  on  other  analytes. 


2.4.4  Detection  of  Calcium  at  Varying  Concentrations 

Ca  was  used  to  determine  whether  increased  pressure  affects  the  limit  of  detection. 
Five  pressures  ranging  from  7  x  105  Pa  to  2.76  x  10'  Pa  were  investigated  at  con¬ 
centrations  ranging  from  10  to  5,000  ppm  in  a  solution  containing  several  analytes 
(69  ppm  Br,  10,828  ppm  Na,  89  ppm  Fe,  958  ppm  K,  46  ppm  Mu,  18,932  ppm  Cl 
in  DI  water)  found  in  hydrothermal  vent  fluids  at  representative  concentrations  (10 
ppm,  25  ppm,  50  ppm,  100  ppm,  500  ppm,  1000  ppm,  5000  ppm),  using  the  optical 
configuration  of  Figure  2-2(b)  (E^  -  31  mJ;  E2  =  15  mJ;  AT  =  72  11s;  t.^  =  700  11s;  t 
=  1  /ts,  =  35  /an).  Figure  2-11  shows  that  varying  concentrations  of  Ca  are  detectable 
at  pressure  and  with  a  minimal  change  in  intensity,  suggesting  that  detection  limits 
to  the  ppm  level  will  be  obtainable  at  high  pressure. 


2.4.5  Solution  Temperature  Effects  on  Calcium  Spectra 

To  characterize  the  temperature  effect  for  Ca  spectra,  the  sample  cell  was  placed  in  a 
sand  bath  heated  by  a  hot  plate.  The  drainage  port  was  removed  and  a  thermocouple 
was  inserted  to  record  the  temperature  of  the  aqueous  solution.  We  investigated  the 
effect  of  temperature  on  the  peak  intensity  of  Ca(I)  (422.673  11m)  over  the  range  27 
-  99°C.  Once  the  solution  reached  99°C,  the  hot  plate  was  turned  off  and  allowed 
to  cool.  Spectra  were  taken  repeatedly  using  single-pulse  LIBS  as  the  temperature 
dropped.  Ca  line  intensities  were  measured  for  a  solution  of  1,000  ppm  Ca  and  2,540 
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ppm  NaCl  at  atmospheric  pressure  using  a  single  laser  pulse  of  37  mJ,  (4  =  100  ns, 
4  =  1  /.is,  W  =  35  /rm),  and  the  optical  configuration  shown  in  Figure  2-2(c).  Over 
this  range,  temperature  had  no  effect  on  peak  intensity. 

2.5  Conclusions 

An  optimal  range  of  low  laser  pulse  energies  exists  for  the  detection  of  Li,  Ca,  Mn,  K, 
and  Na  in  bulk  aqueous  solutions  at  both  low  and  high  pressures.  No  pressure  effect 
was  seen  on  the  emission  intensity  for  Ca  and  Na,  and  an  increase  in  intensity  with 
increased  pressure  was  seen  for  Mn.  No  line  broadening  due  to  pressure  was  observed 
for  Ca,  Na,  or  Mn  emission.  A  low  energy  pulse  may  create  a  smaller,  more  tightly 
focused  plasma  that  forms  only  at  the  focal  spot.  However,  for  a  high  energy  pulse, 
the  high  energy  density  may  cause  breakdown  even  before  the  pulse  reaches  the  focal 
spot.  This  may  allow  breakdown  to  occur  over  a  longer  distance.  In  addition,  plasma 
shielding  may  occur.  Further  studies  using  imaging  techniques  will  help  to  elucidate 
the  relationship  between  the  laser  pulse  energy  and  the  subsequent  plasma  that  is 
formed.  Using  the  dual  pulse  technique  for  several  analytes,  a  very  short  interpulse 
delay  resulted  in  the  greatest  emission  intensity.  Since  this  condition  approaches  sin¬ 
gle  pulse  conditions,  dual  pulse  LIBS  may  not  be  advantageous  for  some  elements  at 
high  pressure.  For  different  gate  delays  at  fixed  pressure,  laser  pulse  energy  affects 
peak  intensity.  The  addition  of  NaCl  enhanced  the  emission  intensity  for  Ca,  but  had 
no  effect  on  the  intensity  of  Mn  or  K  peaks.  Ca  was  detectable  over  a  wide  range 
of  concentrations  and  pressures.  In  addition,  temperature  changes  below  99°C  had 
no  noticeable  effect  on  the  emission  intensity  of  Ca.  Overall,  increased  pressure,  the 
addition  of  NaCl  to  a  solution,  and  temperature  did  not  inhibit  detection  of  analytes 
in  solution.  The  results  presented  here  suggest  that  LIBS  is  a  viable  technique  for 
in  situ  chemical  analysis  in  the  deep  ocean  and  further  work  should  be  carried  out  to 
develop  LIBS  into  an  in  situ  oceanographic  sensor. 
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Chapter  3 

Analysis  of  laser-induced 
breakdown  spectroscopy  (LIBS) 
spectra:  The  case  for  extreme 
value  statistics 


The  work  in  this  chapter  is  currently  in  press  in  Spectrochimica  Acta  Part  D  as,  A. 
P.  M.  Michel  and  A.  D.  Chave,  “Analysis  of  laser-induced  breakdown  spectroscopy 
(LIBS)  spectra:  The  case  for  extreme  value  statistics.” 


3.1  Abstract 

In  most  instances,  laser-induced  breakdown  spectroscopy  (LIBS)  spectra  are  obtained 
through  analog  accumulation  of  multiple  shots  in  the  spectrometer  CCD.  The  aver¬ 
age  acquired  in  the  CCD  at  a  given  wavelength  is  assumed  to  be  a  good  represen¬ 
tation  of  the  population  mean,  which  in  turn  is  implicitly  regarded  to  be  the  best 
estimator  for  the  central  value  of  the  distribution  of  the  spectrum  at  the  same  wave¬ 
length.  Multiple  analog  accumulated  spectra  are  taken  and  then  in  turn  averaged 
wavelength-by-wavelength  to  represent  the  final  spectrum.  In  this  paper,  the  statis¬ 
tics  of  single-shot  and  analog  accumulated  LIBS  spectra  of  both  solids  and  liquids 
were  examined  to  evaluate  whether  the  typically  used  spectrum  averaging  approach 
is  statistically  defensible.  At  a  given  wavelength,  LIBS  spectra  are  typically  drawn 
from  a  Frechet  extreme  value  distribution,  and  hence  the  mean  of  an  ensemble  of 
LIBS  spectra  is  not  necessarily  an  optimal  summary  statistic.  Under  circumstances 
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that  are  broadly  general,  the  sample  mean  for  LIBS  data  is  statistically  inconsistent 
and  the  central  limit  theorem  does  not  apply.  This  result,  appears  to  be  due  to  very 
high  shot-to-shot  plasma  variability  in  which  a  very  small  number  of  spectra  are  high 
in  intensity  while  the  majority  are  very  weak,  yielding  the  extreme  value  form  of  the 
distribution.  The  extreme  value  behavior  persists  when  individual  shots  are  analog 
accumulated.  An  optimal  estimator  in  a  well-defined  sense  for  the  spectral  average 
at  a  given  wavelength  follows  from  the  maximum  likelihood  method  for  the  extreme 
value  distribution.  Example  spectra  taken  with  both  an  Echelle  and  a  Czerny-Turner 
spectrometer  are  processed  with  this  scheme  to  create  smooth,  high  signal-to-noise 
summary  spectra.  Plasma  imaging  was  used  in  an  attempt  to  visually  understand 
the  observed  variability  and  to  validate  the  use  of  extreme  value  statistics.  The  data 
processing  approach  presented  in  this  paper  is  statistically  reliable  and  should  be  used 
for  accurate  comparisons  of  LIBS  spectra  instead  of  arithmetic  averaging  on  either 
complete  or  censored  data  sets. 


3.2  Introduction 

Laser-induced  breakdown  spectroscopy  (LIBS)  is  a  spectrochemical  technique  that 
has  been  successfully  used  for  elemental  analysis  of  solids,  liquids,  gases,  and  aerosols, 
and  is  finding  increasing  application  in  basic  and  applied  research.  However,  one  of 
the  major  problems  that  precludes  more  quantitative  use  is  a  lack  of  reproducibility 
of  spectra  at  a  given  wavelength  on  a  shot-to-shot.  basis. 

Dramatic  peak  intensity  fluctuations  at  the  shot-to-shot  level,  suggesting  the  pres¬ 
ence  of  high  random  variability,  has  been  noted  by  many  LIBS  researchers  [1  15]. 
There  are  numerous  potential  causes  for  this,  including  repetitive  laser  pulse  insta¬ 
bility,  unstable  laser  pulse  characteristics,  laser  pulse-plasma  interactions,  lens-to- 
sample  distance  variation  (which  in  turn  changes  the  distance  from  the  plasma  to  the 
collection  fiber),  laser- material  coupling,  variable  sample  ablation,  plasma  position 
instability,  matrix  effects,  perturbations  of  the  plasma  due  to  physical  and  chemical 
characteristics  of  the  sample  (i.e.,  composition,  homogeneity,  roughness,  color,  and 
moisture  content),  scattering  of  light,  atmospheric  conditions,  weak  ionization  of  the 
plasma,  and  non-optimal  collection  of  plasma  emission  [1,  3-11,  16].  Optical  insta¬ 
bility  affects  the  ablation  process,  the  plasma  profile,  and  the  plasma  volume  [7]. 
Carranza  and  Hahn  [16]  suggest  that  above  a  threshold  value,  absorption  of  pulse  en¬ 
ergy  by  the  plasma  saturates,  reducing  variability  at  higher  levels.  A  given  plasma  is 
not  completely  homogeneous,  and  there  may  be  property  gradients  due  to  boundary 


58 


effects  and  its  transient  nature.  Spatial  variation  in  the  position  of  the  plasma  changes 
the  coupling  of  the  plasma  light  into  the  collection  optics.  Carranza  and  Hahn  fur¬ 
ther  suggest  that  shot-to-shot  variation  may  be  reduced  by  using  sufficient  laser  pulse 
energy  to  achieve  saturation  and  a  suitable  collection  geometry  (backscatter  mode) 
to  minimize  spatial  variability.  The  intensity  of  the  laser  itself  can  fluctuate  by  1-5% 
[8].  However,  Castle  et  al.  simultaneously  measured  the  analyte  signal  and  the  laser 
pulse  energy,  and  found  no  significant  correlation,  suggesting  laser  pulse  variance  has 
only  a  minor  influence  on  overall  variability  [13]. 

The  sample  type  also  influences  variability.  For  aqueous  samples,  additional  fluc¬ 
tuation  can  be  caused  by  “moving  breakdown”  that  changes  the  distance  between 
the  spark  and  the  collection  fiber  as  the  plasma  moves  in  the  solution.  The  plasma 
typically  expands  along  the  beam  path  toward  the  laser,  inducing  elongate  plasmas 
that  cavitate  radially  [17].  Variability  in  aqueous  solutions  can  also  be  caused  by 
suspended  ablated  particles  [2,  18]  and  bubbles  [19]  both  by  reducing  the  energy 
delivered  to  a  sample  and  the  light  transmitted  to  the  collection  optics.  Bubbles, 
formed  when  breakdown  occurs,  and  dissolved  gases  can  scatter  or  absorb  incident 
laser  radiation  [17,  19]. 

Lazic  et  al.  [2]  report  high  variability  of  the  plasma  intensity  for  both  aqueous 
solutions  (including  bulk  water)  and  solid  samples  immersed  in  water.  Significant 
variability  was  not  observed  when  high  laser  pulse  energies  were  used  to  measure 
the  elemental  composition  of  flat  homogeneous  solid  samples.  In  aqueous  solution, 
LIBS  emission  was  sometimes  not  detectable  even  when  the  maximum  laser  energy 
was  used.  The  lack  of  emission  was  also  observed  for  rough  inhomogeneous  solids. 
For  flat  samples,  the  only  time  no  breakdown  occurred  was  when  low  laser  energies 
were  used;  yet,  shot-to-shot  signal  variability  was  always  present.  Lazic  et  al.  [10] 
reported  peak  intensity  histograms.  The  distributions  of  these  data  sets  are  clearly 
not  Gaussian,  and  show  that  a  very  high  intensity  peak  is  a  rare  event,  with  very  low 
intensity  occurring  for  the  bulk  of  the  trials. 

For  a  solid  sample,  inhomogeneity,  porosity,  or  surface  roughness  can  change  the 
distance  between  the  focusing  optics  and  the  sample,  either  from  prior  crater  for¬ 
mation  or  by  changing  the  location  of  ablation.  Panne  et  al.  [12]  report  significant 
pulse-to- pulse  variation  of  the  plasma  electronic  excitation  temperature  and  electron 
density  from  material-laser  interaction  for  homogeneous  glass  samples. 

Laser  ablation  is  highly  nonlinear,  and  even  more  so  in  aerosol  samples  as  the 
plasma  may  form  at  different  positions  along  the  beam  [8].  For  aerosols,  the  location 
of  the  particles  within  the  plasma  volume  and  the  focal  volume  of  the  optics  con- 
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tribute  to  variability  [14].  Schechter’s  [15]  analysis  of  spectral  fluctuations  of  aerosols 
showed  large  shot-to-shot  variability  possibly  caused  by  laser  pulses  hitting  different, 
numbers  of  particles,  particle  characteristic  variation  (size,  mass,  and  location),  and 
location  variation  of  the  plasma.  Whatever  the  cause  of  observed  LIBS  intensity 
variability,  analog  averaging  of  multiple  plasma  emissions,  where  light  from  numerous 
laser  shots  is  accumulated  on  a  CCD  to  create  a  single  spectrum,  is  an  often  used  ex¬ 
perimental  approach  in  order  to  increase  the  signal  and  the  signal-to-noise  ratio  in  the 
presence  of  shot-to-shot  variability  [20].  Analog  averaged  spectra  are  replicated  and 
the  ensemble  of  replicates  are  in  turn  averaged  to  create  a  representative  spectrum. 
However,  this  implicitly  assumes  that  the  sample  mean  is  a  reasonable  estimator  for 
the  statistical  average,  and  this  condition  may  not  hold  for  some  non-Gaussian  dis¬ 
tributions. 

LIBS  researchers  have  recognized  the  potential  impact  of  intensity  variability, 
and  have  devised  a  variety  of  methods  to  reduce  the  effect  of  spectral  variability. 
Schechter  used  a  rejection  algorithm  to  eliminate  anomalous  spectra  (e.g.,  spectra 
with  no  elemental  lines,  spectra  with  a  too  weak  or  too  intense  baseline  due  to  laser 
fluctuations,  and  weak  spectra)  from  the  ensemble.  This  typically  removed  75%  of 
measured  spectra  [15].  Carranza  and  Hahn  [16]  used  a  sorting  algorithm  to  remove 
irregular  spectra,  eliminating  60-70%  of  single  shot  data.  Lazic  et  al.  [2]  removed 
spectra  below  a  threshold  value  to  increase  the  signal-to-noise  ratio  and  make  ('mis¬ 
sion  lines  more  readily  visible. 

In  the  present  work,  rather  than  using  an  ad  hoc  approach,  the  statistical  vari¬ 
ability  of  LIBS  spectral  intensity  has  been  quantified  and  a  data  processing  scheme 
based  on  the  observed  statistics  has  been  devised.  It  will  be  demonstrated  that.  LIBS 
intensity  (whether  single  shot  or  analog  averaged)  typically  has  a  Frechet  extreme 
value  distribution,  and  that,  for  the  characteristic  range  of  statistical  parameters,  the 
distribution  may  not  possess  a  variance.  As  a  consequence,  the  sample  mean  is  not 
an  appropriate  estimator  for  the  average  intensity,  the  central  limit  theorem  does  not 
apply,  and  Gaussian-based  inference  will  be  in  error.  Instead,  a  maximum  likelihood 
estimator  for  the  extreme  value  distribution  is  advocated  as  an  alternative.  The  result 
is  illustrated  using  single  shot,  and  analog  averaged  LIBS  spectra  for  a  solid  target, 
using  one  experimental  set-up  and  for  bulk  aqueous  solutions  using  two  experimental 
set-ups. 
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Figure  3-1:  Laboratory  set-up  of  LIBS  using  an  Echelle  spectrometer. 

3.3  Experimental 

3.3.1  Echelle  Spectrometer  Set-up 

To  examine  the  variability  of  LIBS  at  shot-to-shot  and  analog  averaged  levels,  the 
peak  intensity  was  examined  over  the  580  -  600  nin  range  for  Na  using  both  solid 
NaCl  (halite)  obtained  from  Fisher  Scientific  and  an  aqueous  NaCl  solution.  Solutions 
were  made  using  de-ionized  water  and  NaCl  to  yield  a  Na  concentration  of  100  parts 
per  million  (ppm,  wt./vol.).  The  variability  of  intensity  was  measured  for  the  solid 
sample  using  an  Echelle  spectrometer  and  for  aqueous  solutions  using  both  Echelle 
and  Czerny- Turner  units.  A  dark  background  spectrum  was  initially  subtracted  from 
all  raw  spectra.  For  the  halite  specimen,  100  single  shot  and  10  shot  analog  averaged 
spectra  were  obtained.  For  the  aqueous  specimen,  100  single  shot  and  100  shot  analog 
averaged  spectra  were  collected. 

The  first  experimental  set-up  utilizing  an  Echelle  spectrometer  (LLA  Echelle  ESA 
3000)  is  shown  in  Figure  3-1.  The  spectrometer  is  capable  of  detecting  elements  over 
the  200  -  780  nm  range  with  a  spectral  resolution  of  10  to  50  pm.  A  Big  Sky  CFR- 
200  Nd:YAG  laser  (7.5-ns  pulse  width)  operated  at  the  fundamental  wavelength  of 
1064  nm  with  a  repetition  rate  of  5  Hz  was  used  for  plasma  excitation.  The  laser  is 
equipped  with  a  variable  attenuator  controlled  by  a  computer  that  allows  laser  pulse 
energy  to  range  from  0  to  200  mJ  in  increments  of  <1  mJ.  A  timing  box  (Berkeley 
Nucleonics  Corporation  Model  565)  was  used  to  accurately  control  firing  of  the  laser 
in  relation  to  turn-on  of  the  spectrometer. 

For  liquid  samples,  a  cubic  titanium  sample  chamber  (8.89  cm  x  8.89  cm  x  8.89 
cm)  equipped  with  two  sapphire  windows  (Meller  Optics  -  2.54  cm  diameter  x  0.64 
cm  thick,  AR  coated  at  1064  nm,  custom  part)  that  allows  laser  pulses  to  enter  the 
cell  and  the  plasma  to  be  imaged  from  the  side  of  the  cell  (orthogonal  to  the  entering 
laser  beam)  was  used.  AR-coated  optics  focus  the  laser  beam  into  the  chamber.  For 
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solid  samples,  the  chamber  was  removed  and  the  final  focusing  lens  was  placed  in 
front  of  the  sample.  Additional  optics  were  used  to  focus  the  plasma  light  onto  an 
optical  fiber  that  delivers  it  to  the  spectrometer.  Data  were  collected  using  ESAWIN 
software.  All  spectra  were  taken  with  a  pulse  energy  of  80  mJ  and  the  maximum 
MCP  amplification  of  4000.  For  aqueous  NaCl  solutions,  the  gate  delay  =  75  ns  and 
gate  width  =  200  ns.  For  halite  samples,  the  gate  delay  =  10,000  ns  and  gate  width 
=  100  ns. 

Plasma  images  were  taken  using  a  Pixelfiy  camera  with  a  microscope  lens  and 
an  iris  diaphragm.  The  images  were  taken  through  the  sapphire  window  on  the 
pressure  chamber,  orthogonal  to  the  incoming  laser  pulses  (80  mJ/pulse).  The  shutter 
remained  open  for  5  /xs  and  was  externally  synched  to  the  Q-switch  of  the  laser. 


3.3.2  Czcrny-Turner  Spectrometer  Set-up 


The  second  set-up  used  a  Czerny-Turner  spectrometer  and  is  shown  in  Figure  3-2.  A 
Continuum  Surelite  III  laser  (5-ns  pulse  width,  1004  nm,  1  Hz  repetition  rate)  was 
used  for  plasma  excitation  with  a  pulse  energy  of  81  mJ.  Laser  pulses  were  focused 
into  a  chamber  constructed  of  stainless  steel  Swagelok  fittings  with  six  2.54  cm-ID 
and  3.18  cm-OD  ports.  Two  ports  were  fitted  with  2.54  cm  diameter,  0.32  cm  thick 
circular  sapphire  windows  (MSW100/125,  Meller  Optics  Incorporated)  held  in  place 
by  hex  nuts  and  sealed  with  rubber  washers,  allowing  1.91  cm  of  each  window  to 
be  visible  outside  the  cell.  The  plasma  emission  was  focused  onto  a  2-mm-core- 
diameter,  0.51-N.A.  light  guide  (Edmund  Scientific  Co.  Model  02551).  The  light 
guide  was  connected  to  a  0.25-m,  //4  spectrograph  (Chromex  model  250is/RF)  with 
a  1200-groove/mm  grating  blazed  at  500  nm.  Data  were  collected  on  an  intensified 
CCD  detector  (Princeton  Instruments,  I-Max  1024E)  and  acquired  with  a  computer 
running  WinSpec/32  software. 

All  spectra  were  taken  at  the  maximum  gain  setting  of  255,  with  gate  delay  =  175 
ns  and  gate  width  =  200  ns.  Solutions  were  made  using  de-ionized  water,  M11SO4, 
ZnBr2,  and  NaSO.i.  The  solution  contained  5000  parts  per  million  (ppm,  wt./vol.) 
Mn,  5000  ppm  Zn,  and  2000  ppm  Na. 
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Figure  3-2:  Laboratory  set-up  of  LIBS  using  a  Czerny-Turner  spectrometer. 

3.4  Results  and  Discussion 

3.4.1  The  Generalized  Extreme  Value  Distribution 

Extreme  value  distributions  describe  the  stochastic  behavior  of  the  maximum  or  min¬ 
imum  of  independent  and  identically  distributed  random  variables  drawn  from  some 
parent  distribution.  There  are  three  types  of  extreme  value  distributions:  Weibull. 
Gumbel,  and  Frechet.  The  von  Mises-Jenkinson  or  generalized  extreme  value  dis¬ 
tribution  (GEVD)  combines  the  three  into  a  single  functional  form  [21,  22].  The 
probability  density  function  for  the  GEVD  is  given  by 


f(x\k,fi,cr)  =  (-)exp(-(l  +  k— — — )  *)(1  +  k— — — )(  1  (3.1) 

a  o  o 

where  the  distribution  is  Weibull,  Gumbel  and  Frechet  for  the  shape  parameter  k< 0, 
A;=0,  and  A:>0,  respectively.  The  remaining  parameters  in  Equation  3.1  are  the  loca¬ 
tion  parameter  /z  (analogous  to  the  mean  for  the  Gaussian  distribution)  and  the  scale 
parameter  a  (analogous  to  the  standard  deviation  for  the  Gaussian  distribution).  The 
Gumbel  distribution  is  obtained  through  an  appropriate  limiting  process,  but  is  not 
of  interest  in  the  present  work.  For  the  Weibull  and  Frechet  distributions,  the  range 
of  the  variate  is  respectively  -oo<x<fi-a/k  and  [i-a  jk<x< oo.  The  Weibull  distribu¬ 
tion  has  a  finite  upper  endpoint  and  hence  corresponds  to  a  short-tailed  parent.  The 
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Frechet  distribution  has  a  polynomially  decreasing  upper  tail  and  corresponds  to  a 
long-tailed  parent. 

As  will  subsequently  be  shown,  LIBS  intensities  are  typically  distributed  as  Frechet 
extreme  value.  It  is  important  to  note  that  the  second  and  higher  order  moments 
(and  hence  the  variance)  do  not  exist  for  A; >  1  / 2 ,  and  the  first  moment  (the  mean) 
does  not  exist  for  A:>1.  This  implies  that  the  standard  estimator  for  the  ensemble 
average,  the  sample  mean,  will  either  be  inconsistent  (i.e.,  will  not  display  a  reduced 
variance  as  the  size  of  the  sample  increases)  for  A:>  1/2  or  will  not  exist  at  all  for  k>  1. 
A  different  formulation  is  required  to  obtain  the  three  parameters  in  Equation  3.1  so 
that  defensible  statistical  inferences  about  LIBS  intensities  can  be  made.  A  standard 
approach  is  the  maximum  likelihood  method  that  seeks  the  solutions  for  k,  ft  and  o 
that  maximize  the  joint  distribution  of  a  given  set  of  data,  or  the  likelihood  function. 
For  independent  samples,  the  joint  distribution  of  N  data  is  the  product  of  Equation 
3.1  for  each  datum  with  common  shape,  location  and  scale  parameters 


L{k,n,a\xt)  =  TT(  — )ea:p(— (1  +  k 

A  A  G 


X, 
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)  *  )(1  +  k 
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■I—, 


(3.2) 


This  is  maximized  by  setting  the  first  derivative  of  Equation  3.2,  or  its  logarithm,  for 
each  parameter  to  zero,  yielding  a  set  of  three  equations 


-Y -(1  +  *:(*»-/*))-(»u)  + 
a 

X  1 


k  +  1 


1  +  kx-^ 

(7 


=  0 


(3.3) 


p£S>g(i  +  *^^)-*«i  +  *o 

i  l 

kJ2( i  +  ^))-i('og(1+^)  - 


N 

“<r(  1  +  ) 


)  + 


ka(  \  { 


l£gbK)]=0 


(3.4) 


N 


--  +  EH!£t£ 

g  L ^  a1 

i— 1 


)(1 + *  -  (1  + 1  )(T-2_r)(^  -  *,)]  =  o 

x  (7  ' 

(3.5) 


These  are  coupled  and  nonlinear,  and  must  be  solved  iteratively  for  the  maximum 
likelihood  estimators  (mlcs)  k ,  fi  and  a  using  an  appropriate  algorithm.  For  the 
GEVD,  the  rules  are  asymptotically  efficient  (loosely  speaking,  highly  concentrated 
about  the  true  value  for  large  numbers  of  data),  normal  and  unbiased,  but  are  nci- 
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ther  unbiased  nor  fully  efficient  for  finite  samples.  However,  neither  is  the  sample 
mean  when  the  second  and  higher  moments  do  not  exist.  The  mles  have  the  distinct 
advantage  of  being  defined  when  the  moments  of  the  distribution  do  not  exist,  and 
are  relatively  easy  to  compute.  Presuming  that  the  extreme  value  distribution  is  a 
good  fit  to  LIBS  intensity  data,  the  maximum  likelihood  estimate  for  the  location 
parameter  /t  is  a  good  representation  of  the  peak  intensity  at  a  given  wavelength,  and 
the  scale  parameter  a  and  asymptotic  normality  can  be  used  to  compute  approximate 
confidence  intervals  on  /}. 

3.4.2  Applicability  of  Extreme  Value  Statistics 

Quantile-quantile  (q-q)  plots  will  be  used  to  demonstrate  that  LIBS  intensities  are 
typically  distributed  as  the  Frechet  extreme  value  distribution.  The  N  quantiles  of 
a  target  distribution  are  the  abscissa  values  that  divide  the  area  under  the  pdf  into 
N+l  equal  probability  intervals.  They  are  easily  obtained  from  the  pdf  by  solving 
for  Qj  in 


(3.6) 


where  j=l,...,N.  The  order  statistics  of  the  intensity  data  are  obtained  by  sorting 
them  into  ascending  order.  The  order  statistics  divide  the  area  under  the  target  pdf 
into  intervals  that  will  correspond  to  equal  probability  if  the  data  are  drawn  from 
it,  and  hence  a  plot  of  the  quantiles  against  the  order  statistics  will  be  a  straight 
line.  Systematic  departures  of  the  data  from  the  distribution  are  visible  as  changes 
in  slope,  and  anomalous  values  or  outliers  are  manifest  at  the  extremes,  and  hence  a 
q-q  plot  is  a  useful  qualitative  tool  to  assess  the  suitability  of  the  target  distribution 
as  a  statistical  model. 

The  fit  may  be  quantified  by  testing  the  null  hypothesis  that  the  data  are  drawn 
from  the  target  distribution  against  the  alternate  hypothesis  that  they  are  not  using 
the  nonparametric  Kolmogorov-Smirnov  statistic  [23]  that  compares  the  empirical 
and  target  cumulative  distribution  functions.  The  Komogorov-Smirnov  test  statistic 
may  be  assessed  at  the  standard  0.95  level  for  which  the  critical  value  is  0.134  for  100 
realizations,  as  were  used  throughout  this  work.  The  null  hypothesis  is  rejected  if  the 
test  statistic  exceeds  the  critical  value. 

Figure  3-3  shows  q-q  plots  for  the  588.9953  nm  Na  I  peak  for  single  shot  and  10 
shot  analog  accumulations  on  halite  using  the  Echelle  set-up.  Both  are  approximately 
straight,  and  both  accept  the  null  hypothesis  that  the  extreme  value  distribution  is 
correct  (Kolmogorov-Smirnov  test  statistics  of  0.067  and  0.066,  respectively).  Figure 
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3-4  shows  q-q  plots  at  the  same  wavelength  for  single  shot  and  100  shot  analog 
accumulations  for  Na  in  bulk  aqueous  solution  using  the  Echelle  set-up.  The  single 
shot  data  exhibit  a  shift  in  slope  that  reflects  the  fact  that  the  vast  majority  of  the 
data  correspond  to  no  signal,  and  hence  are  only  instrument  noise.  The  100  shot 
analog  accumulations  are  slightly  short-tailed  at  the  upper  end.  Nevertheless,  both 
pass  the  Kolmogorov-Smirnov  test  with  test  statistics  of  0.113  and  0.098,  respectively. 
Figure  3-5  compares  q-q  plots  at  the  same  peak  for  single  shot  and  100  shot  analog 
accumulations  using  the  Czerny-Turner  set  up  with  a  bulk  aqueous  target.  Both  are 
slightly  short  tailed  at  the  top  of  the  distribution,  but  both  accept  the  null  hypothesis 
that  the  data  are  extreme  value  (Kolmogorov-Smirnov  test  statistics  of  0.084  and 
0.058,  respectively).  Similar  results  are  observed  for  the  589.5923  nm  Na  peak  (not 
shown),  or  at  other  wavelengths  where  signal  is  present.  Further,  the  shape  parameter 
persistently  lies  in  the  region  corresponding  to  the  Frechet.  extreme  value  distribution, 
and  in  many  instances  exceeds  0.5  so  that  the  variance  does  not  exist. 
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(b)  10  Shot  Accumulations 


Figure  3-3:  q-q  plots  for  the  588.9953  nm  Na  I  peak  for  halite  using  the  Echelle  set-up 
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Figure  3-4:  q-q  plots  for  the  588.9953  nm  Na  I  peak  for  bulk  aqueous  solution  using 
the  Echelle  set-up 
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Figure  3-5:  q-q  plots  for  the  588.9953  nm  Na  I  peak  for  bulk  aqueous  solution  using 
the  Czerny-Turner  set-up 
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3.4.3  Extreme  Value  Statistical  Parameters 


As  a  demonstration  of  the  importance  of  using  an  appropriate  set  of  statistical  esti¬ 
mators  for  LIBS  intensity  data,  the  extreme  value  mles  and  the  sample  mean  were 
computed  over  the  wavelength  band  578.5  -  605.9  nm  for  the  halite  and  aqueous  Na 
samples  using  both  the  Eelielle  and  Czerny-Turner  set-ups.  At  each  wavelength,  100 
realizations  of  10  (halite)  or  100  (aqueous  samples)  analog-accumulated  shots  serve 
as  the  data  from  which  the  mles  (Equations  3.3  -  3.5)  and  the  sample  mean  are 
estimated. 

Figure  3-6  compares  the  shape  parameter,  the  location  parameter,  and  the  sample 
mean  for  halite  using  the  Echelle  set-up.  While  there  is  wavelength- by- wavelength 
statistical  variability,  the  shape  parameter  is  persistently  above  the  0.5  threshold 
beyond  which  the  variance  does  not  exist,  and  frequently  exceeds  the  1.0  threshold 
beyond  which  the  mean  does  not  exist,  especially  in  the  vicinity  of  the  589  nm  Na 
doublet  and  over  two  bands  slightly  above  600  nm.  The  extreme  value  location  param¬ 
eter  produces  a  smooth  representation  of  the  LIBS  spectrum  with  limited  statistical 
variability  that  is  consistent  with  the  number  of  samples.  By  contrast,  the  sample 
mean  displays  substantial  statistical  variability  with  two  anomalous  nulls  amid  the  Na 
doublet  band.  These  correspond  to  the  wavelengths  where  the  shape  parameter  dips 
below  1.0,  so  that  the  mean  becomes  defined.  For  data  that  are  long-tailed,  such  as 
those  drawn  from  a  Frechet  extreme  value  distribution,  the  sample  mean  will  be  dom¬ 
inated  by  a  few  large  values.  This  results  both  in  the  large  wavelength-to-wavelength 
variability  that  is  apparent  in  the  sample  mean  and  the  substantial  difference  in  am¬ 
plitude  between  the  sample  mean  and  the  location  parameter.  In  addition,  two  large 
peaks  are  observed  above  600  nm  that  are  barely  visible  in  the  extreme  value  location 
parameter.  These  correspond  to  peaks  well  above  1.0  in  the  shape  parameter  where 
the  extreme  value  distribution  is  very  long-tailed,  and  serve  as  graphic  illustration  of 
the  sort  of  erroneous  conclusions  that  can  be  derived  through  use  of  inappropriate1 
statistical  estimators. 

Figure  3-7  compares  the  shape  parameter,  the  location  parameter,  and  the  sample 
mean  for  Na  in  bulk  aqueous  solution  using  the  Echelle  set-up.  The  shape  parameter 
is  much  more  uniform  with  wavelength  than  for  halite  (Figure  3-6),  but  persistent  ly 
lies  around  0.75  where  the  variance  does  not  exist.  As  a  consequence,  the  extreme 
value  location  parameter  displays  much  less  variability  than  the  sample  mean,  as  in 
Figure  3-6.  Figure  3-7  is  another  example  where  the  use  of  a  standard  sample  mean 
estimator  may  lead  to  incorrect  inferences. 

Figure  3-8  compares  the  shape  parameter,  the  location  parameter,  and  the  sample 
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mean  for  Na  in  bulk  aqueous  solution  using  the  Czerny-Turner  set-up.  In  contrast 
to  the  results  with  the  Echelle  set-up,  the  shape  parameter  is  persistently  below  0.5, 
and  displays  occasional  excursions  below  0  where  the  distribution  is  Weibull.  Since 
both  the  mean  and  variance  exist  throughout  the  wavelength  domain,  the  extreme 
value  location  parameter  and  the  sample  mean  yield  qualitatively  similar  results. 
However,  the  distribution  remains  extreme  value  rather  than  simple  Gaussian,  and 
the  uncertainty  inferred  for  the  sample  mean  will  be  inaccurate. 

While  LIBS  intensity  data  (whether  single  shot  or  analog  accumulated)  empiri¬ 
cally  appear  to  persistently  be  drawn  from  an  extreme  value  distribution,  systematic 
differences  are  observed  between  different  experimental  set-ups.  Whether  this  is  due 
to  the  spectrometer  or  laser  design,  the  experimental  geometry,  the  element  under 
study  or  some  other  factor  remains  unknown.  Since  incorrect  conclusions  might  be 
drawn  from  the  use  of  an  inappropriate  statistical  model,  it  is  strongly  urged  that 
LIBS  practitioners  examine  their  data  to  determine  the  correct  approach  for  each 
set-up  and  sample. 
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Figure  3-6:  Shape  parameter,  location  parameter,  and  sample  mean  for  halite  using 
the  Echelle  set-up 
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Figure  3-7:  Shape  parameter,  location  parameter,  and  sample  mean  for  bulk  aqueous 
solution  using  the  Echelle  set-up 
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Figure  3-8:  Sliape  parameter,  location  parameter,  and  sample  mean  for  bulk  aqueous 
solution  using  the  Czerny-Turner  set-up 
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3.4.4  Variability 


In  an  attempt  to  understand  the  source  of  data  variability,  images  of  plasmas  in  liquids 
were  taken  and  the  energies  of  numerous  single  shot  laser  pulses  were  compared. 
Figure  3-9  shows  three  sample  plasma  images  taken  under  the  same  conditions  that 
visually  demonstrates  the  significant  variability  of  plasma  formation.  The  images  are 
time  averaged  over  5  ps.  Although  hundreds  of  plasma  images  were  taken,  these 
images  were  selected  to  illustrate  the  variation  in  intensity,  size,  and  location  of 
formation  of  the  plasma  within  a  bulk  liquid.  The  plasma  shown  in  Figure  3-9(c) 
displays  significantly  greater  emission  intensity  and  size  than  those  in  (a)  and  (b), 
illustrating  the  extreme  nature  of  some  plasmas.  In  contrast,  the  plasma  in  Figure  3- 
9(b)  is  very  weak,  with  very  little  emission  produced.  The  distinct  differences  shown 
between  these  three  plasmas  is  a  clear  indication  that  on  a  shot-to-shot  basis  that 
plasma  variability  exists  and  that  extreme  plasmas  are  formed  which  could  account 
for  the  extreme  intensities  recorded. 

The  energy  of  500  laser  pulses  was  measured  to  examine  the  contribution  of  laser 
pulse  energy  fluctuations  to  plasma  variability.  Figure  3-10  shows  that  shot-to-shot 
pulse  energy  fluctuations  do  exist;  however,  the  variability  is  not  extreme  in  nature, 
suggesting  that  laser  pulse  energy  variation  is  not  the  dominant  cause  of  plasma  and 
peak  intensity  variations.  This  suggests  that  the  variability  of  the  plasma  formation 
is  due  to  other  effects. 
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Figure  3-9:  Images  of  plasmas  formed  in  bulk  aqueous  solution  that  illustrate  the 
shot-to-shot  variability  of  formation.  Images  were  taken  orthogonal  to  the  incoming 
laser  beam.  In  the  images  shown,  the  beam  enters  from  the  left.  Plasmas  were  formed 
using  80  mJ  of  laser  pulse  energy. 
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Figure  3-10:  Comparison  of  laser  energies  measured  for  500  individual  laser  pulses. 
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3.5  Conclusions 


Examination  of  the  variability  of  peak  intensity  for  both  single  shot  and  analog  ac¬ 
cumulated  LIDS  spectra  reveals  that  such  data  are  drawn  from  an  extreme  value 
distribution.  In  many  instances,  the  distribution  has  no  variance,  and  in  some  cases 
the  mean  is  also  undefined.  Under  either  circumstance,  the  use  of  the  sample  mean  or 
variants  that  include  censoring  will  be  statistically  inconsistent  and  the  central  limit 
theorem  will  not  apply.  A  maximum  likelihood  estimator  data  processing  scheme  is 
presented  that  accurately  deals  with  the  extreme  value  nature  of  laser-induced  plasma 
formation.  It  is  strongly  urged  that  this  approach  be  used  to  ensure  accurate  scien¬ 
tific  inference  from  LIBS  data,  and  that  use  of  estimators  based  on  the  sample  mean 
be  discontinued.  Plasma  images  reveal  large  spatial  and  intensity  differences  on  a 
shot-to-shot  basis.  Laser  pulse  energy  fluctuations  are  shown  to  contribute  to  the 
variability  but  are  not  the  primary  source. 
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Chapter  4 


Single  pulse  laser-induced 
breakdown  spectroscopy  of  bulk 
aqueous  solutions  at  oceanic 
pressures:  Interrelationship  of  gate 
delay  and  pulse  energy 

4.1  Abstract 

The  ability  to  make  sustained  measurements  of  ocean  processes  is  limited  by  the 
number  of  sensors  that  are  usable  for  long-term  in  situ  analysis.  In  recent  years, 
laser-induced  breakdown  spectroscopy  (LIBS)  has  been  identified  as  a  viable  technique 
for  development  into  an  oceanic  chemical  sensor.  In  this  paper,  single  pulsed  laser- 
induced  breakdown  spectroscopy  of  high  pressure  bulk  aqueous  solutions  is  used  to 
detect  three  analytes  (sodium,  manganese,  and  calcium)  which  are  of  key  importance 
in  hydrothermal  vent  fluids,  an  ocean  environment  that  would  greatly  benefit  from  the 
development  of  an  oceanic  LIBS  sensor.  The  interrelationship  of  the  key  experimental 
parameters,  pulse  energy  and  gate  delay,  for  a  range  of  pressures  up  to  2.76  x  10'  Pa, 
are  studied.  A  minimal  effect  of  pressure  on  the  peak  intensity  is  observed.  A  short 
gate  delay  (less  than  200  ns)  must  be  used  at  all  pressures.  The  need  for  a  relatively 
low  laser  pulse  energy  (less  than  «  60  mJ)  for  optimal  detection  of  analytes  at  high 
pressure  is  also  established.  Na,  Mn,  and  Ca  are  detectable  at  pressures  up  to  2.76  x 
107  Pa  at  50  ppm,  500  ppm,  and  50  ppm,  respectively,  using  an  Echellc  spectrometer. 
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4.2  Introduction 


New  chemical  sensors  are  needed  for  use  in  process  studies,  and  are  of  critical  im¬ 
portance  as  oceanography  shifts  to  a  new  operational  mode  using  permanent  ocean 
observatories.  New  sensors  take  significant  time  to  develop  and  transform  from  bench- 
top  laboratory  prototypes  to  ocean-going  systems.  The  development  phase  requires 
validation  that  an  analytical  technique  will  work  under  in  situ  conditions. 

LIBS  is  a  type  of  atomic  emission  spectroscopy  that  has  been  identified  as  a  viable 
technique  for  use  as  a  field-going  sensor  for  geochemical  and  environmental  sensing 
[1],  For  example,  a  new  mobile  instrument  has  been  developed  for  evaluating  polluted 
soils  [2,  3].  Palanco  et  al.  have  proposed  a  field  deployable  laser-induced  breakdown 
spectrometer  system  for  stand-off  measurements  at  hundreds  of  meters  range  [4]. 
Several  groups  have  investigated  LIBS  for  space  exploration  [5  10] .  Courreges-Lacoste 
et  al.  [7]  are  developing  a  combined  Raman/LIBS  instrument  for  investigating  past 
and  present  life  on  Mars.  Arp  et  al.  [6]  have  investigated  the  use  of  LIBS  in  the  high 
temperature  (>700K),  high  pressure  (order  of  9  x  106  Pa)  environment  of  Venus. 
Another  proposed  in  situ  application  of  LIBS  is  its  development  into  an  oceanic 
chemical  sensor. 

One  ocean  environment  that  would  benefit  greatly  from  the  development  of  such  a 
sensor  is  hydrothermal  vents  that  occur  at  mid-ocean  ridges  where  seawater  circulates 
through  the  permeable  ocean  crust.  As  seawater  moves  through  the  crust,  the  fluid 
interacts  with  the  surrounding  rock,  inducing  major  chemical  changes  to  the  rock  and 
the  fluid.  At  vent  orifices,  exit  temperatures  reach  200  -  405°C  at  ambient  pressures 
of  8.1  x  10f)  Pa  to  3.6  x  10'  Pa  corresponding  to  ocean  depths  of  800  m  to  3600 
m  [11],  In  situ  chemical  measurement  of  the  fluids  is  difficult  due  to  the  corrosive 
nature  of  the  vent  environment  and  irreversible  changes  in  composition  that  occur 
when  they  are  removed  to  the  surface.  Three  elements  of  importance  in  vent  fluids 
are  sodium,  calcium,  and  manganese.  Sodium  is  the  most  abundant  cation  in  vent 
fluids  and  can  be  studied  to  understand  phase  separation  processes  [12],  Manganese 
exists  as  a  trace  metal  in  seawater  but  is  leached  from  the  host  rock  making  it  present 
at  higher  concentrations  in  vent  fluids  [12].  When  measured  simultaneously  with  Fe, 
Mil  can  be  used  as  an  indication  of  subsurface  deposition  as  Fe  precipitates  out  while 
Mil  stays  in  solution.  Calcium  is  the  second  most  abundant  cation  in  vent  fluids,  and 
is  typically  enriched  in  vent  fluids,  when  compared  to  seawater  [13].  Ca  is  released 
into  vent  fluids  when  sodium  is  taken  up  during  albitization  reactions  with  the  host 
rock.  In  vent  fluids,  concentrations  range  from  approximately  250  -  23,163  ppm  for 


82 


Na,  0.6  -  399  ppm  for  Mn,  and  -54  -  4477  ppm  for  Ca  [11].  In  seawater,  concentrations 
are  approximately  10,933  ppm  Na,  <0.001  ppm  Mn,  and  419  ppm  Ca  [11]. 

There  is  limited  prior  work  on  the  study  of  dissolved  analytes  within  bulk  aqueous 
solutions  [14  20]  due  to  inherent  difficulties  in  detection.  The  plasma  formed  in  a 
bulk  liquid  displays  reduced  light  intensity  and  a  reduction  in  emission  lifetime  due 
to  quenching  [14,  15,  21,  22].  In  addition,  spectral  lines  are  broadened  through  the 
Stark  effect  [14].  Furthermore,  “moving  breakdown”  occurs  that  randomly  changes 
the  distance  between  the  plasma  and  the  collection  optics,  a  phenomenon  that  is  not 
important  for  solids  in  air.  The  plasma  expands  along  the  beam  path  of  the  laser, 
resulting  in  an  elongated  plasma  that  cavitates  cylindrically  [23].  For  many  aqueous 
applications,  these  issues  can  be  avoided  by  analysis  on  a  liquid  surface,  jet.  or  film; 
however,  for  the  development  of  an  in  situ  oceanic  system,  it  is  necessary  to  work 
directly  with  bulk  liquids. 

Although  we  have  previously  reported  on  the  successful  use  of  LIBS  for  detection  of 
bulk  aqueous  analytes  at  high  pressure  (up  to  2.76  x  107  Pa)  [17  19],  the  development 
of  LIBS  into  an  oceanic  chemical  sensor  requires  the  optimization  of  the  experimental 
system  to  maximize  the  signal-to-background  ratio  (SBR)  of  the  spectra  and  improve 
the  limit  of  detection.  In  this  work,  a  comprehensive  study  of  the  effect  of  the  two 
key  parameters  for  single  pulse  LIBS  on  SBR  was  completed.  Peak  intensities  were 
measured  to  determine  optimal  conditions  for  the  detection  of  Na,  Ca,  and  Mn  at  high 
pressure.  Subsequently,  calibration  curves  were  constructed  to  estimate  the  limits  of 
detection  using  an  Echelle  spectrometer. 


4.3  Experimental 

The  laboratory  set-up  for  simulating  a  LIBS  sensor  in  the  deep  ocean  is  depicted 
in  Figure  4-1.  Plasma  formation  is  induced  with  a  Big  Sky  CFR-200  Nd:YAG  laser 
operated  at  1064  nm  with  a  5  Hz  repetition  rate.  The  laser  is  equipped  with  a 
motorized  variable  attenuator,  serially  controlled  by  a  computer,  enabling  the  laser 
pulse  energy  (E)  to  be  varied  from  0  mJ  to  200  mJ  in  increments  of  approximately 
1  mJ.  Plasma  emission  is  collected  with  an  Echelle  spectrometer  (LLA  Echelle  ESA 
3000)  capable  of  detecting  wavelengths  of  200  -  780  nm  with  a  spectral  resolution  of 
10  -  50  pm.  Timing  control  of  the  laser  and  turn-on  of  the  spectrometer  is  managed 
by  a  timing  box  (Berkeley  Nucleonics  Corporation  Model  565).  The  LIBS  timing 
parameters  are  gate  delay  (4,  the  time  between  the  laser  pulse  and  turn-on  of  the 
spectrometer)  and  gate  width  (4,  the  integration  time  of  the  spectrometer),  both 
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Figure  4-1:  LIBS  laboratory  set-up 


shown  in  Figure  4-2. 

An  8.89  cm  x  8.89  cm  x  8.89  cm  titanium  sample  chamber  that  holds  27  ml 
of  liquid  is  connected  to  a  high  pressure  metering  pump  (Eldex  Model  A-30-S)  that 
is  used  to  pressurize  samples  to  4.1  x  10'  Pa.  The  sample  chamber  is  equipped 
with  a  sapphire  window  (Meller  Optics,  2.54  cm  diameter  x  0.G4  cm  thickness.  AR- 
coated  at  532  nm/1064  run,  custom  part)  that  enables  laser  pulses  to  enter  the  sample 
chamber.  A  series  of  AR-coated  optics  are  used  to  focus  the  laser  beam  into  the  sample 
chamber  (Figure  4-3),  with  the  final  focusing  lens  fit  into  the  sample  chamber.  An 
additional  lens  is  used  to  focus  the  plasma  light  onto  an  optical  fiber  and  delivers  it 
to  the  Eehelle  spectrometer  (Figure  4-3).  The  plasma  light  is  collected  collinear  to 
the  incoming  laser  beam.  This  optical  geometry  was  selected  because  it  is  the  only 
practical  configuration  for  an  ocean-going  LIBS  system.  Data  were  collected  using 
ESAWIN  software.  Laser  energy  is  measured  using  a  laser  energy  sensor  (Coherent 
J25LP-MB)  combined  with  an  energy  meter  (Coherent  FieldMaxII-Top). 

To  determine  optimal  conditions  for  the  detection  of  the  three  analytes,  spectral 
intensities  were  measured  over  a  range  of  LIBS  system  parameters.  For  Na  and  Mn, 
detailed  studies  were  conducted  at  five  pressures  (1  x  105  Pa,  6.89  x  10*’  Pa,  1.38  x 
10'  Pa,  2.07  x  10'  Pa,  and  2.76  x  10'  Pa).  For  Ca,  the  studies  were  conducted  at 
three  pressures  (1  x  105  Pa,  1.38  x  107  Pa,  and  2.76  x  107  Pa).  The  gate  delay  and 
laser  pulse  energy  were  systematically  varied  to  determine  their  effect  on  both  plasma 
intensity  and  SBR.  Five  spectra  were  taken  for  each  parameter  pair,  each  composed 
of  100  accumulated  shots.  The  laser  pulse  energy  ranged  from  10  to  170  mJ  in  10  mJ 
increments.  Laser  beam  waist  width  dCTo  can  be  estimated  from 


(L 


4/AM2 
7 xD 


(4.1) 


where  f  is  the  focal  length  of  the  focusing  lens  (35  mm),  A  is  the  laser  wavelength 


84 


E 


E  =  Laser  pulse  energy 

td  =  Gate  delay  (length  of  time  between  the  laser  pulse 
and  the  spectrometer  turning  on) 
tb=  Gate  width  (length  of  time  the  spectrometer 
remains  on) 


Figure  4-2:  Timing  parameters 


Figure  4-3:  Optical  Configuration.  Mi  =25  mm  diameter  1064  nm  Nd:YAG  mirror; 
L)  =  12  mm  x  -12  mm  lens;  L2  =  25  mm  x  50  mm  lens;  M2  =  50  mm  diameter  1064 
nm  Nd:YAG  mirror;  L3  =  25  mm  x  35  mm  lens;  FO  =  Fiber  Optic 
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(10G4  nni),  M2  is  the  beam  propagation  ratio  which  is  typically  2  -  10  for  Nd:YAG 
lasers  (we  therefore  choose  a  value  of  6),  and  D  is  the  diameter  of  the  illuminated 
aperture  of  the  focusing  lens  («  25  mm)  [24].  The  beam  waist  width  for  the  system 
is  approximately  0.07  mm.  The  average  irradiance  (I/)  at  the  beam  waist  is 


7 tElI)2 

f  4rL/2A2M4 


(4.2) 


where  is  the  laser  pulse  energy  and  rp  (=  7.5  ns)  is  the  pulse  duration  at  the'  full 
peak  width  at  half  of  the  maximum  intensity  (FWHM)  [24].  The  pulse  energy  was 
varied  between  10  -  170  mJ  resulting  in  variation  of  irradiance  at  the  beam  waist  from 
«  1.31  x  1012  to  2.23  x  10u  W/cm2.  The  gate  delay  was  50,  75,  100.  125,  150.  200, 
300,  and  500  ns.  Each  combination  of  energy  and  gate  delay  was  tested,  resulting  in 
13G  different  conditions  for  the  optimization  studies. 

The  signal-to-background  ratio  is 


SBR  =  20/o.9io 


Amplihulv^gk 
,4  [II  pi  i  t  udt’barkgroujid 


(4.3) 


where  the  amplitude  of  the  background  is  defined  as  the  spectral  average  over  a  region 
where  no  peaks  are  expected.  The  background  was  calculated  for  Na  by  using  the 
spectral  region  200  nrn  to  500  nm  and  for  Mn  and  Ca  by  using  the  spectral  region 
430  nm  to  530  nm. 

Calibration  curves  were  made  for  Na,  Ca,  and  Mn  with  ten  spectra  being  taken  at 
each  concentration,  each  composed  of  100  accumulated  shots.  The  experimentation 
conditions  used  for  the  calibration  curves  are  detailed  in  Table  4.1. 


Table  4.1:  Calibration  curve  conditions 


Analyte 

Concentrations  Tested  (ppm) 

E  (mJ) 

U  (ns) 

Na 

0.05,  0.1,  0.5,  1,  5,  10,  50,  100,  500.  1000 

40 

50 

Mn 

1,  5,  10,  50.  100,  500.  1000 

30 

50 

Ca 

1,  5,  10,  50,  100,  500,  1000 

30 

50 

For  all  experiments,  the  gate  width  was  held  constant  at  200  ns.  In  addition,  the 
amplification  of  the  Echelle  spectrometer  was  set  to  the  maximum  value  of  4000. 

All  raw  spectra  were  processed  using  extreme  value  distribution  statistics  detailed 
in  a  paper  by  Michel  and  Chave  [25].  Data  from  9  wavelengths  were  grouped  for 
processing.  Where  shown,  error  bars  represent  the  double  sided  95%  confidence  limits 
[25]. 
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Solutions  were  made  from  NaCl,  MnSOj-HjO.  and  CaCl2-2H20  dissolved  in  de¬ 
ionized  water  for  the  Na,  Mn,  and  Ca  studies,  respectively.  All  concentrations  are 
given  in  parts  per  million  (ppm,  wt./vol.). 

4.4  Results  and  Discussion 

4.4.1  Sodium 

The  interrelationship  of  gate  delay  and  laser  pulse  energy  for  sodium  was  studied  using 
a  concentration  of  100  ppm.  The  intensities  of  the  588.995  nm  and  589.6  mn  Na  peaks 
were  measured,  and  are  shown  in  Figures  4-4  and  4-5.  As  pressure  rises,  an  increase 
in  signal  intensity  is  observed  with  the  maximum  peak  intensity  present  at  2.76  x  10' 
Pa.  An  examination  of  these  figures  shows  that  the  greatest  peak  intensity  exists  at 
the  shortest  gate  delay.  As  t^  increases,  peak  intensity  decreases,  independent  of  both 
the  laser  pulse  energy  and  ambient  pressure.  The  data  variability  present  in  the  plots 
is  indicative  of  significant  plasma  variability.  Overall,  there  appears  to  be  little  effect 
of  energy  on  intensity.  Examination  of  the  SBR  provides  important  information  for 
selecting  optimal  parameters.  Figures  4-6  and  4-7  detail  the  interrelationship  of  t 
energy,  and  SBR.  A  smaller  t,i  tends  to  exhibit  a  higher  SBR  due  to  a  stronger  signal, 
as  seen  in  Figures  4-4  and  4-5.  Furthermore,  Figure  4-4  suggests  that  a  lower  energy 
pulse  consistently  provides  a  higher  SBR.  As  pressure  increases,  the  SBR  again  tends 
to  increase.  The  data  suggest  that  the  highest  SBR  exists  when  a  low  energy  pulse 
(20  -  60  mJ)  and  a  relatively  small  t^  are  used  (50  -  150  ns).  As  first  reported  in 
Michel  et  al.  [17],  this  suggests  that  an  optimal  range  of  laser  energies  exists  that 
tend  to  be  relatively  low. 

For  the  best  SBR,  a  pulse  energy  of  40  mJ  and  a  gate  delay  of  50  ns  were  identified 
for  detection  of  Na  over  a  range  of  pressures  and  spectra  of  this  condition  are  plotted  in 
Figure  4-8.  Calibration  curves  for  sodium  were  constructed  to  determine  the  limit  of 
detection  (Figure  4-9).  These  suggest  that  sodium  can  be  detected  at  a  concentration 
of  approximately  50  ppm  using  the  present  apparatus.  Spectra  of  the  Na  calibration 
data  are  shown  in  Figure  4-10  that  illustrate  that  50  ppm  is  the  lowest  concentration 
where  Na  is  detectable. 
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Figure  4-4:  Interrelationship  of  gate  delay,  laser  pulse  energy,  and  peak  intensity  for 
Na  (588.995  nra)  over  a  range  of  pressures  (a)  1  x  105  Pa  (b)  6.89  x  106  Pa  (c)  1.38 
x  107  Pa  (d)  2.07  x  107  Pa  (e)  2.76  x  107  Pa 
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Figure  4-5:  Interrelationship  of  gate  delay,  laser  pulse  energy,  and  peak  intensity  for 
Na  (589.6  nm)  over  a  range  of  pressures  (a)  1  x  10°  Pa  (b)  6.89  x  106  Pa  (c)  1.38  x 
107  Pa  (d)  2.07  x  107  Pa  (e)  2.76  x  107  Pa 
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Figure  4-6:  Interrelationship  of  pressure,  gate  delay,  energy,  and  signal-to-background 
for  Na  (588.995  nm)  (a)  1  x  105  Pa  (b)  6.89  x  106  Pa  (c)  1.38  x  107  Pa  (d)  2.07  x 
107  Pa  (e)  2.76  x  107  Pa 
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Figure  4-7:  Interrelationship  of  pressure,  gate  delay,  energy,  and  signal-to-background 
for  Na  (589.6  nm)  (a)  1  x  105  Pa  (b)  6.89  x  106  Pa  (c)  1.38  x  107  Pa  (d)  2.07  x  107 
Pa  (e)  2.76  x  107  Pa 
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Figure  4-8:  Spectra  of  Na  (588.995  nm  and  589.6  nm)  taken  with  a  pulse  energy  of 
40  rnJ  and  a  gate  delay  of  50  ns.  From  bottom  to  top,  the  spectra  were  taken  at  1 
x  105  Pa,  6.89  x  10fi  Pa,  1.38  x  107  Pa,  2.07  x  107,  and  2.76  x  107  Pa.  For  clarity, 
the  spectra  have  been  offset  from  each  other  by  8000  a.u. 


92 


10” 


10 


1  io2 


C/3 

I  10° 

c 


10 


-2 


10 


10 


-1 


10" 


B 

6 


# 


10  10 
Concentration  (ppm) 


10 


10 


(a) 


10  10 
Concentration  (ppm) 

(b) 


Figure  4-9:  Sodium  calibration  curves  of  the  (a)  588.995  nm  peak  and  the  (b)  589.6 
nm  peak.  0=  1  x  105  Pa,  □  =  1.38  x  107  Pa,  A  =  2.76  x  107  Pa 
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Figure  4-10:  Spectra  of  sodium  (588.995  nm  and  589.6  nm)  at  2.76  x  1()7  Pa  made 
over  a  range  of  NaCl  concentrations  The  concentrations  from  bottom  to  top  are  10 
ppm,  50  ppm.  100  ppm,  500  ppm,  and  1000  ppm.  For  clarity,  the  spectra  have  been 
offset  from  each  other  by  1000  a.u. 
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4.4.2  Manganese 

The  effect  of  laser  pulse  energy  and  gate  delay  on  spectra  for  manganese  were  studied 
using  a  concentration  of  1000  ppm  Mn  over  a  range  of  pressures.  Although  a  Mn 
triplet  exists  at  403  mn,  peak  broadening  in  liquids  causes  it  to  be  unresolvable,  and 
therefore  we  report  a  single  403  mil  peak.  The  interrelationship  of  pressure,  gate 
delay,  energy,  and  intensity  are  shown  in  Figure  4-11.  As  pressure  increases,  peak 
intensity  also  rises.  A  similar  finding  was  reported  by  Michel  et  al.  under  a  single 
tested  condition  (single  E  and  t^)  [17].  Figure  4-11  shows  that  irrespective  of  gate 
delay  and  energy  the  peak  intensity  of  Mn  increases  with  pressure.  Mn  also  exhibits 
a  higher  peak  intensity  at  a  lower  laser  pulse  energy  at  a  shorter  trf.  When  the 
corresponding  SDR  plots  are  examined  (Figure  4-12),  SDR  is  shown  to  be  smallest 
at  the  lowest  pressure  (1  x  10°  Pa).  Again,  the  need  for  a  short  t d  and  a  low  E  is 
evident. 

From  the  optimization  studies,  an  energy  of  30  mJ  with  a  gate  delay  of  50  ns  was 
selected  as  a  condition  that  would  provide  good  detection  of  Mn  over  a  broad  range 
of  pressures.  The  selected  condition  is  plotted  at  all  five  pressure  conditions  in  Figure 
4-13.  To  determine  the  limit  of  detection  of  Mn,  a  calibration  curve  was  constructed 
using  30  mJ  and  t of  50  ns  (Figure  4-14).  Figure  4-15  shows  spectra  made  at  these 
conditions  at  2.76  x  107  Pa  over  a  range  of  concentrations.  The  limit  of  detection 
was  found  to  be  500  ppm  which  is  higher  than  the  concentration  found  in  vent  fluids. 
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Figure  4-11:  Interrelationship  of  pressure,  gate  delay,  energy,  and  intensity  for  Mil 
(403.076  nm)  (a)  1  x  105  Pa  (b)  6.89  x  106  Pa  (c)  1.38  x  107  Pa  (d)  2.07  x  107  Pa 
(e)  2.76  x  107  Pa 
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Figure  4-12:  Interrelationship  of  pressure,  t^,  E,  and  SBR  for  Mn  (403.076  inn)  (a)  1 
x  105  Pa  (b)  6.89  x  106  Pa  (c)  1.38  x  107  Pa  (d)  2.07  x  107  Pa  (e)  2.76  x  107  Pa 
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Figure  4-13:  Manganese  spectra  using  a  30  mJ  energy  pulse  and  a  gate  delay  of  50 
ns.  (a)  Full  spectral  region,  (b)  Manganese  (403  nm  peak)  In  (a)  and  (b)  the  spectra 
from  bottom  to  top  are  at  1  x  105  Pa,  6.89  x  106  Pa,  1.38  x  107  Pa,  2.07  x  107  Pa, 
2.76  x  107,  respectively.  For  clarity,  the  spectra  have  been  offset  from  each  other  by 
2000  a.u. 
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Figure  4-14:  Mn  (403  nm)  calibration  curve.  0=  1  x  105  Pa,  □  =  1.38  x  107  Pa, 
A  =  2.76  x  107  Pa. 
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Figure  4-15:  Spectra  of  manganese  at  2.76  x  107  Pa  made  at  a  range  of  Mn  concen¬ 
trations  (a)  Full  spectral  region  (b)  Manganese.  In  (a)  and  (b)  the  concentrations  of 
the  spectra  from  bottom  to  top  are  100  ppm,  500  ppm,  and  1000  ppm.  For  clarity, 
the  spectra  have  been  offset  from  each  other  by  100  a.u. 
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4.4.3  Calcium 


Three  calcium  peaks,  393  nm  (ionic),  396  nm  (ionic)  and  422  nm  (atomic),  were 
studied.  The  interrelationship  of  the  measurement  parameters  of  Ca  were  determined 
and  are  shown  in  Figures  4-16  to  4-21.  The  importance  of  a  short  t^  is  evident  when 
SBR  is  examined  (Figures  4-19  to  4-21).  Although  three  calcium  peaks  are  detectable, 
the  422  nm  peak  is  the  strongest,  and  therefore  the  selection  of  an  optimal  condition 
was  based  on  this  peak.  A  laser  energy  pulse  of  30  mJ  with  a  gate  delay  of  50  ns 
was  selected  as  the  optimal  condition  for  detection  of  Ca,  and  spectra  illustrating 
this  condition  are  shown  in  Figure  4-22.  Calcium  calibration  curves  were  constructed 
using  these  conditions  and  are  shown  in  Figure  4-23  for  both  the  393  nm  and  the  422 
nm  peaks.  These  suggest  the  limit  of  detection  for  Ca  is  50  ppm  using  the  present 
apparatus.  Spectra  for  selected  concentrations  are  illustrated  in  Figure  4-24. 


Figure  4-16:  Interrelationship  of  pressure,  t^,  E,  and  intensity  for  Ca  (393  nm)  (a)  1 
x  10s  Pa  (b)  1.38  x  107  Pa  (c)  2.76  x  107  Pa 
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Figure  4-17:  Interrelationship  of  pressure,  t^,  energy,  and  intensity  for  Ca  (390  nm) 
(a)  1  x  105  Pa  (b)  1.38  x  107  Pa  (c)  2.76  x  107  Pa 
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Figure  4-18:  Interrelationship  of  pressure,  t^,  E,  and  intensity  for  Ca  (422  nm)  (a)  1 
x  105  Pa  (b)  1.38  x  107  Pa  (c)  2.76  x  107  Pa 
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Figure  4-19:  Interrelationship  of  pressure,  t,*,  E,  and  signal-to-background  for  Ca  (393 
nm)  (a)  1  x  105  Pa  (b)  1.38  x  107  Pa  (c)  2.76  x  107  Pa 
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Figure  4-20:  Interrelationship  of  pressure,  t^,  E,  and  signal-to-background  for  Ca  (396 
nm)  (a)  1  x  105  Pa  (b)  1.38  x  107  Pa  (c)  2.76  x  107  Pa 
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Figure  4-21:  Interrelationship  of  pressure,  tj,  E,  and  signal-to-background  for  Ca  (422 
nm)  (a)  1  x  105  Pa  (b)  1.38  x  107  Pa  (c)  2.76  x  107  Pa 
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Figure  4-22:  Calcium  spectra  using  30  mJ  and  a  50  ns  gate  delay,  (a)  Full  spectrum 
(b)  Calcium  peaks  (393  nm,  396  nm,  and  422  nm).  Both  (a)  and  (b)  spectra  from 
bottom  to  top:  1  x  105  Pa,  1.38  x  107  Pa,  and  2.76  x  107  Pa.  For  clarity,  the  spectra 
have  been  offset  from  each  other  by  15,000  a.u. 


107 


Intensity  (a.u.)  Intensity  (a.u.) 


(a)  Calcium  (3911  mu)  calibration  curve 


10  10' 
Concentration  (ppm) 

(b)  Calcium  (422  nm)  calibration  curve 


Figure  4-23:  Calcium  calibration  curves  (a)  393  nm  (b)  422  nm,  0=  1  x  105  Pa,  □ 
=  1.38  x  107  Pa,  A  =  2.76  x  107  Pa. 
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Figure  4-24:  Calcium  spectra  at  2.76  x  107  Pa,  Calcium  peaks  are  present  at  393 
nm,  396  nm,  and  422  nm.  (a)  Spectra  from  bottom  to  top  are  10  ppm,  50  ppm,  100 
ppm,  500  ppm,  1000  ppm  respectively.  For  clarity,  the  spectra  have  been  offset  from 
each  other  by  2000  a.u.  (b)  Spectra  from  bottom  to  top  are  10  ppm,  50  ppm,  and 
100  ppm  respectively.  For  clarity,  the  spectra  have  been  offset  from  each  other  by  50 
a.u. 
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4.5  Conclusions 


Sodium,  manganese,  and  calcium  are  all  detectable  in  high  pressure,  bulk  aqueous 
solutions  using  single  pulse  LIBS  with  little  effect  of  pressure  on  the  spectra.  This 
comprehensive  study  of  the  interrelationship  of  gate  delay  and  energy  for  selecting 
the  optimal  condition  for  detection  of  these  analytes  has  shown  that,  irrespective  of 
the  laser  pulse  energy  selected,  the  gate  delay  should  be  very  short  (less  than  200  ns) 
for  the  detection  of  analytes  in  bulk  aqueous  solutions.  The  need  for  a  short  gate 
delay  is  independent  of  pressure.  This  study  has  also  shown  that  a  low  energy  pulse 
(less  than  ~  60  mJ)  is  optimal. 

Calibration  curves  were  made  to  determine  limits  of  detection  using  the  current 
system  set-up  and  further  work  is  needed  to  look  at  reproducibility  of  the  actual 
curves.  Calibration  curves  show  that  Na,  Mil,  and  Ca  can  be  detected  at  50  ppm,  500 
ppm,  and  50  ppm,  respectively.  The  calibration  curves  also  demonstrate  a  minimal 
effect  of  pressure  on  spectra.  However,  the  limits  of  detection  were  higher  than 
expected.  With  the  current  LIBS  set-up,  the  detection  limits  of  Na  and  Ca  are  below 
the  levels  found  in  vent  fluids.  However.  Mn  would  not  be  detectable.  This  can  be 
attributed  to  the  low  light  throughput  of  the  // 10  Echelle  system.  To  significant  ly 
improve  the  light  throughput,  it  would  be  advisable  to  use  a  spectrometer  with  a 
smaller  /  number.  For  example,  using  a  spectrometer  with  an  /  number  of  2  could 
improve  the  throughput  by  a  factor  of  approximately  25  and  therefore  improve  the 
ability  to  detect  Na,  Mn,  and  Ca  in  high  pressure  aqueous  environments.  For  example, 
with  the  use  of  a  PMT  as  the  detector,  Cremers  et  nl.  [14]  showed  substantially 
improved  detection  limits  for  bulk  liquids  for  Na  I  (589.00  mil)  at  a  level  of  0.014 
ppm  and  Ca  II  (393.37  mil)  at  a  level  of  0.8  ppm.  Therefore,  additional  work  is 
necessary  to  optimize  the  light  collection  by  changing  the  system  components. 
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Chapter  5 


Double  pulse  laser-induced 
breakdown  spectroscopy  of  bulk 
aqueous  solutions  at  oceanic 
pressures:  Interrelationship  of  gate 
delay,  pulse  energies,  interpulse 
delay,  and  pressure 

5.1  Abstract 

Laser-induced  breakdown  spectroscopy  (LIBS)  has  been  identified  as  an  analytical 
chemistry  technique  suitable  for  field  use.  Mid-ocean  ridge  hydrothermal  vents  would 
greatly  benefit  from  the  development  of  an  in  situ  LIBS  sensor.  In  this  paper,  double 
pulse  laser-induced  breakdown  spectroscopy  is  used  to  detect  five  analytes  (sodium, 
manganese,  calcium,  magnesium,  and  potassium)  that  are  of  key  importance  in  un¬ 
derstanding  the  chemistry  of  hydrothermal  vent  fluids,  and  of  mixtures  of  vent  fluids 
and  seawater.  The  high  pressure  aqueous  environment  of  the  deep  ocean  is  simulated 
in  the  laboratory  and  the  key  double  pulse  experimental  parameters  (laser  pulse  en¬ 
ergies,  gate  delay,  and  interpulse  delay)  are  studied  at  pressures  up  to  2.76  x  10'  Pa. 
Each  element  is  found  to  have  a  unique  optimal  set  of  parameters  for  detection,  and 
the  elements  are  not  detectable  outside  of  the  set.  For  all  pressures  and  energies,  a 
short  (<100  ns)  gate  delay  is  necessary.  As  pressure  increases,  a  shorter  interpulse 
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delay  is  needed  and  the  conditions  effectively  become  single  pulse.  Calibration  curves 
reveal  the  limits  of  detection  of  the  elements  (5000  ppm  Mg,  500  ppm  K,  500  ppm 
Ca,  1000  ppm  Mn,  and  50  ppm  Na).  When  compared  to  our  previous  single  pulse 
work  for  Ca,  Mn,  and  Na,  double  pulse  LIBS  for  high  pressure  aqueous  solutions  did 
not  improve  the  limits  of  detection. 

5.2  Introduction 

Laser-induced  breakdown  spectroscopy  (LIBS)  has  recently  been  identified  as  an  an¬ 
alytical  chemistry  technique  suitable  for  Held  deployment  for  the  analysis  of  environ¬ 
mental  and  geochemical  samples  [1].  Several  groups  are  evaluating  the  use  of  LIBS  for 
space  exploration  [2-7].  Another  environmental  area  in  which  there  is  a  critical  need 
for  new  sensors  is  the  ocean.  Sensor  development  for  use  with  underwater  vehicles  is 
ongoing  and  a  more  requirement  is  for  sensors  to  be  developed  for  use  on  permanent 
ocean  observatories. 

One  environment  identified  as  potent  ially  benefitting  from  applicat  ion  of  an  oceanic 
LIBS  sensor  is  the  deep-sea  hydrothermal  vent  environment.  Hydrothermal  vents  oc¬ 
cur  at  mid-ocean  ridges  where  seawater  circulates  through  the  permeable  ocean  crust . 
allowing  the  fluid  to  interact  with  the  surrounding  rock  and  resulting  in  major  fluid 
chemical  changes.  At  vent  orifices,  exit  temperatures  reach  200  -  405"C  at  ambient 
pressures  of  8.1  x  10u  Pa  to  3.6  x  10'  Pa  corresponding  to  ocean  depths  of  800  m  to 
3600  m  [8]. 

At  hydrothermal  vents,  as  the  hot  fluids  mix  with  seawater,  rapid  chemical  changes 
occur  and  some  elements  precipitate  out  (e.g.,  sulfate  minerals  [9]).  Obtaining  in  situ 
chemical  measurements  of  the  fluid  is  difficult  due  to  the  corrosive  nature  and  high 
temperature  of  the  fluid.  Collection  of  the  fluid  for  analysis  shipboard  or  in  a  lab¬ 
oratory  introduces  chemical  changes  as  the  temperature  and  pressure  of  the  fluid  is 
changed  during  sample  recovery.  Five  critical  elements  at  hydrothermal  vents  are 
sodium,  calcium,  manganese,  magnesium,  and  potassium.  Sodium  is  the  dominant 
cation  in  vent  fluids,  and  provides  insight  into  phase  separation  processes  [10].  Cal¬ 
cium  is  the  second  most  dominant  cation  in  vent  fluids,  and  is  usually  found  at  a 
greater  concentration  in  vent  fluids  than  seawater  [11].  Ca  is  released  into  vent  fluids 
when  Na  is  taken  up  in  albitization  reactions  with  the  host  rock  [11].  Manganese 
exists  as  a  trace  metal  in  seawater  but  has  a  higher  concentration  in  vent  fluids  due 
to  leaching  from  the  host  rock  [10].  Magnesium  is  practically  nonexistent  in  hy¬ 
drothermal  vent  fluids;  however,  if  any  is  detected  in  vent  fluids,  contamination  by 
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entrainment  of  ambient  seawater  is  indicated  [10].  Potassium  is  typically  highly  en¬ 
riched  in  vent  fluids  due  to  leaching  from  basalts  [10].  In  vent  fluids,  concentrations 
range  from  approximately  250  -  23,163  ppm  for  Na,  0.6  -  399  ppm  for  Mn,  -54  -  4477 
ppm  for  Ca,  -47  -  3166  ppm  for  K,  and  0  ppm  for  Mg  [8].  In  seawater,  concentrations 
are  approximately  10,933  ppm  Na,  <0.001  Mn,  419  ppm  Ca,  405  ppm  K,  and  1300 
ppm  Mg  [8]. 

Development  of  an  oceanic  LIBS  sensor  necessitates  laboratory  investigations  into 
the  system  parameters  for  the  detection  of  analytes  under  high  pressure  bulk  aqueous 
conditions.  Although  LIBS  analysis  of  liquids  is  more  difficult  than  the  analysis 
of  solid  or  air  samples,  a  few  studies  have  focused  on  dissolved  analytes  in  bulk 
solutions  [12-18].  Some  issues  for  LIBS  in  aqueous  solutions  include  a  reduction  in 
plasma  light  intensity  and  emission  lifetime  due  to  quenching,  the  Stark  effect  may 
cause  spectral  lines  to  be  broadened,  and  moving  breakdown  can  change  the  distance 
between  the  plasma  and  the  collection  fiber  optic  [12,  13,  19-21].  The  breakdown 
threshold  in  water  is  also  significantly  greater  than  for  solids  [22].  Laser-induced 
plasmas  are  weak  in  water,  as  they  are  cooled  by  Bremsstrahlung  and  shockwave 
emissions,  and  by  thermal  conduction.  Rapid  cooling  of  the  plasma  also  increases 
electron-ion  recombination  and  plasma  emissions  last  only  a  few  hundred  nanoseconds 

[23]- 

In  1984,  Cremers  et  al.  showed  that  several  elements,  including  Na,  K,  Mg,  and 
Ca,  could  be  identified  in  bulk  aqueous  solution  using  double  pulse  LIBS.  This  work, 
carried  out  at  atmospheric  pressure,  showed  that  double  pulse  LIBS  improved  the 
detection  limit  for  metals  and  ions  in  bulk  aqueous  solution  [12].  In  double  pulse 
LIBS  for  dissolved  analytes,  the  first  laser  pulse  produces  a  plasma  that,  creates  a 
laser-induced  cavitation  bubble.  The  second  laser  pulse  produces  a  plasma  within 
the  bubble  [23,  24].  This  is  in  contrast  to  single  pulse  LIBS,  where  the  plasma  is 
simply  formed  in  a  liquid  environment.  The  duration  of  the  laser-induced  plasma  in 
the  bubble  is  on  the  order  of  a  few  microseconds  and  the  bubble  lifetime  is  on  the 
order  of  a  few  hundred  microseconds.  Therefore,  it  can  be  assumed  that  the  plasma 
from  the  second  pulse  is  expanding  in  a  quasi-stationary  environment  induced  by 
the  first  laser  pulse  [24,  25].  When  the  bubble  is  first  formed,  its  pressure  is  greater 
than  that  of  the  surrounding  liquid  and  the  bubble  begins  to  expand  which  leads  to 
a  pressure  drop.  At  the  point  of  maximum  expansion,  the  bubble  pressure  is  less 
than  the  pressure  of  the  surrounding  fluid  and  the  bubble  begins  to  collapse.  During 
this  collapsing  phase,  the  temperature  and  pressure  in  the  bubble  again  increases 
and  if  there  is  enough  energy  stored  within  the  bubble,  it  can  re-expand.  Many 
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such  oscillations  of  expansion  and  compression  are  possible  [23].  By  adjusting  the 
interpulse  delay  between  the  laser  pulses,  it  is  possible  to  select  the  phase  that  the 
bubble  is  in  such  that  bubble  pressure  is  at  a  nadir.  If  the  interpulse  delay  is  too 
long  or  too  short,  the  plasma  will  expand  in  a  high  pressure  environment  so  that  the 
spectral  lines  will  be  broadened  due  to  collisions  [23].  Selection  of  the  interpulse  delay 
time  therefore  is  critical  for  bulk  liquids. 

Since  Cremer’s  pioneering  work,  several  groups  have  used  dual  pulse  LIBS  for  the 
detection  of  analytes  in  bulk  liquids  [12,  14,  18,  25-27]  and  more  recently  for  bulk 
solutions  at  high  pressures  [15,  16].  In  Michel  et  al.  [15],  it  was  reported  that  analyte 
detection  in  high  pressure  bulk  solutions  was  highly  dependent  on  the  interpulse 
delay.  When  a  shorter  interpulse  delay  time  was  used  (-C  1  //s) ,  the  signal  intensity 
was  enhanced  compared  to  when  longer  delay  times  were  employed.  However,  it  was 
noted  that  this  time  may  not  be  long  enough  for  cavitation  bubbles  to  fully  form  and 
expand,  as  occurs  at  low  pressure.  Michel  et  al.  also  found  that  the  optimal  energy 
levels  needed  for  emission  seemed  to  vary  by  analyte  [15].  Lawrence-Snyder  et  al.  [16] 
report  that  increasing  solution  pressure  reduces  double  pulse  emission  enhancement 
so  that  little  improvement  was  noted  over  single  pulse  above  1  x  10'  Pa. 

In  Chapter  4,  optimization  of  single  pulse  conditions  for  three  analytes  (Na.  Mu 
and  Ca)  was  reported  and  the  limits  of  detection  were  not  at  the  levels  that  we  had 
hoped  to  achieve.  Through  optimization  of  the  double  pulse  set-up,  the  goal  is  to 
improve  the  limits  of  detection  and,  in  addition,  establish  limits  of  detection  for  other 
elements.  Although  Lawrence-Snyder  et  al.  previously  reported  on  the  use  of  double 
pulses  for  the  detection  of  analytes  at  high  pressures  and  found  no  enhancements 
above  1  x  10'  Pa  [16],  several  differences  exist  between  these  studies.  Lawrence- 
Snyder  et  al.  report  their  findings  using  an  orthogonal  beam  geometry  which  would 
be  impractical  for  the  development  of  an  oceanic  sensor  in  many  applications,  the 
work  does  not  extend  to  the  ambient  pressures  for  most  hydrothermal  vents,  and 
optimization  of  laser  pulse  energies  was  not  carried  out.  Through  the  completion  of  a 
thorough  optimization,  the  use  of  double  pulse  LIBS  at  pressures  up  to  2.76  x  10'  Pa 
is  now  investigated  and  the  limits  of  detection  for  five  key  elements  are  determined. 


5.3  Experimental 

Double  pulse  high  pressure  aqueous  LIBS  experiments  were  completed  using  the  labo¬ 
ratory  set-up  detailed  in  Figure  5-1.  Two  Big  Sky  CFR-200  Nd:YAG  lasers  configured 
with  the  beams  co-linearly  aligned  prior  to  exiting  the  aperture  (Particle  Image  Ve- 
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Figure  5-1:  Laboratory  set-up  for  high  pressure  aqueous  double  pulse  LIBS  experi¬ 
ments 

locimetry  or  PIV  configuration)  were  operated  at  1064  nm  with  a  5  Hz  repetition  rate. 
Each  laser  is  equipped  with  a  motorized  variable  attenuator  serially  controlled  by  a 
computer,  enabling  the  laser  pulse  energies  (Ei  and  E2)  to  be  varied  independently 
from  0  to  200  mJ,  in  increments  of  approximately  1  mJ.  Plasma  emission  is  collected 
with  an  Echelle  spectrometer  (LLA  Echelle  ESA  3000).  The  spectrometer  is  capable 
of  detecting  elements  with  wavelengths  of  200  -  780  nm  at  a  spectral  resolution  of 
10  -  50  pm.  Accurate  timing  of  the  laser  firing,  the  interpulse  delay  time  (AT,  the 
time  between  the  firing  of  laser  pulse  1  and  laser  pulse  2)  and  the  gate  delay  (t^.  the 
time  between  the  firing  of  the  second  laser  pulse  and  the  turn-on  of  the  spectrometer) 
was  controlled  by  a  timing  box  (Berkeley  Nucleonics  Corporation  Model  565).  The 
integration  time  of  the  spectrometer,  gate  width  (t^),  is  also  controllable.  The  double 
pulse  parameters  are  clarified  in  Figure  5-2. 

An  8.89  cm  x  8.89  cm  x  8.89  cm  titanium  sample  chamber  that  holds  27  ml 
of  liquid  connects  to  a  high  pressure  metering  pump  (Eldex  Model  A-30-S)  using 
Swagelok  fittings  to  pressurize  samples  up  to  4.1  x  107  Pa.  The  sample  chamber 
is  equipped  with  a  sapphire  window  (Meller  Optics,  2.54  cm  diameter  x  0.64  cm 
thickness,  AR-coated  at  532  mn/1064  nm,  custom  part)  that  allows  laser  pulses  to 
enter  the  sample  chamber.  A  series  of  AR-coated  optics  are  used  to  focus  the  laser 
beams  into  the  sample  chamber  and  to  focus  the  plasma  light  onto  an  optical  fiber  for 
delivery  to  the  spectrometer.  The  plasma  light  is  collected  collinear  to  the  incoming 
laser  beam  to  simulate  the  design  that  would  be  most  practical  for  an  ocean-going 
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E,  =  Laser  1  pulse  energy 
E2  =  Laser  2  pulse  energy 

AT  =  Interpulse  delay  time  (length  of  time  between  the 
firing  of  laser  pulse  1  and  laser  pulse  2) 
td  =  Gate  delay  (length  of  time  between  laser  pulse  2 
and  the  spectrometer  turning  on) 
tb  =  Gate  width  (length  of  time  the  spectrometer 
remains  on) 


Figure  5-2:  Double  pulse  LIBS  timing  parameters 

LIBS  system.  Data  are  collected  using  ESAWIN  software.  Laser  energy  is  measured 
using  a  laser  energy  sensor  (Coherent  J25LP-MB)  combined  with  an  energy  meter 
(Coherent  FieldMaxII-Top) . 

For  double  pulse  optimization  studies,  peak  intensities  for  five  analytes  (Mg.  K. 
Na,  Ca,  and  Mn)  were  measured  over  a  range  of  system  parameters:  Ej ,  E2,  t,i,  and 
AT.  Ei  and  E2  were  each  tested  at  20  mJ,  60  mJ,  100  111J,  and  140  mJ.  The  average 
irradiance  (If)  at  the  beam  waist  can  be  estimated  from 

7 rELD2 

1  4r/J'2A2A/4  ^ 

where  E 1  is  the  laser  pulse  energy,  D  is  the  is  the  diameter  of  the  illuminated  aperture 
of  the  focusing  lens  («  25  mm),  and  tl  is  the  pulse  duration  at  the  full  peak  width 
at  half  of  the  maximum  intensity  (FWHM;  tl,  =  7.5  11s),  f  is  the  focal  length  of  the 
focusing  lens  (35  mm),  A  is  the  laser  wavelength  (1064  nm),  and  M2  is  the  beam 
propagation  ratio  that  is  typically  2-10  for  Nd:YAG  lasers  (we  estimate  a  value 
of  6)  [28].  The  pulse  energies  of  20  mJ,  60  mJ,  100  mJ,  and  140  mJ  correspond  to 
irradiance  at  the  beam  waist  of  «  2.62  x  1012  W/cm2,  7.87  x  1012  W/cm2,  1.31  x 
1013  W/cm2,  1.84  x  1013  W/cm2,  respectively. 

For  optimization  of  experimental  parameters,  studies  were  conducted  at  three 
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pressures,  1  x  10°  Pa,  1.38  x  107  Pa,  and  2.76  x  107  Pa,  which  correspond  to  ocean 
depths  of  approximately  0  in,  1362  m,  and  2724  m  respectively.  For  Mg  and  K.  the 
gate  delay  was  tested  at  10  ns,  50  ns,  100  ns,  and  500  ns  and  the  interpulse  delay 
was  tested  at  50  ns,  100  ns,  500  ns,  1000  ns,  and  5000  ns.  For  Na,  Ca,  and  Mn, 
the  gate  delay  was  held  constant  at  50  ns  and  the  interpulse  delay  was  tested  at 
50  ns,  500  ns,  and  5000  ns.  Each  combination  of  these  conditions  was  evaluated 
for  a  total  of  320  conditions  for  Mg  and  K  and  48  conditions  for  Mn,  Na,  and  Ca. 
For  all  studies,  the  gate  width  was  held  constant  at  200  ns  and  the  amplification 
of  the  Echelle  spectrometer  was  set  to  the  maximum  value  of  4000.  Five  spectra 
were  taken  at  each  optimization  condition  and  for  calibration  curves  ten  spectra  were 
taken.  All  spectra  were  comprised  of  100  accumulations.  To  account  for  the  high 
resolution  of  the  spectrometer  and  the  peak  broadening  that  occurs  from  liquids,  all 
data  were  grouped  into  sets  of  9  wavelengths.  Each  set  of  datapoints  (9x5  for  the 
optimization  or  9  x  10  for  the  calibration  curves)  was  then  processed  using  extreme 
value  distribution  statistics  described  by  Michel  and  Chave  [29]. 

All  calibration  curves  were  made  at  three  pressures  (1  x  10;>  Pa,  1.38  x  107  Pa,  and 
2.76  x  10"  Pa).  Where  shown,  error  bars  represent  the  double-sided  95%  confidence 
limits  for  the  extreme  value  parameters  defined  in  [29]. 

Solutions  were  made  using  MgCl2-6H20,  NaCl.  KC1,  M11SO4H2O.  and  CaCl2-H2() 
dissolved  in  DI  water.  All  concentrations  are  given  in  parts  per  million  (ppm. 
wt./vol.).  For  the  optimization  studies  the  concentrations  used  were  5000  ppm  Mg. 
1000  ppm  K.  100  ppm  Na,  1000  ppm  Ca,  and  1000  ppm  M11. 


5.4  Results  and  Discussion 

For  double  pulse  LIBS  of  high  pressure  aqueous  solutions,  optimizing  the  key  pa¬ 
rameters  (Ej,  E2,  id,  and  AT)  individually  for  each  element  of  interest  is  essential 
for  identifying  the  conditions  under  which  each  can  be  detected.  The  optimization 
studies  presented  here  show  that  these  conditions  are  pressure  dependent  for  double 
pulse  LIBS.  Outside  the  range  of  these  conditions,  some  of  the  elements  prove  to 
be  undetectable.  Through  the  optimization  of  the  parameters,  a  set  of  conditions 
are  established  that  allow  calibration  curves  to  be  made  for  determining  the  limits 
of  detection  for  high  pressure  bulk  aqueous  solutions  using  the  current  LIBS  system 
setup. 
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5.4.1  Magnesium 


Results  for  the  518.4  nm  Mg  (I)  peak  are  presented  in  Figures  5-3  to  5-5.  From 
these  figures,  the  need  for  a  short  gate  delay  (<  100  ns)  is  evident,  irrespective  of 
the  pressure  or  the  laser  pulse  energies.  At  all  pressures,  when  t^  =  500  ns.  Mg  (I) 
is  not  detectable.  At  1  x  10°  Pa.  the  intensity  was  greatest  with  a  longer  interpulse 
delay  (1000  ns  to  5000  ns)  (Figure  5-3).  It  was  also  best  to  use  a  smaller  E|  than 
E2.  For  example,  two  favorable  conditions  were  E|  =  60  mJ,  E2  =  100  mJ  and  E| 

20  mJ,  E2  =  100  mJ.  At  2.70  x  107  Pa,  the  greatest  intensity  peak  exists  when  AT 
had  the  smallest  value  (50  ns)  (Figure  5-5).  The  Ei  and  E2  that  gave  the  greatest 
intensity  were  both  in  the  range  of  60  mJ  -  140  mJ.  At  1.38  x  10'  Pa.  the  optimal 
AT  was  intermediate  between  that  for  1  x  10°  Pa  and  2.76  x  10'  Pa.  At  2.76  x  10' 
Pa,  a  very  intense  peak  was  obtained  for  the  condition  E|  =  60  mJ,  E2  =  60  mJ,  t,y 
=  50  ns,  and  AT=  50  ns.  Therefore,  a  calibration  curve  was  constructed  using  these 
conditions  as  shown  in  Figure  5-6.  From  this  calibration  curve  and  from  examining 
spectra  of  the  calibration  data  (Figure  5-7),  it  is  clear  that  Mg  (I)  can  be  detected  to 
only  5000  ppm. 
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Figure  5-3:  Mg  (I)  (518.4  nm  peak)  optimization  at  1  x  105  Pa 
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figure  5-4:  Mg  (I)  (518.4  nm  peak)  optimization  at  1.38  x  107  Pa 


i  '  •  >  '  -  i 

2500 

1  1  1 

• 

• 

2000 

•  t 

d 

1 

5  1500 

•£ 

• 

•  • 

g  1000 
jj: 

• 

• 

• 

500 

f 

:  1 

•  1 

•  •  i  • 

0 

•  •II 

80  100 
E,  (mJ) 


60  80  100 
E  (mJ) 


Interpul.se  Delay  (ns) 


Gate  Delay  (ns) 


Figure  5-5:  Mg  (I)  (518.4  nm  peak)  optimization  at  2.76  x  10'  Pa 


123 


10 


10 


3 

3 

.5  io° 

C/3 

fi 

<u 


f 

B 

e 


$ 


$ 

$ 


$ 

§  5 


$ 


e 


10 


-2 


10 


10 


-1 


10 


10  10 
Concentration  (ppm) 


10 


10 


Figure  5-6:  Mg  (I)  calibration  curve  of  the  518.4  nm  peak,  Q  =  1  x  10’  Pa,  □ 
1.38  x  1()7  Pa,  A  =  2.76  x  107  Pa,  (Ei  =  60  mJ,  E2  =  60  mJ,  td  =  50  ns,  and  AT 
50  ns) 


1400 


200 

?05  510  515  520  525  530 

Wavelength  (nm) 


Figure  5-7:  Spectra  of  Mg  (I)  peak  (518.4  nm)  at  2.76  x  107  Pa.  The  concentrations 
from  bottom  to  top  are  1000  ppm  and  5000  ppm.  (E]  =  60  mJ,  E2  =  60  mJ,  t,* 

50  ns,  and  AT=  50  ns).  For  clarity,  the  spectra  have  been  offset  from  each  other  by 
500  a.u. 
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5.4.2  Potassium 


The  results  of  the  optimization  studies  for  the  769.9  nm  K  (I)  peak  are  shown  in 
Figures  5-8  -  5-10.  Potassium  was  detectable  over  an  unusually  wide  range  of  condi¬ 
tions.  Potassium  has  the  lowest  ionization  energy  (4.31  eV)  of  the  elements  that  were 
studied  which  contributes  to  ease  of  detection.  In  particular,  a  wide  range  of  inter¬ 
pulse  delay  times  was  suitable  at  all  pressures.  At  1  x  10°  Pa  and  1.38  x  107  Pa,  the 
use  of  two  high  energy  pulses  resulted  in  the  greatest  peak  intensity.  At  the  highest 
pressure  condition  (2.76  x  10'  Pa),  a  lower  energy  pulse  followed  by  a  higher  energy 
pulse  was  advantageous.  For  the  calibration  curve,  the  condition  Ei  =  100  mJ,  E2  = 
140  mJ,  trf  =  1000  ns,  and  AT=  50  ns  was  selected  (Figure  5-11).  The  calibration 
curve  reveals  that  K  (I)  is  detectable  to  500  ppm.  Spectra  of  the  calibration  data 
confirm  this  and  are  shown  in  Figure  (5-12). 
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Figure  5-8:  K  (I)  (769.9  nm  peak)  optimization  at  1  x  10r>  Pa 


125 


Intensity  (a.u.)  Intensity  <a.u.)  Intensity  (a.u.)  Intensity  (a.u.) 


II !  I 


2000  3000 

Interpulse  Delay  (ns) 


140 

140 

• 

120 

• 

120 

100 

—  100 

•  • 

d 

80 

t 

*  80 

>> 

60 

• 

•  1 

!  “ 

40  ^ 

• 

1 

-  40 

• 

| 

20  § 

0  1 

| 

X 

20 

20  44) 

60 

84)  100  120  140 

° 

EJ  (mj) 

140 

150 

• 

120* 

104) 

•  • 

3  100 

80  •  # 

3 

c  50 


,1  I  I 


i  I 


80  100 
E,  (mj) 


200  300 

Cate  Delay  (ns) 


t 


I 


I 


144) 


500 


Figure  5-9:  K  (I)  (769.9  nm  peak)  optimization  at  1.38  x  107  Pa 
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Figure  5-10:  K  (I)  (769.9  nm  peak)  optimization  at  2.76  x  1()7  Pa 


126 


10 


i<F 


3 

>>  IA0 
.-a  10 

c/3 

C 

4> 


10 


-2 


10 


10 


-1 


10" 


B 

e 


B 

A 


10 


e 

B 


10 


Concentration  (ppm) 

(a) 


4 


10 


10 


Figure  5-11:  K  (I)  calibration  curve  of  the  769.9  nm  peak,  Q  =  1  x  105  Pa,  □  = 
1.38  x  107  Pa,  A  =  2.76  x  107  Pa,  (Ej  =  100  mJ,  E2  =  140  mJ,  t^  =  1000  ns,  and 
AT=  50  ns) 
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Figure  5-12:  Spectra  of  769.9  nm  K  (I)  peak  at  2.76  x  107  Pa.  Concentrations  of 
spectra  from  bottom  to  top  are  100  ppm,  500  ppm,  1000  ppm.  (Ei  =  100  mJ.  E2  = 
140  mJ,  t d  =  1000  ns,  and  AT=  50  ns).  For  clarity,  the  spectra  have  been  offset  from 
each  other  by  200  a.u. 
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5.4.3  Calcium 


The  results  for  the  422  nm  Ca  (I)  peak  are  shown  in  Figures  5-13  to  5-15.  At  1  x 
10;>  Pa,  the  greatest  intensity  peak  was  detected  when  a  long  interpulse  delay  time 
was  used  (500  to  5000  ns)  with  a  low  energy  pulse  (typically  20  mJ)  followed  by  a 
higher  energy  pulse  (60  mJ  -  140  mJ).  At  both  1.38  x  107  Pa  and  2.76  x  107  Pa, 
the  greatest  intensity  peak  was  detected  when  two  high  energy  pulses  (typically  100 
mJ  -  140  mJ)  were  separated  by  50  ns.  As  a  result,  the  condition  selected  for  the 
calibration  curve  was  E]  =  100  mJ,  E2  =  100  mJ,  t^  =  50  ns,  and  AT=  50  ns.  The 
calcium  calibration  curves  for  both  the  393  nm  Ca  (II)  and  422  nm  Ca  (I)  peaks  are 
shown  in  Figure  5-16.  The  limit  of  detection  for  calcium  using  the  conditions  selected 
is  500  ppm.  Spectra  at  these  conditions  are  shown  in  Figure  5-17. 
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Figure  5-13:  Ca  (I)  (422  nm  peak)  optimization  at  1  x  105  Pa 
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Figure  5-14:  Ca  (I)  (422  nm  peak)  optimization  at  1.38  x  107  Pa 
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Figure  5-15:  Ca  (I)  (422  nm  peak)  optimization  at  2.76  x  107  Pa 
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Figure  5-16:  Calcium  calibration  curves,  O  =  1  x  105  Pa,  □  =  1.38  x  1()7  Pa,  A 
2.76  x  10'  Pa,  (Ei  =  100  mJ,  E2  =  100  mJ,  =  50  ns,  and  AT=  50  ns) 
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Figure  5-17:  Spectra  of  calcium  (393  nm  Ca  (II),  396  nm  Ca  (II),  and  422  nm  Ca  (I) 
peaks)  at  2.76  x  107  Pa.  Concentrations  of  spectra  from  bottom  to  top  are  100  ppm, 
500  ppm,  and  1000  ppm.  (E!  =  100  mJ,  E2  =  100  mJ,  t^  =  50  ns,  and  AT=  50  ns). 
For  clarity,  the  spectra  have  been  offset  from  each  other  by  15,000  a.u. 
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5.4.4  Manganese 


The  optimization  study  results  for  the  403  nm  Mn  (I)  peak  are  shown  in  Figures  5-18 
to  5-20.  Although  a  Mn  (I)  triplet  exists  at  403  nm.  peak  broadening  in  liquids  causes 
it  to  be  unresolvable,  and  hence  a  single  403  nm  peak  is  used  for  this  study.  At  1 
x  10'’  Pa,  a  low  first  energy  pulse  (20  mJ)  followed  by  a  higher  second  energy  pulse 
(GO  mJ  -  140  mJ)  with  a  long  interpulse  delay  time  (5000  ns)  gave  the  highest  peak 
intensity  (Figure  5-18).  When  the  pressure  was  increased,  the  need  for  a  significantly 
shorter  interpulse  delay  time  was  evident  (Figures  5-18  -  5-19).  At  the  highest  pressure 
condition,  two  high  energy  pulses  gave  the  most  intense  peak  (Figure  5-20).  Using  a 
low  second  energy  pulse  (20  mJ)  was  not  beneficial  at  this  pressure  condition.  At  flu* 
intermediate  pressure,  1.38  x  10'  Pa,  several  parameter  combinations  can  be  used. 
Either  a  low  energy  first  pulse  (20  mJ)  followed  by  a  higher  pulse  (GO  mJ  -  140  mJ)  or 
two  higher  energy  pulses  (GO  mJ  -  140  mJ)  were  suitable  characteristics  of  both  the 
low  pressure  and  high  pressure  conditions  are  workable  and  hence  this  is  possibly  a 
transition  pressure.  For  Mn  (I),  the  conditions  that  were  selected  for  the  calibration 
curve  were  Ei  =  100  mJ,  E2  =  GO  mJ,  t.^  =  50  ns,  AT=  50  ns.  with  results  shown 
in  Figure  5-21.  Spectra  at  these  conditions  are  shown  in  Figure  5-22.  From  both  the 
calibration  curve  (Figure  5-21)  and  the  subsequent  spectra  (Figure  5-22),  the  limit  of 
detection  is  1000  ppm. 
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Figure  5-18:  Mn  (I)  (403  nm  peak)  optimization  at  1  x  10r>  Pa 
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Figure  5-19:  Mil  (I)  (403  nm  peak)  optimization  at  1.38  x  107  Pa 
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Figure  5-20:  Mil  (I)  (403  nm  peak)  optimization  at  2.7G  x  107  Pa 
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Figure  5-21:  Mn  (I)  (403  nm  peak)  calibration  curve.  Q  =  1  x  105  Pa,  □  =  1.38  x 
107  Pa,  A  =  2.76  x  107  Pa,  (Et  =  100  mJ,  E2  =  60  mJ,  td  =  50  ns,  and  AT=  50  ns) 
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Figure  5-22:  Spectra  of  Mn  (I)  peak  (403  run)  at  2.76  x  107  Pa.  Concentrations  from 
bottom  to  top  are  100  ppm,  500  ppm,  and  1000  ppm.  (E!  =  100  mJ,  E2  =  60  mJ,  t^ 
=  50  ns,  and  AT=  50  ns).  For  clarity,  the  spectra  have  been  offset  from  each  other 
by  1000  a.u. 
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5.4.5  Sodium 


The  588.995  inn  Na  (I)  peak  from  the  sodium  doublet  was  used  for  optimization 
studies  (Figures  5-23  -  5-25).  The  highest  intensity  at  1  x  10'’  Pa  was  recorded  when 
a  low  energy  pulse  was  followed  by  a  high  energy  pulse  and  a  long  interpulse  delay  time 
was  used  (Figure  5-23).  For  example,  the  greatest  intensity  was  recorded  for  E]  =  20 
mJ,  E‘2  =  140  mJ,  and  AT  =  5000  ns.  At  1.38  x  107  Pa,  the  greatest  intensity  peaks 
were  recorded  when  two  high  energy  pulses  (60  mJ  to  140  mJ)  were  fired  in  rapid 
succession  in  either  order  and  separated  by  50  ns  (Figure  5-24).  Since  the  interpulse 
delay  time  is  very  small,  these  conditions  are  close  to  that  for  single  pulse  operation 
with  a  very  high  energy  pulse.  At  2.76  x  107  Pa,  the  greatest  intensity  peak  again 
exists  when  two  high  energy  pulses  (60  mJ  and  140  mJ)  are  rapidly  fired  in  either 
order  separated  by  50  ns  (Figure  5-25).  The  lowest  intensity  peaks  were  recorded  at 
all  pressures  when  the  second  energy  pulse  was  20  mJ,  suggesting  that  the  second 
pulse  must  be  of  sufficient  energy  or  irradiance  to  excite  or  re-excite  plasma  emission. 
The  sodium  calibration  curve  was  therefore  made  at  Ei  =  60  mJ,  E2  =  140  mJ,  AT 
50  ns,  and  t d  —  50  ns  and  is  shown  in  Figure  5-26.  Spectra  made  at  the  high  pressure 
(2.76  x  107  Pa)  condition  which  clearly  indicate  the  limit  of  detection  of  50  ppm  are 
shown  in  Figure  5-27. 
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Figure  5-23:  Na  (1)  (588.995  nm)  optimization  at  1  x  105  Pa 
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Figure  5-24:  Na  (I)  optimization  at  1.38  x  107  Pa 
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Figure  5-25:  Na  (I)  (588.995  nm)  optimization  at  2.76  x  107  Pa 
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Figure  5-2G:  Na  (I)  calibration  curves  of  the  (a)  588.995  nm  peak  and  the  (b)  589. G 
nm  peak.  Q  =  1  x  105  Pa,  □  =  1.38  x  107  Pa,  A  =  2.7G  x  107  Pa. 
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Figure  5-27:  Spectra  of  the  Na  (I)  doublet  peaks  (588.995  nm  and  589.6  nm)  at  2.76 
x  107  Pa.  Concentrations  from  bottom  to  top  are  10  ppm,  50  ppm,  100  ppm,  500 
ppm,  and  1000  ppm.  (Ej  =  60  mJ,  E2  =  140  mJ.  AT  =  50  ns,  and  t^  =  50  ns).  For 
clarity,  the  spectra  have  been  offset  from  each  other  by  2000  a.u.,  except  for  the  1000 
ppm  spectrum  which  has  been  offset  from  the  500  ppm  spectrum  by  8000  a.u. 
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The  limits  of  detection  for  sodium  as  well  as  the  other  analytes  reported  here  are 
much  higher  than  expected  based  on  earlier  work  [12.  14].  In  an  effort  to  determine 
if  the  low  throughput  (// 1 0)  Echelle  spectrometer  is  the  limiting  factor  an  additional 
calibration  curve  was  completed  at  1  x  10°  Pa  at  an  optimal  atmospheric  pressure 
condition.  Lawrence-Snyder  et  al.  report  with  the  use  of  a  Chromex  Czerny- Turner 
spectrometer  coupled  to  an  ICCD  camera  that  for  Na  at  3.4  x  10(’  Pa  using  E] 

7  m J .  E2  =  48  mJ,  and  t,;  =  1  ps,  the  optimal  interpulse  delay  time  is  between 
approximately  20  and  50  ps  [16].  In  our  work,  first,  we  selected  a  low  energy  pulse 
(20  mJ)  followed  by  a  high  energy  pulse  (140  mJ)  and  the  optimal  interpulse  delay 
time  between  them  was  determined.  This  interpulse  delay  was  then  used  to  create  a 
calibration  curve.  In  the  optimization  study  for  high  pressure,  the  interpulse  timing 
was  only  carried  out  to  5  ps.  The  first  peak  intensity  was  measured  at  1  ps,  then 
at  5  /is,  then  at  5  /is  increments  to  a  maximum  of  170  ps.  Figure  5-28  details  the 
effect  of  interpulse  timing  on  the  intensity  on  the  588.995  nm  Na  (I)  peak.  The 
intensity  is  fairly  uniform  from  10  ps  until  140  ps.  After  140  ps,  the  intensity  drops 
off.  From  this  plot,  it  appears  valid  to  select  an  interpulse  delay  time  between  1-140 
ps;  therefore,  we  select  70  ps  to  construct  a  calibration  curve.  Comparing  this  time 
scale  to  that  shown  in  Lawrence-Snyder  et  al.  [16],  this  seems  consistent  due  to  our 
use  of  a  significantly  higher  second  energy  pulse.  Using  the  70  /is  interpulse  delay 
time,  a  calibration  curve  was  constructed  and  is  shown  in  Figure  5-29.  Although  we 
have  now  used  the  optimal  condition  for  1  x  10°  Pa  with  a  long  interpulse  delay 
time,  the  detection  limit  is  again  only  50  ppm.  This  can  also  be  verified  by  looking  at 
the  spectra  (Figure  5-30).  This  suggests  that  the  Echelle  spectrometer  is  the  limiting 
detection  factor. 
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Figure  5-28:  Effect  of  interpulse  delay  on  intensity  on  the  588  nm  Na  peak  at  1  x 
10°  Pa.  (1000  ppm) 


Figure  5-29:  Na  calibration  curve  at  1  x  105  Pa  made  using  E]  =  20  mJ,  E2  =  140 
mJ,  and  t d  =  50  ns.  □  =  588  nm,  Q  =  589  nm 
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Figure  5-30:  Na  (I)  spectra  at  1  x  105  Pa.  (Ei  =  20  in J ,  E2  =  140  mJ,  and  =  50 
ns) 
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5.5  Conclusions 


Double  pulse  LIBS  was  used  to  detect  five  analytes  in  bulk  aqueous  solutions  impor¬ 
tant  in  hydrothermal  vent  chemistry  at  pressures  up  to  2.76  x  107  Pa.  The  key  double 
pulse  parameters  were  optimized  for  each  of  the  elements  at  three  pressures  (1  x  105 
Pa,  1.38  x  107  Pa,  and  2.76  x  10'  Pa).  The  parameters  needed  for  detection  were 
found  to  be  both  element  and  pressure  dependent.  The  use  of  the  optimal  parame¬ 
ters  is  essential  because  outside  of  this  set  of  parameters,  the  elements  may  not  be 
detectable.  Potassium  and  sodium  were  detectable  over  a  wide  range  of  conditions. 
In  general,  for  all  elements,  as  pressure  was  increased,  the  use  of  a  shorter  interpulse 
delay  was  necessary  and  at  2.76  x  107  Pa,  an  interpulse  delay  time  on  the  order  of 
50  ns  should  be  used.  For  all  conditions  studied,  a  short  gate  delay  (usually  <  100 
ns)  was  required.  This  is  similar  to  the  results  from  single  pulse  LIBS  (Chapter  4). 
The  need  for  a  short  gate  delay  suggests  that  in  bulk  liquids,  the  plasma  lifetime  is 
short,  possibly  lasting  only  on  the  order  of  500  ns. 

Using  each  of  the  optimally  established  conditions,  calibration  curves  were  made 
at  three  pressures  (1  x  10r>  Pa,  1.38  x  107  Pa,  and  2.76  x  107  Pa).  From  these,  the 
limits  of  detection  for  the  five  analytes  were  found  to  be  5000  ppm  Mg,  500  ppm  K, 
500  ppm  Ca,  1000  ppm  Mn,  and  50  ppm  Na  using  the  current  system  set-up.  The 
limits  of  detection  were  the  same  for  all  three  pressures  tested.  The  main  reason  that 
LIBS  researchers  choose  to  use  double  pulse  LIBS  instead  of  single  pulse  LIBS  is  to 
achieve  improved  sensitivity  and  improved  signal.  When  optimization  of  conditions 
was  completed  using  single  pulse  LIBS  and  then  subsequently  calibration  curves  were 
made  using  the  same  system  set-up  for  the  present  work  for  Mn,  Ca,  and  Na  (Chapter 
4),  the  resulting  limits  of  detection  were  found  to  be  500  ppm  Mn,  50  ppm  Ca,  and 
50  ppm  Na.  This  suggests  that  the  use  of  double  pulse  LIBS  in  high  pressure  aqueous 
solutions  may  not  be  advantageous.  Lawrence-Snyder  et  al.  reported  no  emission 
enhancements  using  double  pulse  LIBS  above  1  x  107  Pa  [16].  DeGiacomo  et  al. 
emphasize  the  need  for  selecting  an  appropriate  interpulse  delay  time  for  underwater 
LIBS  [23].  DeGiacomo  et  al.  stress  the  need  to  find  the  interpulse  delay  time  that 
allows  for  the  second  pulse  to  form  a  plasma  in  a  bubble  that  is  maximally  expanded. 
The  high  pressure  environment  of  our  experiments  may  cause  the  bubble  to  collapse 
too  rapidly.  As  a  result,  the  highest  intensity  peaks  that  are  observed  occur  when 
two  pulses  are  close  together,  similar  to  a  single  pulse.  Lawrence-Snyder  et  al.  [30) 
suggest  that  at  higher  solution  pressures  (8  x  107  Pa),  the  bubble  formed  by  the 
first  laser  pulse  is  confined  by  its  surrounding  pressure.  As  a  result,  the  bubble  never 
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expands  to  the  maximum  volume  that  is  observed  at  lower  pressures.  Therefore,  the 
emission  enhancements  are  not  observed. 

Although  the  use  of  double  pulse  LIBS  proved  less  favorable  than  expected,  it 
should  be  noted  that  one  major  contributing  factor  was  the  spectrometer  used  in 
these  studies.  The  Echelle  spectrometer  has  a  very  high  resolution  but  a  very  low  light 
throughput  and  poor  sensitivity,  with  an  /  number  of  10.  In  an  effort  to  maximize 
the  light  throughput,  it  would  be  advisable  to  use  a  spectrometer  with  a  smaller  / 
number.  As  an  example,  the  use  of  a  spectrometer  with  an  /  number  of  2  could 
improve  the  throughput  by  a  factor  of  25  and  thus  improve  the  limits  of  detection. 
Furthermore,  the  use  of  a  PMT  as  the  detector,  may  further  improve  detection  limits 
for  bulk  liquids  as  demonstrated  by  Cremers  et  al.  who  measured  at  atmospheric 
pressure  Na  I  (589.00  nm)  at  0.014  ppm,  K  I  (766.49  nm)  at  1.2  ppm,  Mg  II  (279.55 
nm)  at  100  ppm,  and  Ca  II  (393.37  nm)  at  0.8  ppm  [12].  Further  work  is  needed 
to  maximize  the  light  collection  by  changing  system  components  for  bulk  aqueous 
solution  experiments. 
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Chapter  6 


Preliminary  investigations  on 
matrix  effects  of  Na,  K,  and  Ca  for 
bulk  liquids  at  oceanic  pressures 

6.1  Abstract 

Chemical  matrix  effects  occur  when  one  element  present  in  a  sample  affects  the  emis¬ 
sion  of  another  element  also  present  in  the  sample.  In  this  study,  the  effect  of  Nad 
on  the  detection  of  K  and  Ca  in  bulk  aqueous  solutions  at  pressures  up  to  2.76  x 
107  Pa  is  explored.  In  addition,  the  effect  of  the  background  matrix  (chloride  versus 
sulfate)  on  the  detection  of  Na  and  K  is  examined.  While  the  investigations  into  the 
effect  of  NaCl  on  K  and  Ca  proved  inconclusive,  the  background  matrix  has  no  effect 
on  the  ability  to  detect  Na  or  K. 


6.2  Introduction 

When  multiple  elements  are  present  in  a  sample,  chemical  matrix  effects  can  oc¬ 
cur  in  which  the  presence  of  one  element  affects  the  emission  of  another  element. 
Changing  the  concentration  of  one  element  can  affect  the  signal  intensity  of  one  or 
more  other  elements  even  when  the  concentrations  of  those  elements  are  not  altered. 
This  can  adversely  affect  the  ability  to  make  quantitative  measurements  and  lead  to 
complications  in  calibration.  However,  if  matrix  effects  can  be  understood,  then  the 
concentration  effect  can  be  quantified  [1]. 

In  a  plasma,  the  addition  of  an  easily  ionizable  element  can  shift  the  neutral-ion 
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equilibrium  concentration  of  an  analyte  toward  the  neutral  species  (Lc  Chateliers 
principle).  For  example,  addition  of  an  easily  ionizable  species  to  a  plasma  can  reduce 
the  emission  of  an  ionized  species.  The  easily  ionized  species  increases  the  electron 
density  which  decreases  the  concentration  of  the  ionized  species  [1].  Crerners  et  al. 
[2]  reported  a  matrix  effect  in  bulk  liquids  which  showed  that  the  intensity  ratio  of  Ca 
Il/Ca  I  decreased  with  the  addition  of  NaCl.  A  similar  effect  was  reported  by  Michel 
et  al.  [3]  (Chapter  2)  using  double  pulse  LIBS  over  a  range  of  pressures.  Michel  et 
al.  also  reported  no  effect  on  K  or  Mil  emission  intensity  with  the  addition  of  Nad. 

Several  other  matrix  effects  have  been  reported.  For  example,  Charfi  and  Haritli 
[4]  showed  that  when  Mg  and  Na  were  measured  in  pure  solution  versus  in  mixed 
solutions  (both  Mg  and  Na  present),  the  limits  of  detection  were  lower  for  the  pure 
solutions.  Eppler  et  al.  [5]  investigated  matrix  effects  of  Ba  and  Pb  in  sand  and 
soil,  and  saw  that  Ba  (II)  and  Pb  (I)  emission  was  dependent  on  analyte  speciation. 
As  samples  were  varied  from  a  pure  sand  to  a  pure  soil  composition,  the  Ba  (II)/C 
(I)  signal  decreased.  The  decrease  correlated  with  an  increase  in  electron  density, 
possibly  due  to  a  change  in  ionized  species. 

Matrix  effects  can  be  corrected  by  applying  calibration  curves  for  each  element 
contained  in  the  substrate  of  interest,  analogous  to  estimating  the  partial  derivatives 
for  a  multivariate  process,  as  opposed  to  that  for  a  single  one  as  in  most  calibration 
curves.  This  requires  many  reference  samples  and  is  very  time  consuming.  In  addition, 
it  is  not  feasible  for  in  situ  measurements  with  unknown,  mixed  component  samples 
[0].  Several  research  groups  have  developed  methods  to  compensate  for  matrix  effects. 
St-Onge  et  al.  have  used  an  internal  standardization  method  to  reduce  matrix  effects 
[7].  Barrette  and  Turmel  showed  that  matrix  effects  could  be  partially  overcome  by 
the  use  of  a  multivariable  calibration  curve  [8].  Ciucci  et  al.  developed  an  algorithm- 
based  procedure  for  calibration-free  quantitative  elemental  analysis  of  materials  (9]. 

Understanding  matrix  effects  is  essential  for  quantitative  LIBS.  In  this  preliminary 
investigation,  matrix  effects  of  three  elements  (Na.  K,  and  Ca).  are  explored  in  bulk 
aqueous  solutions  at  pressures  up  to  2.76  x  10'  Pa.  In  addition,  the  effect  of  a 
background  matrix  (chloride  versus  sulfate)  on  the  detection  of  Na  and  K  is  examined. 


6.3  Experimental 

The  laboratory  set-up  for  simulating  the  use  of  LIBS  in  the  deep  ocean  has  been 
previously  described  in  Chapter  4.  For  all  experiments  described  here,  single  pulse 
LIBS  experiments  were  run  with  the  laser  pulse  energy  kept  constant  at  40  mJ  (which 
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corresponds  to  an  irradiance  of  the  beam  at  the  beam  waist  of  ss  5.25  x  1012  W/cm2 
for  this  system).  For  all  experiments,  t a  —  50  ns,  and  t <,  =  200  ns  and  the  Echelle 
spectrometer  amplification  was  set  to  the  maximum  value  of  4000.  All  spectra  taken 
were  composed  of  100  accumulated  shots.  All  raw  spectral  data  were  processed  using 
generalized  extreme  value  distribution  statistics  detailed  in  a  paper  by  Michel  and 
Chave  [10]  after  binning  the  data  over  20  wavelengths.  Where  shown,  error  bars 
represent  95%  confidence  intervals  of  the  GEVD. 

Solutions  were  made  from  NaCl.  KC1,  NaS04-H20,  MnS04-H20,  and  CaCl2-2H20 
dissolved  in  de-ionized  water.  All  concentrations  are  given  in  parts  per  million  (ppm. 
wt./vol.).  All  experiments  were  completed  at  three  pressures,  1  x  105  Pa,  1.38  x  107 
Pa,  and  2.76  x  107  Pa. 

6.4  Results  and  Discussion 

6.4.1  Matrix  Effects  of  K  and  Na 

Effect  of  Sodium  Chloride  on  the  Emission  from  Potassium 

Four  different  conditions  (Table  6.1)  were  used  to  study  the  effect  of  NaCl  on  the 
emission  intensity  of  the  769.9  nm  K  peak.  120  spectra  were  taken  at  each  of  Con¬ 
ditions  1  and  3  to  establish  a  low  and  high  calibration  point  for  K  alone.  Three 
replicate  measurements  each  consisting  of  60  spectra  were  taken  at  Conditions  2  and 
4  to  assess  the  effect  of  NaCl  on  K.  Figure  6-1  shows  the  results.  When  NaCl  was 
added  to  low  concentration  (102  ppm)  K,  no  peak  intensity  change  was  seen  at  all 
three  pressures.  In  contrast,  in  high  concentration  (1022  ppm)  K,  when  NaCl  was 
added,  the  K  peak  intensity  decreased  at  all  three  pressures.  This  is  in  contrast  to 
prior  work  [3]  (Chapter  2)  where  no  effect  of  NaCl  was  seen  on  the  peak  intensity  of 
K. 


Table  6.1:  Conditions  used  to  study  the  K-Na  matrix  effect 


Condition 

K 

Cl 

Na 

1 

102  ppm 

93  ppm 

0  ppm 

2 

102  ppm 

4959  ppm 

3155  ppm 

3 

1022  ppm 

927  ppm 

0  ppm 

4 

1022  ppm 

4959  ppm 

2614  ppm 
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Figure  6-1:  Effect  of  presence  of  Na  on  peak  intensity  of  K  (769.9  nm)  for  two 
concentrations  ofK.  0  =  K,IH  =  K  +  NaCl 
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Effect  of  Potassium  Chloride  on  Emission  from  Sodium 


Four  different  conditions  (Table  6.2)  were  used  to  study  the  effect  of  KC1  on  the 
detection  of  the  588  nm  Na  peak.  120  spectra  were  taken  of  Na  at  Conditions  1 
and  3  to  create  low  and  high  concentration  calibration  points  for  Na  alone.  Three 
replicate  measurements  each  consisting  of  60  spectra  were  taken  at  Conditions  2  and 
4  to  assess  the  effect  of  KC1  on  Na.  The  results  are  shown  in  Figure  6-2.  The  high 
variability  from  low  to  high  pressure  without  any  discernible  trend  suggests  that  more 
laboratory  work  is  needed  to  determine  if  a  matrix  effect  exists. 


Table  6.2: 


Conditions  used  to  study 


the  Na-K  matrix  effect 


Condition 

K 

Cl 

Na 

1 

0  ppm 

4865  ppm 

3155  ppm 

2 

102  ppm 

4959  ppm 

3155  ppm 

3 

0  ppm 

4032  ppm 

2614  ppm 

4 

1022  ppm 

4959  ppm 

2614  ppm 
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Figure  6-2:  Effect  of  presence  of  K  on  peak  intensity  of  Na  (588  nm).  □  -  Na,  O  - 
Na  +102  ppm  K,  A  =  Na  +  1022  ppm  K 
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6.4.2  Matrix  Effects  of  Ca  and  Na 


Effect  of  Sodium  Chloride  on  Emission  from  Calcium 

Four  different  conditions  (Table  6.3)  were  used  study  the  effect  of  NaCl  on  the  emission 
intensity  of  Ca  [393  Ca  (II)  and  422  nm  Ca  (I)].  120  spectra  were  taken  at  Conditions 
1  and  3  to  establish  low  and  high  concentration  calibration  points  for  Ca  alone.  Three 
replicate  measurements  each  consisting  of  60  spectra  were  taken  at  Conditions  2  and 
4  to  assess  the  effect  of  NaCl  on  the  393  nm  Ca  (II)  ionic  peak  and  the  422  nm  (I) 
atomic  peak.  Figure  6-3  shows  the  effect  of  NaCl  on  the  393  nm  Ca  (II)  ionic  peak. 
For  106  ppm  Ca,  the  peak  intensity  was  not  measurable  at  any  pressure,  both  in  the 
presence  and  absence  of  NaCl.  When  NaCl  is  added  to  a  1055  ppm  Ca  solution,  the 
emission  intensity  is  significantly  lower.  It  should  be  noted  that  at  1  x  105  Pa  and 
1.38  x  107  Pa,  the  measured  intensity  under  all  conditions  is  relatively  low.  Figure 
6-4  shows  the  results  for  the  effect  of  NaCl  on  the  422  nm  Ca  (I)  atomic  peak.  The 
results  are  inconclusive.  To  assess  possible  drift  in  the  spectrometer  between  the  pure 
Ca  data  and  the  Ca  plus  NaCl  data,  the  ratio  of  the  intensity  of  the  393  nm  peak  to 
the  422  nm  peak  is  plotted  in  Figure  6-5.  The  result  suggests  that  the  variability  is 
not  due  to  simple  drift,  but  to  another  cause  of  systematic  error. 

Table  6.3:  Conditions  used  to  study  the  Ca-Na  Matrix  Effect 


Condition 

Ca 

Cl 

Na 

1 

106  ppm 

187  ppm 

0  ppm 

2 

106  ppm 

4959  ppm 

3094  ppm 

3 

1055  ppm 

1867  ppm 

0  ppm 

4 

1055  ppm 

4959  ppm 

2005  ppm 
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Figure  6-3:  Effect  of  presence  of  Na  on  peak  intensity  of  Ca  (II)  (393  nm)  for  two 
concentrations  of  Ca.  Q  =  Ca,  □  =  Ca  +  NaCl 
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Figure  6-4:  Effect  of  presence  of  Na  on  peak  intensity  on  Ca  (I)  (422  nm)  for  two 
concentrations  of  Ca.  Q  =  Ca,  □  =  Ca  +  NaCl 
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Figure  6-5:  Ratio  of  393  inn  Ca  (II)  peak  to  422  nm  Ca  (I)  peak  Q  =  Ca,  □  =  Ca 
+  NaCl 
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Effect  of  Calcium  Chloride  on  the  Emission  of  Sodium 

A  sodium  calibration  curve  was  made  with  120  spectra  taken  at  three  concentrations 
(Conditions  1-3  shown  in  Table  6.4).  Three  replicates  of  60  spectra  were  taken 
at  two  concentrations  of  Ca  plus  NaCl.  The  data  shown  in  Figure  6-6  reveal  high 
variability  in  the  LIBS  system  measurements.  In  general,  when  Ca  was  present,  the 
Na  peak  intensity  was  reduced.  However,  the  systematic  error  denoted  by  the  scatter 
in  the  data  suggests  that  more  data  collection  is  needed  to  confirm  the  result . 


Table  6.4: 


Conditions  used  to  study  the  effect  of  CaC 


•2H20  on  emission  of  Na 


Condition 

Ca 

Cl 

Na 

1 

0  ppm 

3091  ppm 

2005  ppm 

2 

0  ppm 

4032  ppm 

2614  ppm 

3 

0  ppm 

4865  ppm 

3155  ppm 

4 

1055  ppm 

4959  ppm 

2005  ppm 

5 

106  ppm 

4959  ppm 

3094  ppm 
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Figure  6-6:  Effect  of  presence  of  Ca  on  peak  intensity  on  Na  (588  urn),  □  =  Na,  A 
=  Na  +  1055  ppm  Ca,  Q  =  106  ppm  Ca  +  Na 
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6.4.3  Detection  of  Na  and  K  in  a  Chloride  Versus  Sulfate 
Matrix 

Two  experiments  were  carried  out  to  determine  if  the  presence  of  a  chloride  or  sulfate 
substrate  affects  the  detection  of  an  element.  The  intensity  of  potassium  (1002  ppm 
K)  in  two  matrices,  KC1  and  K2S04,  and  the  intensity  of  sodium  (2356  ppm  Na), 
in  two  matrices,  NaCl  and  NaSO.j,  were  measured.  Figures  6-7  and  6-8  contrast  the 
ability  to  detect  K  and  Na,  respectively,  in  a  chloride  matrix  versus  a  sulfate  matrix. 
For  K,  no  matrix  effect  is  seen  at  any  of  the  pressures  (lx  10;>  Pa,  1.38  x  107  Pa, 
2.76  x  107  Pa)  tested.  For  Na  (Figure  6-8),  there  is  a  possible  matrix  effect,  although 
additional  experiments  are  necessary  to  determine  this  due  to  scatter  in  the  data. 
In  the  development  of  an  oceanic  LIBS  sensor  operable  in  the  vent  environment,  the 
sensors  must  be  capable  of  detecting  both  K  and  Na  in  the  presence  of  both  sulfate 
and  chloride,  which  these  data  support. 


(a) 


Figure  6-7:  Detection  of  769  nm  potassium  (1002  ppm)  in  a  chloride  versus  sulfate 
matrix.  Q  =  KC1,  A  =  K2S04 
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Figure  6-8:  Detection  of  sodium  (2356  ppm  Na)  in  a  chloride  versus  sulfate  matrix. 
O  =  NaCl.  A  =  NaSQ4 
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6.5  Conclusions 


Matrix  effects  of  Na,  K,  and  Ca  were  examined  in  a  preliminary  study.  Although 
the  main  goal  was  determining  how  the  presence  of  NaCl  affects  the  peak  intensities 
of  K  and  Ca,  the  data  showed  inconsistencies  and  were  therefore  inconclusive.  In 
addition,  the  results  were  not  consistent  with  previous  preliminary  findings  reported 
by  Cremers  et  al.  [1]  and  Michel  et  al.  [3].  The  data  suggest  that  matrix  effects  are 
smaller  than  the  systematic  variability  of  the  measurements,  and  suggest  that  further 
laboratory  experimentation  is  needed  to  elucidate  if  significant  matrix  effects  exist. 

A  comparative  study  to  investigate  the  ability  to  detect  analytes  in  a  chloride 
versus  a  sulfate  matrix  showed  no  significant  difference.  The  results  indicate  that  Na 
and  K  are  both  detectable  in  chloride  and  sulfate  matrices.  The  K  peak  intensities 
were  not  influenced  by  the  substrate  present  (sulfate  and  chloride)  at  pressures  up  to 
2.76  x  107  Pa.  The  Na  peak  intensities  were  possibly  influenced  by  the  background 
matrix  (sulfate  and  chloride);  yet,  more  experiments  are  needed  to  determine  this 
conclusively. 
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Chapter  7 


Conclusions  and  Future  Directions 


7.1  Conclusions 

In  this  thesis,  laser-induced  breakdown  spectroscopy  was  evaluated  for  its  potential 
for  development  into  a  new  in  situ  chemical  sensor  for  the  deep  ocean.  The  elements  of 
focus  for  this  work  were  selected  for  their  importance  in  hydrothermal  vent  chemistry. 
Although  it  is  potentially  and  theoretically  possible  to  use  LIBS  to  detect  all  elements, 
it  was  found  during  this  laboratory  investigation  that  it  is  more  difficult  to  detect 
elements  in  bulk  liquids  at  the  concentrations  desired  than  was  anticipated.  However, 
one  important  finding  was  that  if  an  element  could  be  detected  in  a  bulk  liquid 
at  atmospheric  pressure,  then  it  could  also  be  detected  in  the  same  liquid  at  high 
pressures.  This  thesis  shows  that  the  LIBS  technique  can  be  used  successfully  to 
detect  dissolved  Li,  Na,  K,  Mg,  Ca,  and  Mn  at  pressures  up  to  2.76  x  107  Pa, 
although  not  at  the  levels  observed  naturally  in  seawater  or  vent  fluids  for  all  of  these 
elements.  The  work  in  this  thesis  was  completed  using  two  different  LIBS  systems, 
one  with  a  Czerny-Turner  spectrometer  and  one  with  an  Echelle  spcctometcr,  and 
with  two  different  pressure  cells. 

7.1.1  Development  of  a  New  Data  Processing  Scheme 

A  new  data  processing  scheme  for  LIBS  spectra  was  developed  in  this  thesis  and  was 
applied  to  all  data,  except  that  presented  in  Chapter  2.  In  an  examination  of  the 
variability  of  peak  intensities  for  both  single  shot  and  ensemble-averaged  LIBS  spec¬ 
tra,  LIBS  data  were  found  to  have  a  dramatically  non-normal  statistical  distribution. 
The  distribution  of  the  peak  intensities  was  instead  found  to  follow  the  generalized 
extreme  value  distribution.  A  preliminary  investigation  into  the  sources  of  the  vari- 
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ability  was  carried  out.  Laser  pulse  fluctuations,  while  identified  as  a  contributing 
source  of  variability,  were  ruled  out  as  the  the  primary  source.  Plasma  images  re¬ 
vealed  large  spatial  and  intensity  differences  on  a  shot-to-shot  basis.  This  analysis 
led  to  the  development  of  a  data  processing  scheme  that  accurately  deals  with  the 
extreme  nature  of  laser-induced  plasma  formation  and  should  be  used  for  statistically 
accurate  comparisons  of  LIBS  spectra  instead  of  simply  averaging  spectra,  the  stan¬ 
dard  method  used  by  LIBS  researchers.  This  scheme  was  found  to  be  applicable  for 
both  solid  and  liquid  samples  and  for  use  with  data  taken  with  both  a  Czerny-Turner 
and  an  Echelle  spectrometer. 

7.1.2  Single  Pulse  LIBS 

Results  Using  a  Czerny-Turner  Spectrometer 

Preliminary  investigations  wore  carried  out  using  a  Czerny-Turner  spectrometer  that 
showed  that  Li.  Ca,  Mn,  K,  and  Na  were  detectable  in  high  pressure  bulk  aqueous 
solutions.  This  work  focused  on  the  energy  levels  needed  for  the  detection  of  these 
analytes  and  revealed  that  an  optimal  range  of  low  laser  pulse  energies  exists  for  their 
detection  in  both  low  and  high  pressure  solutions.  It  was  hypothesized  that  a  low- 
energy  pulse  could  create  a  smaller,  more  tightly  focused  plasma  that  forms  only  at 
the  focal  spot.  However,  for  a  high  energy  pulse,  the  high  energy  density  may  cause 
breakdown  even  before  the  pulse  reaches  the  focal  spot  causing  breakdown  to  occur 
over  a  longer  distance.  The  effect  of  pressure  on  emission  intensity  was  investigated 
for  Ca,  Na,  and  M11  with  no  pressure  effect  seen  for  Ca  and  Na,  yet.  with  an  increase 
in  intensity  with  increased  pressure  observed  for  Mn. 


Results  Using  an  Echelle  Spectrometer 

Using  an  Echelle  spectrometer,  Na,  Mn,  and  Ca  were  again  shown  to  be  detectable 
with  little  effect  of  pressure  on  the  spectra.  Irrespective  of  the  laser  pulse  energy 
selected  or  the  solution  pressure,  the  need  for  a  short  (less  than  200  11s)  gate  delay 
was  found  for  the  detection  of  analytes.  This  study  confirmed  that  a  low  energy  pulse 
(less  than  «  GO  mJ)  is  optimal.  Calibration  curves  showed  that  Na,  M11,  and  Ca 
are  detectable  at  50  ppm,  500  ppm,  and  50  ppm,  respectively.  However,  the  limits  of 
detection  were  higher  than  expected  which  was  attributed  to  both  the  poor  sensitivity 
and  the  low  light  throughput  of  the  system  used. 
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7.1.3  Double  Pulse  LIBS 


Results  Using  a  Czerny-Turner  Spectrometer 

In  this  preliminary  work,  double  pulse  LIBS  was  used  to  detect  four  analytes  (Ca, 
Li.  Na,  and  Mn)  and  the  energy  levels  for  maximum  emission  were  found  to  vary 
by  analyte.  Analyte  detection  was  found  to  be  highly  dependent  on  the  interpulse 
delay.  If  the  interpulse  delay  is  short  (<^C  1  /xs),  signal  intensity  is  greatly  enhanced 
when  compared  to  that  measured  using  longer  delay  times.  Two  pulses  separated  by 
a  short  AT  approaches  single  pulse  conditions.  Such  a  small  inter  pulse  delay  may 
not  be  sufficient  for  a  cavitation  bubble  to  fully  form  before  the  second  laser  pulse 
creates  a  spark.  This  was  an  early  suggestion  that  dual  pulse  LIBS  might  not  be 
advantageous  at  elevated  pressure. 

Results  Using  an  Echelle  Spectrometer 

Double  pulse  LIBS  was  used  to  detect  five  analytes  (Mg,  K,  Mn,  Na,  and  Ca)  at 
pressures  up  to  2.76  x  107  Pa.  The  key  double  pulse  parameters  were  optimized 
for  each  of  the  elements  at  three  pressures  (1  x  10°  Pa,  1.38  x  10'  Pa,  and  2.76  x 
107  Pa).  The  parameters  needed  for  detection  were  found  to  be  both  element  and 
pressure  dependent.  The  use  of  the  optimal  parameters  is  essential  because  outside 
this  set  of  parameters,  some  of  the  elements  were  not  detectable,  although  K  and 
Na  were  detectable  over  a  wide  range  of  conditions.  In  general,  for  all  elements,  as 
pressure  was  increased,  the  use  of  a  shorter  interpulse  delay  was  necessary  and  at  2.76 
x  107  Pa  an  interpulse  delay  time  on  the  order  of  50  ns  was  ideal.  For  all  conditions 
studied,  a  short  gate  delay  (usually  <  100  ns)  was  required,  suggesting  that  in  bulk 
liquids,  the  plasma  lifetime  is  short,  possibly  lasting  only  on  the  order  of  500  ns. 

Calibration  curves  were  made  that  established  the  limits  of  detection  for  the  five 
analytes:  5000  ppm  Mg,  500  ppm  K,  500  ppm  Ca,  1000  ppm  Mn,  and  50  ppm  Na 
using  the  current  system  set-up.  Using  single  pulse  LIBS,  limits  of  detection  were 
found  to  be  500  ppm  Mn,  50  ppm  Ca,  and  50  ppm  Na  which  suggests  that  the  use 
of  double  pulse  LIBS  in  high  pressure  aqueous  solutions  may  not  be  advantageous. 
The  high  pressure  environment  may  cause  the  bubble  to  collapse  too  rapidly  and 
as  a  result  the  bubble  does  not  expand  to  the  maximum  volume  observed  at  lower 
pressures  [1].  As  a  result,  the  emission  enhancements  seen  at  atmospheric  pressure 
are  not  observed  in  the  high  pressure  environment. 

Although  the  use  of  double  pulse  LIBS  proved  less  favorable  than  expected,  it 
should  be  noted  that  one  major  contributing  factor  was  the  spectrometer  used  in 
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these  studies.  The  Echelle  spectrometer  has  a  very  high  resolution  but  a  very  low 
light  throughput  and  poor  sensitivity,  with  an  /  number  of  10.  Further  work  is 
needed  to  maximize  the  light  collection  by  changing  system  components  for  bulk 
aqueous  solution  experiments. 

7.1.4  Matrix 

Results  Using  a  Czerny-Turner  Spectrometer 

The  addition  of  NaCl  was  found  to  enhance  the  emission  intensity  for  Ca. 

Results  Using  an  Echelle  Spectrometer 

The  matrix  effect  of  Na,  K.  and  Ca  was  examined  to  determine  how  the  presence  of 
NaCl  affects  the  peak  intensity  of  K  and  Ca,  although  the  data  showed  inconsistencies 
and  were  therefore  inconclusive.  Na  and  K  were  found  to  be  detectable  in  both 
chloride  and  sulfate  matrices.  For  Na  the  peak  intensities  were  possibly  influenced 
and  for  K  the  peak  intensities  were  not  influenced  by  the  background  matrix  (sulfate 
versus  chloride)  at  pressures  up  to  2.76  x  107  Pa. 

7.2  Future  Work 

7.2.1  Laboratory  Work 

Significantly  more  laboratory  work  is  needed  before  a  robust  ocean-going  system  can 
be  built  and  used  successfully  in  the  field.  Although  this  thesis  showed  that  many 
elements  could  be  detected  in  a  high  pressure  aqueous  environment .  and  several  were 
detectable  at  levels  found  at  or  below  the  concentrations  present  in  vent  fluids,  there 
are  many  issues  that  need  to  be  addressed  to  improve  upon  the  current  LIBS  system. 

Due  to  the  attenuation  of  light  in  water,  532  nm  or  green  laser  light  will  travel 
substantially  further  than  1064  nm  laser  light  in  the  ocean.  Therefore,  it  will  be 
advantageous  to  develop  a  LIBS  system  that  uses  532  nm  laser  pulses  as  the  excitation 
source.  The  disadvantage  to  using  532  nm  laser  pulses  as  the  excitation  source  is  the 
need  for  a  filter  before  the  collection  fiber  optic.  Since  532  nm  is  in  the  range  that 
can  be  detected  by  the  spectrometer,  it  is  important  to  never  collect  the  532  nm  laser 
pulse  light.  To  prevent  capturing  the  laser  light,  a  532  nm  notch  filter  is  added  before 
the  collection  fiber  optic  to  filter  out  the  light  from  the  laser  pulses.  The  addition 
of  this  filter  will  minimally  reduce  the  light  available  to  the  spectrometer.  LIBS  of 
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liquids  using  532  nni  should  be  feasible;  therefore,  the  next  phase  of  laboratory  work 
should  explore  this  option. 

Improvements  to  the  limits  of  detection  could  be  made  by  using  a  spectrometer 
with  a  smaller  /  number  than  the  Echelle  spectrometer  that  was  used  for  much  of 
this  thesis  work.  Although  an  Echelle  spectrometer  has  the  advantage  of  being  able 
to  simultaneously  detect  multiple  elements  due  to  its  broad  wavelength  coverage,  its 
high  resolution  yields  low  throughput.  The  Echelle  spectrometer  also  had  a  small  slit 
width  which  limited  light  input.  In  addition,  the  use  of  a  PMT  instead  of  an  ICCD 
as  the  detector  could  be  beneficial  for  increasing  sensitivity.  Changes  therefore  to  the 
detection  system  may  allow  for  the  elements  presented  in  this  thesis  to  be  detected 
at  much  lower  limits  of  detection  and  for  additional  elements  to  be  detected.  For 
example,  the  detection  of  Cu,  Fe,  Cl,  and  Br  was  attempted  as  part  of  this  research 
without  success.  It  is  recommended  that  a  new  spectrometer  with  a  smaller  /  number 
be  selected  for  continued  work  on  this  project.  The  Echelle  spectrometer  used  in  this 
work  has  an  /  number  of  10  and  it  is  strongly  recommended  to  use  a  spectrometer 
with  an  /  number  of  4  or  2,  which  would  have  greater  throughput. 

Future  work  should  include  significant  studies  on  the  variability  of  plasma  forma¬ 
tion.  This  can  include  imaging  of  both  laser-induced  plasmas  and  bubbles.  Plasma 
dynamics  studies  should  be  carried  out  in  an  effort  to  find  ways  to  minimize  shot- 
to-shot  plasma  variability.  Although  this  thesis  included  some  plasma  images  taken 
orthogonal  from  the  incoming  laser  beam,  plasma  images  could  also  be  taken  from 
multiple  angles  to  gain  a  multi-dimensional  understanding  of  plasma  formation.  Then 
by  looking  at  what  factors  influence  plasma  dynamics  (e.g.,  the  shape  of  the  plasma, 
the  length  of  emission,  emission  intensity),  more  control  over  the  plasma  variability 
may  be  obtainable.  The  factors  that  influence  the  plasma  dynamics  include  the  laser 
(e.g.,  beam  quality,  wavelength),  focusing  optics,  and  sample  pressure  and  should  also 
be  studied  for  their  impact  on  variability. 

Calibration  curves  were  used  in  this  thesis  to  determine  the  limits  of  detection. 
In  future  work,  much  work  should  be  completed  on  the  calibration  curves  to  both 
decrease  the  variability  of  the  curves  and  to  verify  repeatability  of  the  curves.  The 
calibration  curves  should  also  be  extended  to  higher  concentrations  to  verify  that 
the  curves  are  linear.  In  addition,  if  a  new  spectrometer  is  used,  calibration  curves 
can  again  be  used  to  determine  limits  of  detection.  Once  the  variability  issues  are 
resolved  and  the  repeatability  of  calibration  curves  are  established,  more  work  can 
focus  on  matrix  effects.  Matrix  effect  data  showed  inconsistencies  that  are  thought 
to  arise  from  the  variability  issues  and  from  the  spectrometer  used.  Selecting  a  new 
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spectrometer  may  improve  these  studies;  however,  the  reduction  of  plasma  variability 
again  should  first  be  addressed. 

There  are  many  additional  challenges  to  studying  liquids  in  the  ocean  environment 
that  were  not  addressed  in  this  thesis.  For  example,  the  effect  of  particles  in  the 
water  on  the  detection  of  analytes  should  be  addressed.  In  addition,  the  effect  of 
a  flowing  solution  on  LIBS  spectra  and  especially  the  shot-to-shot  variability  in  a 
flowing  system  should  be  examined.  Furthermore,  laboratory  work  should  emphasize 
elements  not  yet  targeted.  In  theory,  all  elements  are  detectable  using  LIBS;  thus, 
significant  work  is  needed  to  determine  all  of  the  elements  which  are  detectable  in 
high  pressure  aqueous  environments. 


7.2.2  Design  of  an  Ocean-Going  Sensor 

Although  there  are  several  issues  that  need  to  be  addressed  in  the  laboratory  to  deal 
with  the  variability  issues  and  to  improve  upon  the  calibration  curves  and  limits  of 
detection,  this  thesis  shows  that  the  use  of  LIBS  in  a  high  pressure  aqueous  envi¬ 
ronment  is  feasible.  Once  these  issues  have  been  addressed,  a  LIBS  sensor  can  be 
designed  for  use  in  the  ocean.  The  natural  extension  of  this  research  therefore  is 
the  development  of  the  first  ever  sea-going  LIBS  system.  This  will  entail  the  de¬ 
sign  of  an  instrument  that  can  be  used  at  hydrothermal  vents  aboard  an  underwater 
vehicle  (Figure  7-1).  The  major  components  of  an  underwater  LIBS  system  are  a 
spectrometer,  source  laser(s),  fiber  optic  links  and  probes  for  these  units,  and  data 
acquisition/control  electronics.  All  of  these  components  will  need  to  be  packaged 
in  underwater  pressure  housings  and  provided  with  power  and  data  communications 
connections.  Several  critical  issues  related  to  integrating  the  system  onto  an  under¬ 
water  vehicle  also  need  to  be  addressed,  including  vehicle  payload,  power  resources, 
and  system  control.  Software  modifications  for  laser  control  and  data  collection  will 
be  a  necessary  component  of  this  work.  The  design  of  an  optical  probe  head  is  a 
critical  element  of  a  sea-going  LIBS  system.  The  fiber  optic  design  must  pass  the 
laser  pulses  into  the  fiber,  image  the  fiber  spot  onto  the  target,  and  image  the  plasma 
emission  onto  a  fiber  bundle  that  couples  to  the  spectrometer.  The  fiber  bundle  must 
be  designed  to  withstand  both  the  pressure  effects  of  the  ocean  and  also  the  high 
energy  laser  pulses.  In  the  ocean  environment,  there  are  other  types  of  samples  that 
ocean  scientists  will  be  interested  in  using  an  in  situ  chemical  LIBS  sensor  for  mak¬ 
ing  measurements  of  in  the  ocean.  For  example,  LIBS  could  be  used  for  analyzing 
sediments,  rocks,  underwater  archaeological  artifacts,  and  the  air-sea  interface.  The 
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successful  design  and  testing  of  an  ocean-going  system  will  pave  the  way  for  its  use 
by  oceanographers  for  numerous  applications. 
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Figure  7-1:  Future  oceanic  LIBS  system.  Illustration  by  E.  Paul  Obcrlander,  WHOI. 
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